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Chapter 1 

Introduction 



Quantum dots have recently attracted experimental interest due to their 'possible' ap- 
plication in future nano devices and for quantum information processing. From the 
theoretical point of view a many-particle method suited to account for the two-particle 
interactions between the electrons in the dot - which strongly affect the physics - is 
needed. In the literature, the numerical renormalization group (NRG) was frequently 
used as a very reliable method to study systems with local Coulomb correlations, but 
its applicability is limited due to the vast computational resources required. Practically, 
only simple geometries can be treated and even there exhaustive scans of the param- 
eter space are impossible. Therefore other accurate methods that allow for a precise 
treatment of the two-particle interactions in quantum dots are needed, but none of the 
recent approaches succeeded in reaching the accuracy of the NRG calculations, most 
times not even on a qualitative level [KEM 2006]. 

Here, we will present a recently-developed renormalization group (RG) scheme based 
on Wilson's general RG idea, the functional renormalization group (fRG). Its starting 
point is the replacement of the free propagator by one depending on an infrared cutoff 
A in the generating functional of the single-particle irreducible vertex functions 7^. 
Differentiating the functional with respect to A one obtains an exact hierarchy of cou- 
pled flow equations for the functions 7^. This infinite set has to be truncated in order 
to render it solvable, and within this thesis we will employ a truncation scheme that 
accounts for the flow of 71 (the self-energy) and the two-particle vertex 72 evaluated 
at zero external frequency (the effectice interaction), but neglects all higher-order func- 
tions. Hence we compute a frequency-independent approximation for 71 which can be 
viewed as a renormalization of the dot's single-particle energies. 

An important transport property of a quantum dot is the linear-response conductance 
G. In a noninteracting model, it is frequently defined using single-particle scattering 
theory, which leads to the well-known Landauer-Biittiker formula. Since the afore- 
mentioned fRG approximation computes effective noninteracting parameters for the 
quantum dot system, one might be tempted to use the ordinary [7 = expression 
to calculate G. This is, however, conceptionally wrong since it neglects the fact that 
this quantity cannot be defined a priori using the ordinary scattering theory approach 
as one is confronted with a system of correlated electrons, whether the approximation 
we employ to describe it can be interpreted as an effective noninteracting one or not. 
Hence a better way to define the conductance is to consider the current-current response 
function within the full interacting problem and to show that the noninteracting ex- 
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pression follows consistently within our fRG framework without need to apply further 
approximations . 

Having an accurate method at hand to describe Coulomb correlations, we will first 
consider the low-energy physics (at zero temperature) of a spin-polarised dot containing 
several levels. This can be viewed as a suitable model for an experimentally realised sin- 
gle quantum dot, assuming that the experiment showed that the spin degree of freedom 
seems to play no role. A series of important experiments on such systems was performed 
in the past few years [YHMS 1995, SBHMUS 1997,AHZMU 2005], and in addition to 
the conductance also the transmission phase and, in the very recent one, the average 
occupancy of the dot were measured. The conductance was found to exhibit resonances 
(Coulomb blockade peaks) each time another electron was added to the system. Surpris- 
ingly, the measurement of the phase yielded universal jumps by vr between consecutive 
peaks, in contradiction to the expected mesoscopic behaviour of jumping between some 
while evolving continuously between others, the actual realisation depending on the dot 
under consideration. Despite a large amount of theoretical work on this issue (for a 
review see [Hackenbroich 2001,Gefen 2002]; the most important works with electron- 
electron correlations as a key ingredient are [BRS 1995, Silvestrov&Imry 2000, SSOD 
2005, Konig&Gefen 2005, Golosov&Gefen 2006]), none succeeded in giving a fully satis- 
factory explanation. Fortunately, the last experiment ( [AHZMU 2005]) provided a clue 
to this phase lapse puzzle. For dots occupied by only a few electrons, the mesoscopic 
behaviour was indeed observed, and only if electrons were successively added to the 
system, the universally appearing phase lapses were recovered. As noted by [AHZMU 
2005], the most fundamental difference between both situations is the single-particle 
level spacing which one would expect to decrease for the electrons occupying the top- 
most states. Exploring this idea in a systematic study of the parameter space of our 
spinless quantum dot model with an accurate treatment of Coulomb correlations is the 
aim of the first part of this thesis. 

Next, we describe a variety of spinful dot geometries where the physics is dominated 
by local correlations between the electrons of both spin directions. We will show how the 
Kondo energy scale can be extracted from the linear conductance in the zero tempera- 
ture limit. We will study two geometries (a short Hubbard chain and two side-coupled 
dots) that have been tackled before by various authors using NRG, and we will demon- 
strate the power of the fRG by scanning a much larger (and in particular generic) region 
of the parameter space of these systems. For the Hubbard chain this will reveal new 
physics, while for the side-coupled geometry it will turn out that the behaviour observed 
within the NRG approach already covers all essential phenomena. We will finally in- 
vestigate the situation of two parallel dots. This geometry allows for a large amount of 
parameters to be varied, and we will report on the physics that is observed in different 
regimes. 

In the last chapter, we will establish the accuracy of the employed fRG approximation 
scheme by comparison to NRG data for every system that has previously been under 
consideration. We will show that the fRG is very reliable up to interaction strengths 
where the physics is dominated by correlations (in the sense that the behaviour of all 
quantities of interest is significantly altered with respect to the noninteracting case). 
Furthermore, we will prove that this method is far superior to simple mean-field ap- 
proaches. 
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The Functional Renormalization 
Group 

The functional renormalization group (fRG) is based on Wilson's general renormaliza- 
tion group (RG) idea. The RG is the appropriate framework to deal with quantum 
many-particle problems that contain a multitude of different energy scales as well as 
with problems that cannot be tackled by ordinary perturbation theory which breaks 
down if certain classes of (low-order) Feynman diagrams diverge. It starts from high 
energy scales (where infrared divergences are cut out) and gradually works its way down 
to the low energy region. 

The fRG provides one way to implement the general renormalization group idea for in- 
teracting quantum many-particle systems [Salmhofer 1998]. It is based on the functional 
integral approach to many-particle physics. Differentiating the generating functional of 
the single-particle irreducible m-particle vertex functions 7m (with the free propagator 
replaced by one that cuts out the low energy modes below a scale A) with respect to 
a cutoff parameter yields an infinite hierarchy of flow equations for the jm- In prac- 
tise, this set has to be truncated to make it solvable. A truncation scheme motivated 
by physical or practical reasons renders the fRG an approximate method, most times 
perturbative in the two-particle interaction. Integrating the resulting finite set of equa- 
tions it turns out that in the limit A — > 0, which corresponds to the original cutoff-free 
problem, one often obtains non-divergent results, and even in problems that are not 
plagued by infrared divergences (as those tackled within this thesis) it shows that the 
fRG produces results far superior to those of perturbation theory. An application of 
the method to the well-known one-dimensional single-impurity Anderson model can be 
found in [HMPS 2004]. In [EMABMS 2005] the fRG is used to successfully derive power 
laws in correlated one-dimensional quantum wires (Luttinger liquids). 

This chapter is organised as follows. In the first section, we recall the functional inte- 
gral approach to many-particle physics and how it is used to compute Green functions. 
In the second section, we replace the noninteracting propagator in the generating func- 
tional of the vertex functions by one containing a cutoff. Differentiating with respect to 
the cutoff parameter A then yields an infinite hierarchy of flow equations for the vertex 
functions. Since this hierarchy needs to be truncated in order to make it solvable, we 
will develop a truncation scheme that becomes exact in the limit of vanishing strength 
of the two-particle interaction. Finally, we will specify the cutoff to be an infrared cutoff 
in frequency space to write down coupled flow equations that perturbatively describe 
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an interacting system. 



2.1 Many-particle Green Functions 

In this section we recall how functional integrals can be used to compute the grand 
canonical partition function of a system of interacting quantum mechanical particles as 
well as Green functions. This will be the appropriate framework to set up the fRG flow 
equations. 

2.1.1 The Partition Function as a Functional Integral 

A Short Reminder of Functional Integrals in Many-Particle Physics 

The fundamental equation for the functional integral approach to many-particle physics 
reads 

e-P"= lim f:e-'3^/^:')'^, 

where is a normal-ordered but otherwise arbitrary many-particle Hamiltonian, and 
the colons denote normal-ordering. We insert the unity operator written in terms of 
the eigenstates of the annihilation operators between each of two factors of : e~^^l^ :. 
Since we are dealing with fermions, we have to introduce Grassmann variables, that is 
anti-commuting numbers, to write down these states. Because of the normal-ordering, 
the annihilation (creation) operators in the exponential can then be replaced by the 
corresponding eigenvalue and we are basically left with an 'ordinary' (Grassmann) in- 
tegral, that becomes a functional integral in the limit N ^ oo. In case that H contains 
only single-particle terms, this integral is purely Gaussian and can be performed 'almost 
trivially'. If H contains more complicated terms (such as a two-particle interaction), 
carrying out the the integral analytically becomes impossible, but one can expand the 
interaction part of the exponential function into a Taylor series (if one notes that we 
are now only dealing with numbers and no longer with operators) to set up a perturba- 
tion theory. The contribution of each order can be evaluated using Wick's theorem, a 
property of Gaussian integrals. 

A more detailed description of the functional integral approach to many-particle 
physics can be found in [Negele&Orland 1988], or, very nicely, in [Schonhammer 2001]. 

The Partition Function 

As motivated above, the grand canonical partition function of a fermionic many-particle 
system can be expressed as a functional integral. 



fl3 

Z = I Vipip exp { — I dr 







where ijj and ijj denote Grassman variables. 



(2.1) 



(-) (-) 
■0; V'z' + V'r ill = 0, 
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and TL is obtained from a normal-ordered (but otherwise arbitrary) Hamilton operator 
H by substituting 

n ({^} , {^}) =h[c\^ Ur + 0), Q ^ Mr)) ■ 
The boundary conditions of the fields read 

UP) = -^/(o), W) = -m, 

with P being the inverse temperature. In the following, we will always choose 7i to 
contain a single-particle term as well as a two-particle interaction, 

ni{tl;} ,{^P}) = ^eiMT + 0)Mr) + \ J2 ^wMr + 0)^^ + 0)Mr)Mr). (2.2) 

I 

Here I denotes a set of quantum numbers in which the single-particle part of H is 
diagonal, ei the one-particle dispersion (including the chemical potential) and v the 
anti-symmetrised matrix elements of the two-particle interaction, respectively. 

For the upcoming calculations it will prove useful to expand the Grassmann fields 
into a Fourier series 

^^^^) = ;^E^'^'"" W(^^n). (2.3) 

The summation extends over all odd (fermionic) Matsubara frequencies. The inverse 
transformation is given by 



(-) I fP ,. (-) 

il^i (ia;„) = — dr e±*""^ V; (r). (2.4) 



VP 

Inserting (|2.3I) into the imaginary-time functional integral (|2.ip yields 



Z 

Zo Zt 



' I ^ ^ ' (2.5) 

4"^^ '^ij,k,l5m+m',n+n''4^i{i^rn)'4^j{i^^m')Ui^ 



with Q'^iiijjn) being the noninteracting single-particle propagator.^ All convergence fac- 
tors have been omitted to keep the notation short and the noninteracting partition 
function Zq was introduced to cancel the functional determinant that appears due to 
the change of the integration variables iIj{t) — > ip{iuj). Introducing the shorthand no- 
tation {'tlj,Xip) := Ylik k' '^kXk,k''4^k' where k = (u;„,Z) allows for rewriting (12. 5p in the 
simple form 



Z 

Zq Zi 



2, 



Vk[,k'2MM^K'^k'^^k2i^k^ > 

M ) (2.6) 

\- j PV'V'exp - Sintj. 



^We have chosen a basis where is diagonal. This is, however, not a vital condition, since all 
properties of the noninteracting problem (like Wick's theorem) follow from the Gaussian nature of 
the functional integral l|2.5p with u = 0, diagonal or not. The shorthand notation introduced in l|2.6p 
therefore explicitly allows for non-diagonal terms in Q'^ . 
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A factor of /? ^ as well as the frequency-conserving (^-function have been absorbed into 
the anti-symmetrised two-particle interaction v. 

2.1.2 Generating Functionals of Green Functions 
Connected Green Functions 

We define the generating functional of the m-particle Green functions as 

mM,W) = ^ y Wexp{5o-5int-(V5,??)-(r?,V)}. (2.7) 

By construction, taking the functional derivative with respect to the external source 
fields rj and f/ yields the m-particle Green function in frequency space^ 



GjniK, ■ ■ ■ ,k'^;ki,. . . ,km) ■■ = (-1)™ (i^k[ ■ ■ ■ i^k'^i^km • • • V'fci) 



-w(M,M) 



^r]k[ ■ ■ ■ ^Vk'^ ^Vk^ ■ ■ ■ Sr]ki 



(2. 



ri=ri=0 



In the noninteracting case (|2.5p is a Gaussian integral that can be solved easily leading 
to the well-known result for the one-particle Green function, 

Gi{l,iuJn) = GfiioJn) = -. , 

tUJn - (-1 

and to Wick's theorem to calculate all higher-order functions diagrammatically using 
Feynman diagrams. 
The functional VV^ that generates the connected Green functions 

G'rn{k'^,...,k'^;ki,...,km) ■■ = (-1)"" (^ipk[---A'^i^k^ ■ --i^ki^^ 



Sim 

-W^(M,{r?}) 



(2.9) 

r?=?7=0 



^'nk{ ■ ■ ■ ^m'^ ■ ■ ■ ^Vki 

is given by (for a proof see [Negele&Orland 1988]) 

W'=({r/},M) = ln[>V(M,{77})]. (2.10) 

Vertex Functions 

Next, we consider so-called m-particle vertex functions which are (in a diagramatical 
fashion) defined to consist of all connected one-particle irreducible diagrams with the 
m external legs amputated. In particular, this implies that they cannot be split into 
two pieces by cutting one single-particle line. 
The generating functional F for these vertex function is given by a Legendre transform 

of 

r({<^} , W) = -miv} , {v}) - {4>, v) - iv, <P) + {4>, </>), (2.11) 



^From j:;^^ = / ^2il^_iw_rf^' ^ p-^/^ j e''^"^' ^^^dr' it follows that by differentiating 
the generating functional W with respect to the Fourier-transformed fields one obtains the Fourier 
transform of the imaginary-time Green function. 
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where the new fields 4) and (j) are defined in the standard way, 



(2.12) 



and the last term in (|2.1ip was added for convenience. Of course it is far from obvious 
that the vertex functions can really be obtained from this functional by taking the 
derivative with respect to the external fields. Anyway, we define 



7m (^1) • • • J ^m> ■ 



(2.13) 



2.1.3 Relations between Vertex and Connected Green Functions 

We will now explicitly show that the one- and two-particle vertex functions are in fact 
obtained by (|2.13|) (or, to say it the other way round, that the so-defined functions 71 
and 72 consists only of connected one-particle irreducible diagrams with their two (four) 
external legs amputated). Therefore we first calculate 



5W^ 5r]g 6W^ 6f}g 



6rjq 64>k drjq Scpk 



+ 



6r]q 6rig 



01-1 

q,k 



+ Vk 



f2l2t 



6W^ Srjq _ 6W 5% 
6T]g 64)h 



6rjq 6r], 



(2.14) 

+ Vk 



JTT2t 



(2.15) 

The additional minus signs appear because we are taking the derivative with respect to 
Grassmann variables. One should be reminded that as we have turned to the Legendre 
transform F, (f) and (j) are independent fields, implying tt = etc. Next, we consider 



arm 



6 5W^ 
Hk' 5f]k 



5^T 



[G 



01-1 



E 



6r]q 



+ 



6f]q 



Hk' 5riq5'qk S(pk' ^Vq^Vk 



1''''/ 6r]g5rik S(j)k'6<j)q 6r]q6r]k 



(2.16) 



The last line follows from differentiating (|2.14l) and (12.151) by cp and (p, respectively. 
Performing similar calculations that start out with = Si. u and 4$^ = tt^ = 0, 
and denoting the resulting equations in a compact form (by interpreting the derivatives 



as matrix indices (fci, k2)), we obtain 



&<t>&<t> 



5^T 



01-1 



SfiSr) 
SrjSr] 



SfjSfj 
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As we will need it later on, we furthermore define 



SfjS'r] 



52 r 



5r)&rj SrjSri 



52 r 



52 r 



T 



(2.17) 



It is now possible to establish a relation between the connected Green functions and 
the vertex functions. Therefore we evaluate p.l7p at vanishing external fields and look 
at the (1,1) - element of the matrix equation which is now diagonal,^ 

-1 



Gi{k'-k) = G\{k'-k) :-- 







\ 5^T 






r]=r)=0 




<j>=(f>=0 



01-1 



(2.18) 



k',k 



The identity (*) is the so-called linked cluster theorem. Since Gi is just the single- 
particle propagator Q (calculated in presence of the interaction term 5int), (|2.18|) is 
nothing else but the Dyson equation. Consequently, the one-particle vertex function is 
related to the self-energy S, which is well-known to be one-particle irreducible,^ by 

^i = £i - =-^- (2-19) 

5(1)54) 4>=<t>=o 

Ignoring the difference in sign, we will call 71 self-energy from now on. 

To prove that the same (one-particle irreducibility) holds for the two-particle vertex, 
we differentiate (|2.16p with respect to cpi, which leads to 



+ 



5^r 



5^T 5^W^ 



')i5(j)k'S(t)q67]q5r]k yScpk' 



[^°]q,k' 



E 



6t]s 5fis 5^W 



+ 



/ 6r]s5r]q6r]k S(t)i 5r]s5r]q6r]k 



Differentiating once again with respect to (pii and setting all external fields to zero 
simplifies this equation enormously, since all terms that do not contain an equal number 
of field derivatives 5^ and 5i vanish. We obtain 



5^W 



6(j)i'54>i54>k'5(t)q 5r]q6r]k 
5^r I 



+ 



5^r 



5r]s 6r]s 



5^W 



Sk,q - 



i^^o.k'i^h.li^^l^s' 



6(1)1 5(j)v 6r]s'57]s6r]q5r]k 
5^W 



4,=(l}=o 



Y,^2{l',q;k',l)g,,q- ^ [[Q]q},,[Q];}[Q]j,l,G'2(.s',k,q,s] 



(/)=0=O 



where we have used 



q,s,s' 



SVs 1231 r^^Ol-l 
A/A. 1-'=' \s,k 

01-1 
\k,s 



5^k ^ 



s,fe 



+ 



(2.20) 



(2.21) 



^The off-diagonal elements are zero unless we are in a phase of broken symmetry. 
■*To be more precise, for a quantity X, being the sum of all connected one-particle irreducible 
diagrams with the two external legs amputated is equivalent to fuUfilling the Dyson equation. 
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Solving (|2.20p for 72 yields 

i2{k',A-MM) = - E [^].t'?it^]S5^[^]^"^!'^^[^]^"ii^2(g'i,9^;9i,92), (2.22) 

which is the prove of the one-particle irreducibility of the two-particle vertex functions. 
Interpreted diagrammatically, (I2.22p just states that the two-particle vertex is obtained 
by cutting a full single-particle line at each external leg of every two-particle connected 
diagram. What is then left are those diagrams that cannot be further split by just 
cutting a single-particle line (since the split-off part would have been amputated before), 
which is precisely the definition of one-particle irreducibility. 

2.2 The fRG Flow Equations 

Having learnt how functional integrals can be used to compute the partition function 
as well as Green functions of a fermionic many-particle system, we have established the 
formalism to set up the fRG flow equations. In order to do so, we introduce a new 
independent variable A in the noninteracting propagator , 

a° ^ a°(A). 

Later on, we will think of A as an infrared cutoff in frequency space, 

g\iuJn,^) = Q{Wn\-^)Q\i^n), 

but for the derivation of the flow equations the actual form of the A-dependence is 
completely irrelevant, so that we will not specify it right now. 

In a very pragmatic sense, the fRG can be viewed as follows. The flow equations 
for any physical quantity X[^°(A)] (such as the generating functional and T) are 
obtained by taking the derivative with respect to A, 

A:(A) = y(A), (2.23) 

where Y{K) remains to be computed. If we now choose two values A'"'*'*^' and A*^°^' such 
that X[^''(A'°'*''^')] is easy-to-calculate, while the propagator evaluated at A*'"^' is just 
the ordinary free propagator Q^{K^^'^^) = ^ we can try to solve (I2.23P with the initial 
condition X(A'°'*'^i) = X[g°(A™"^i)] to compute our desired quantity X[gO(Afi°^i)]. A 
'simple' X[^''(A™'*'^')] is in particular obtained by cutting out all degrees of freedom, 
that is choosing ^0(A'°'*'^') = 0, which will be our future choice. 

If X depends on other independent variables ipi (such as r({(^}, {i;^})), we can expand 
it into a Taylor series with coefficients /ij describing our physical system (such as the 
self-energy), assuming that we have chosen an appropriate expansion point (pi. (|2.23p 
will then provide an infinite hierarchy of flow equations for the expansion coefficients. 
Solving this hierarchy is equivalent to solving the original flow equation for arbitrary 
values of the variables ipi . 

One could immediately argue that in general we are already unable to calculate the 
partition function of an interacting problem, and hence we will totally fail in solving 
infinitely many coupled differential equations whose solution would provide us with 
complete knowledge of our system of correlated fermions. This is certainly true. The 
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fundamental point is that for a suitable choice of X there is physical reasoning to trun- 
cate the infinite hierarchy of flow equations for the expansion coefficients, rendering their 
solution (at least numerically) easy. The generating functional of the vertex functions 
will turn out to be that suitable choice if we want to set up an approximation which 
treats the two-particle interaction perturbatively, justifying our extensive elaboration 
on it. 



2.2.1 Flow Equations of Connected Green Functions 



We will for the moment refrain from setting up a flow equation for the generating 
functional of the vertex functions and turn to the connected Green functions instead. 
The flow equations of the latter will be of no further interest, but they will facilitate 
the computation of the former. 

For convenience, we replace the full partition function in the denominator of (12.71) by 
the noninteracting one, so that we start out with the following functional 



W(M,{7?},A) 



1 



^o(A) 



j Wexp{5o(A) - 5i„t - (^,??) - (2.24) 



This replacement changes and F only by the trivial constant In(Zo), while all higher- 
order functions G^^q and 7m>o remain unaffected. In order to derive the flow equation 
for the generating functional of the connected Green functions, we have to compute the 
derivative 



>V^(A) 



1 



W(A) 



9a 



1 



2o(A) 



T)ip^ exp 



{5o(A)-5i„t-(Vi,ry)-(r?»} 



(2.25) 



Therefore we first calculate 



9a 



2o(A) 



^^o(A)e-(^) 
^^o(A)e-(^) 



^o(A) 



^o(A) 



.50 (A) f 5o(A)e^«(^) 



^o(A) 
e5o{A) 
2o(A) 



(^,aA[g°(A)]-V) -E| / -^Al^^A)]-], 



[(V;,9A[g°(A)]"V)-TV(gO(A)aA[g°(A)]-i)] 



which allows for rewriting (|2.25|) in the simple form 



VV=(A) = -IV (5°(A)8A[ff"(A)]-l~ 
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If we eliminate W in the second term, 



Sr] 
6r] 

6r] 



5ri 

^5W%A) 
6f] 

.i5W^{A) 
5f] 



W=(A) 



+ [±,dA[g'{A)]-'^_]W%A) 
or] or]/ 



k,k' 



6'W{A) 



^6^W%A) 



Tt{dA[g'>{A)] 



6r]Sr] 



we are finally able to write down an analytic expression for the flow equation of the 
generating functional of the connected Green functions, 

_^6^W%Ay 



W^(A) = - IV (gO(A)aA[e?°(A)]"i) + Tr 5A[a°(A)] 



6r]6r] 



(2.27) 



6r] 6r] 
Having solved this differential equation, we can calculate the grand canonical potential,^ 



-rinZ 



-TGl{A = Ai)-T\nZQ. 



Since the flow equation for the generating functional of the Green functions follows 
directly from (|2.26|) . we write it down for reasons of completeness 



W(A) 



T,{g\A)dA[g\A)] 



+ 



f,aA[g°(A)]-^A 

or] or] 



W(A). 



2.2.2 Flow Equations of the Vertex Functions 
Flow Equation of the Generating Functional 

We will now compute the flow equations for our primary quantity of interest, the vertex 
functions. We start out by simply differentiating (12. lip with respect to A, bearing in 
mind that the fields r] and f] have to be expressed in terms of 4> and via (|2.12|) and 
therefore acquire a A-dependence. We obtain 



f({<A},W,A) 



d 

"dA 



W'{{f]{A)],{r]{A)],A) - ilvW) - (^(A), 



+ 



,dA[g 



Oi-l, 



E 



6W , 5W 

m^— + Vk- 



Sr]k 



Svk 



!.,7))-(^,<A) + U.,5A[g^] 



Oi-l, 



arm 
arm 



-W^({^(A)} , {r?(A)} , A) + (<^, 5A[a°]^ V) 



TV (gO(A)aA[a°(A)]-i) - TV ( 5A[g°(A)] 



TV(gO(A)5A[a°(A)] 



-1^ 



6f]6r] 

Tv(9A[a°(A)]-^V(i'^)(<^,<A,A) 



(2.28) 



'Remember our replacement Z — ^ Zo in the denominator of l|2.24p . 
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where the dot on top of VV"^ denotes the outer derivative, and we have defined V^^'^^ as 
the (1,1) - element of the matrix equation (|2.17p that relates the derivatives of to 
those of r (at arbitrary values of the external fields). 



Flow Equations of Vertex Functions 

In order to expand it into a Taylor series, we rewrite (|2.17p as 



V(A) 



01-1 



' ~ [[a(A)]-il^l + l 



1 



^ -U{A) 



-g{A) 

[^(A)]^ 



^ -U{A)) 



-1 



-a(A) 

[g{A)f ' ' 



= :V 



with U{(t),(j), A) being the difference of the one-particle vertex function and the second 
derivative of F evaluated at arbitrary external fields, 

6'^T{A) 6'^T{A) 



U{A): 



(j}=(l>=0 



6^T{A) 



71 (A) 



4£i^-[g(A)]-i + [gO(A)]-i. 



Defining V := iX^'^), (12:28]) can be cast in the following form, 

r(A) = IV (gO(A)9A[a°(A)]"i) - IV (g(A)5A[a°(A)]"iV(A)) . (2.29) 

Of course there is no new physical insight in (|2.29p as it is just a different formulation 
of (|2.28|) with renamed variables. But these new variables allow for a Taylor expan- 
sion, which will be useful to derive ordinary differential equations from the functional 
differential equation (|2.29p . Namely, 



v(A) = i-gu + gugu-g~g^^-^^ + ... 



(2.30) 



where we have written down all terms up to second order of the geometric series. 
Next, we expand r(i^, cp, A) around cp = cp = 0, 



i-iy 



X] J2 7m(fc'i,...,A;^;A:i,...,fcm;A)(/)fc/ 



m=0 ^ ' k[...k'^ki...k„, 

which by substitution into (|2.29p sets up flow equations for the physical meaningful 
expansion coefficients 7m • For 70 we obtain 

70(A) = IV (gO(A)9A[a°(A)]-i) - IV (g(A)5A[gO(A)]-i) , 
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which is the first term in (|2.29p together with the zeroth order contribution of V to the 
second term, since U as well as tttt and tttt are at least of second order in the fields. 



Flow of the Self-Energy 

To set up the fiow equation for 71, we have to compute the part of V that is linear in 
(^(/>. It is given by 

-•15^ 

= 7^ X J Aw^ X] I2{k'i,k'2;ki,k2]k)(t)k'^(t)k'^(t)k2(t>ki 

= - ^'y2{k',q';k,q;A)(j)k'(l)k, 

k',k 

where we have used that by construction it follows that 

72(/c'i,A;2; /ci,fc2; A) = -72(4, /ci, ^2; A), 

and likewise for the interchange of ki and k2. Comparison of the coefficients of the 
term linear in i^t/) on the left- and right hand side of (|2.29p now yields the desired flow 
equation, 

7i(A:';A:;A) = TV(g(A)5A[g°(Arig(A)72(fc', ;k, ;A)), (2.31) 

with 72 (A:', ; /e, ) denoting a matrix with indices [72 (/c', ; /c, )\q>,q'-= ^2{k' ■,q';k^q). The 
self-energy couples directly into this equation via the full propagator, ^(A). 

The Flow of 72 and Higher- Order Functions 

In order to derive a flow equation for the two-particle vertex, we have to find all terms 
on the right hand side of (|2.29p (or rather of (|2.3Qp ) containing four external fields. 
Again, by comparison of coefficients we obtain 

72(A:'i,A;^;fci,fc2;A)=Tr(g(A)5A[a°(A)]"ig(A)73(A;'i,A;^, -MM, ;A)) 

-Tr(g(A)5A[g°(A)]-ig(A)72( , ; A^i, ^2; A)[g(A)f 72( , ;A:i,fc2;A)) 



-Tr 



g(A)5A[g°(A)]-ig(A)72(fc;, -ki, ;A)g(A)72(A;^, ;A:2, ;A) 
— [k'l /C2] ~ [^1 ^2] + [^'1 ^ k2,ki /c2] 



(2.32) 

The coupling to the three-particle vertex 73 emerges from the part that is proportional 
to 4>4>4'4' in the second term of (|2.30p . while the second term results from the quadratic 
part of Oj^G^^-^^- (the fourth term of (I2.30p ). The last four terms arise from 
QUQU (third term of (I2.30p ) after anti-symmetrisation of the coefficient of the quadratic 



®For and ^^^^ this follows because all terms in the Taylor expansion that not containing an 

equal number of fields 4> and 4> vanish (unless we are in a phase of broken symmetry), while for lA all 
zeroth order terms are cancelled by definition. 
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FIG. 2.1: Diagrammatic representation of the flow equation for the self-energy (|2.3ip . 

part (this is necessary since it is not anti-symmetric by itself; or put differently: anti- 
symmetrisation will naturally occur when writing down all possible ways of matching 
the coefficients on both sides of ()2.29|) ). 

We could pursue this game forever and forever, writing down flow equations for the 
vertex functions of arbitrary order that would always turn out to couple to even higher- 
order functions. After having derived all these expressions (which is impossible) we 
would have to solve infinitely many coupled differential equations (which is impossible 
as well). Fortunately, having considered vertex functions will turn out to have been a 
very clever choice, because their fiow becomes neglectible with increasing order if one 
assumes the two-particle interaction to be small. In fact, it will be sufficient to treat 
the infinite hierarchy only up to the second order (setting all other vertex functions to 
zero), and the corresponding fiow equations were derived in this section. 

The general structure of the higher-order equations, however, should have become 
clear by now. The fiow equation for jm contains one term involving "Jm+i as well as 
contributions from all lower orders (except from 70). For example, the fiow of 71 (|2.3ip 
is determined by 72 and by 71 itself (via G). 

Feynman Diagrams 

It is often convenient to denote the matrix summations appearing on the right hand side 
of an arbitrary-order fiow equation in analogy with Feynman diagrams. A m-particle 
vertex is symbolised by a dot with 2m external lines, and there are two (and only two, 
as should be clear from (|2.30|) ) combinations that involve 71 on the right hand side of 
(|2.29l) . namely ^(A) and Q{A)d\[Q{A)]~^Q{A) , and they are represented by a single line 
or a crossed-out line, respectively. The quantity Q{A)dA[G{A)]~^Q{A) =: 5(A) is also 
called single-scale propagator. As for ordinary Feynman diagrams, summation over all 
internal lines is implicitly assumed. An example for the fiow of the self-energy is shown 
in Fig. O 

Initial Conditions 

Up to now we have not thought of the initial conditions that are necessary for the 
solution of all the aforementioned differential equations to be well-defined. It is of 
course possible to derive these conditions analytically [Meden 1998], but here we only 
give a very simple diagrammatic argument. Since at the beginning of the fiow we want 
to cut out all degrees of freedom, that is ^0(A'™*'^^) = 0, all diagrams containing a 
finite number of free propagators are zero. The only nonvanishing diagram is the pure 
interaction v (its external legs have been amputated by definition). Therefore the 
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initial generating functional reads 

i,j,k,l 

and the initial conditions for flow equations of the vertex functions are given by 

72(A;'i,fc^;A;i,fc2;A™*'^') =%i,,^,fc,,fc, (2.33) 

for m = 2, and 

7™(A™"'^i) = (2.34) 

otherwise.^ 

2.2.3 Truncation Schemes 
General Arguments 

As mentioned above, we need to truncate the infinite hierarchy of flow equations for the 
vertex functions in order to render it solvable. In this section, we will set up a truncation 
scheme that is valid in the limit of a small two-particle interaction v. It will turn out, 
however, that despite its perturbative nature the resulting finite set of equations will 
be sufficient in describing the effects of 'fairly large' interactions even on a quantitative 
level in all situations considered in this thesis (we will discuss the limitations of this 
approximation in the last chapter). 

The right hand side of the flow equation for each vertex function contains contri- 
butions from 7m+i as well as from all lower-order vertices. At the beginning of the fRG 
flow, all functions 7m except for the two-particle interaction 72 (A = 00) = v vanish, 
so that all vertices 7^ with m ^ 2 are generated only by v. Now, the fundamental 
point in setting up a truncation scheme perturbative in the interaction is the following. 
Since we are considering vertex functions, an m-particle vertex 7m (A) has to be irre- 
ducible for arbitrary choice of A. In particular, it can only be generated by irreducible 
diagrams on the right hand side of (|2.3Q|) . Since there are no such diagrams with 2m 
external (amputated) legs that contain less than m terms v,^ all vertex functions 7m^2 
are generated by terms that are at least of order m in the interaction. If the latter is 
initially small and stays small for all A it is justified to cut the infinite hierarchy of flow 
equations at a certain order, neglecting the flow of all higher-order vertex functions. 



Approximation Schemes 

The very simplest approximation that emerges from these considerations is to neglect 
the flow of all vertices except the one for the self-energy. This leaves us with one single 
equation (remember that G = +7i]~^)) 

7i(fc';A;;A) =Tr(g(A)9A[a°(A)]-ig(A)72(A:', ;A = 0)) 

= ^ [g(A)5A[g°(A)]-ig(A)]^^^,%,,,fc,„ (2-35) 

^Sometimes one absorbs some additional one-particle terms Vk\k into the interacting part of the 
action (iSint). The initial condition for 71 in then reads 71 (fc'; fc; A'"'*"'') — Vk'^k- 

* Consider for example the three-particle vertex. There is one possibility to write down a diagram 
with two interaction vertices and six external (amputated) legs, but this diagrams contains one single- 
particle line connecting the interaction vertices, rendering it reducible. 
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which is easy to tackle numerically (and sometimes can even be solved analytically). 
Since the bare interaction v is frequency-independent, the self-energy will not acquire 
a frequency dependence during the flow, so that the effect of the interaction is just 
to renormalize the single-particle energies of our system by 7i(A*^°^'), allowing for an 
(effective) one-particle interpretation. 

The next logical step would be to consider the coupled flow equations of the self- 
energy and the two-particle vertex. Unfortunately the resulting problem would not be 
frequency-independent any more, so that we would again have to solve infinitely many 
coupled equations (because we have infinitely many Matsubara frequencies). Therefore 
we implement another approximation: we neglect all frequency-dependencies of the two- 
particle vertex except for frequency conservation, which we have explicitly included in 
V and which is conserved by the flow equation (since is diagonal in frequency space) . 
As the bare interaction is frequency-independent, this will lead to errors only of second 
order for the self-energy, and of third order for the two-particle vertex. Alltogether, 
(|2.31l) and (I2.32p can then be recast as 

7i(fc^ k; A) = ^^[ pjiLOn, A)dA[g\iiOn, A)]"^g(ic^„, A) ],,,,72(fc^ q'; k, q; A) (2.36) 

iuJn q,q' 



::5(ia.)„,A) 



for the flow of the self-energy, and 
i2{k'i,k'2]ki,k2]K) = 



E E 



Sq,q'{iuJnA)l2{(i , s' ] ki,k2] A)^^y (-iw^, A)72 (/c'^ , ; s, g; A) 

^Sq^qi{iuJn,^)l2{k'i,q] ki, S] k)Q s^s' {^^n, ^)l2{k'2, s' ] /C2, q] A) 



[k'l <-> — [ki ^ ^2] + [k'l ^ k2, ki ^ A;2]| 



(2.37) 



for the flow of the effective interaction (as we will call the frequency-independent two- 
particle vertex from now on). Here all the labels ki only denote the single-particle 
quantum numbers, the frequency-dependence has been written out explicitly. It will 
turn out that the remaining summations over the Matsubara frequency iuJn can be 
easily performed if we introduce a sharp cutoff in frequency space. 

Also this more elaborate approximation scheme yields a self-energy 71 
is frequency-independent, which will again facilitate gaining further insights into our 
upcoming results as it allows for an interpretation in a simple effective single-particle 
picture.^ 



Symmetries 

By applying the aforementioned approximations, we have boiled down the infinite hi- 
erarchy of flow equations to N"^ + coupled ordinary differential equations, where 

word of warning is in order. If the whole effect of the interaction is to renormalize the single- 
particle energy of our system, the properties of this effective system can of course be interpreted very 
easily due to its noninteracting nature. This does, however, not answer the question why the interaction 
particularly renormalizes the levels the way it does. 
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N denotes the number of single-particle quantum numbers that define our interacting 
system. For all considered in this thesis {N < 8), this set can be solved very easily by 
direct implementation on a state-of-the-art personal computer. However, if one wants 
to even further speed up the numerics, one should realise that many flow equations are 
still redundant since up to now we have not exploited any symmetries. 

The most obvious symmetry to reduce the number of independent equations is the 
(anti) symmetry of the two-particle vertex. By definition, the bare interaction v is sym- 
metric under the exchange of the first with the last two indices, and it is antisymmetric 
under the exchange of the first or the last indices, respectively. These symmetries 
are conserved by the flow equation (|2.37p . For the first term on the right hand side 
conservation of the symmetry is seen if one considers 

Sg^g'{iujn, A)72(g , s'; k[, k'2; A)^s,s/(-iu;„, A)72(fci, ^2! s, A) 
= Sg^gi{iujn, A)72(s, q; ki,k2] h)Qs,s'{-i^n: A)72(A;i, fcg; q' , s'; A) 
= A)72(g, s; ^1,^2; A)gs' A'^^n, ■^h2ik[, k2; s' , q'; A), 

where we have used Qs,s'{'i'^n) = Ss^si—i^n) as well as (|2.39p . Renaming the sum- 
mation indices {q <-> <-> s',iu!n — > —ioJn) then yields the desired symmetry under 
exchange of the first with the last indices. For the second term this follows directly 
from the definition, and it is obvious that ()2.37|) preserves the anti-symmetry of the 
two-particle vertex as well.^'' 

The second important symmetry of the equations (|2.36p and (|2.37|) is spin conserva- 
tion. If the bare interaction and the free propagator are spin-conserving, no terms that 
are non-diagonal in spin space can be generated by the flow. 

Comparison to Ordinary Perturbation Theory 

In order to render the infinite hierarchy of flow equations solvable, we have applied a 
truncation scheme valid in the limit of vanishing strength of the two-particle interaction. 
One could now immediately ask why such a perturbative approach should be superior 
to ordinary perturbation theory (which does not require 20 pages of elaboration on 
functional integrals to be set up). This, however, can already be seen by considering 
only the very simplest approximation, the flow equation for the self-energy with the 
two-particle vertex set to its initial value v (I2.35|) . 

7i(A;'; A;;iu;„; A) = A)5A[^°(iu;„, A)]~^^(iu;„, A)] vy^q'^k^q 

iuJn q,q' 

[-^A[^°(^'^"'^)]~^]g,g' ^k',q',k,q- 

iun q,q' 

In the second line we have furthermore replaced the full propagator Q by the noninter- 
acting one, neglecting the self-energy. The integration can then be performed trivially, 
leading to 

iuin q,q' 

^°The two-particle vertex is of course antisymmetric for arbitrary choice of A. But since we are only 
computing an approximation for the exact vertex it is, however, not from the beginning obvious that 
this approximate vertex obeys the same antisymmetry relations. 
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which is precisely the expression that would have followed from first order perturbation 
theory, but it is only recovered because we have made a further approximation within 
our original flow equation. To say it the other way round, (|2.35|1 has to correspond 
to a summation of more than first order diagrams. In fact, it will turn out that in all 
situations considered in this thesis, the functional renormalization group will produce 
results far superior to those of perturbation theory. 



2.2.4 Specification of a Cutoff 
Introduction of a Sharp Cutoff 

The only remaining step now is to specify a certain form of the A-dependent propagator 
^''(A). As mentioned above, we choose A to be an infrared cutoff in frequency space, 
namely 

A) = Q{\ujn\ - K)g\iuon). (2.38) 



We have employed a 'sharp' 6 - function cutoff to simplify our calculations and to speed 
up the numerics^^ (in particular to make it possible to carry out the summation over 
the Matsubara frequencies in (|2.36p and (I2.37P analytically) . Since we want to cut out 
all degrees of freedom at the beginning of the flow, we choose A™'*'^' = oo, and we set 
^finai — Q recover the original free propagator. 



Carrying Out the Freqency Integrals: Morris' Lemma 

In order to perform the remaining frequency summation in the flow equations for the 
self-energy and the effective interaction, we consider 



S{iujn,^) = ^(icu„, A)9A[^°(iw„, A)] ^^(zu;„,A) 



1 + a°e7i 
(i + goe7i)2 

6de 



1 + a°e7i 



1 + g^e-fi 



A)dQg{iujn,A). 



(2.39) 



In the zero-temperature limit, on which we will mainly focus within this thesis, we can 
write the summation over Matsubara frequencies as an integral, J2iuj„ ~^ {2-kT)"^ f duj, 



^^Of course the special choice of the cutoff should not influence the underlying physics. It can 
nevertheless effect the quality of the results obtained. We will, however, not elaborate on this issue 
here. 
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and with the help of ()2.39|) . (|2.36l] becomes 

7i(A;';fc; A) = ^ j duj'^Sg^g>{iuj,A)-f2ik',q';k,q;A) 

<i,<i' 

duj ^6{\uj\ - A)deGq,q'{iuJ, A)72(/c', q'; k, q; A) 
q,q' 

(icu^(5(|cj| - A)^ dt {dtgg,q'{iuJ,A;e^t)}-f2{k',q;k,q;A) 



1 

2^ 



2tt 



1 

2^ 



rf^^ ^ (^(1^1 - A)^g^5/(iu;, A)-/2{k',q';k,q;A) 
<i,g' 



2tt 

ll1=±A q,q' 



(2.40) 



The factor of T was cancelled by one implicitly contained in the two-particle vertex 72 
via the initial condition (|2.34p (because we have initially absorbed a factor of (3~^ into 
the bare interaction which we will from now on consequently think of not containing 
the temperature any more), and we have introduced 

^(A) : ^ 



[gO]-i+7i(A)- 

The identity (*) is the so-called Morris' lemma that allows for evaluating the at first 
sight ambiguous product of a S- with a 0-function, assuming that the sharp cutoff 6 is 
implemented as a limit of increasingly sharp broadened cutoff functions Oe^o [Morris 
1994], 

6,ix - A)f[Q,{x - A)] ^ 6ix -A) [ f{t)dt, 



with / being a continuous but otherwise arbitrary function. The functional form of 
for finite e does not affect the result in the limit e — > 0. 

To show how the frequency summation in the flow equation for the effective interaction 
is performed, we examine 

<Sq,q'{iuJn, A)^^,^' (^F^l^n, A) 

iu!„ 

=^ / duj 5{\uo\- A) dt {dt[gq,q,{iu;,A,e ^t)]g,^s'{TiuJ,A,e ^t)} 

=^ J duj 5{\lo\-a) dt |^at[g,,q.(^'^,A,e^t)g,y(^icj,A,e^t)]| 

=— / (itj (5(|u;| - A)^g^g/(iu;, A)^s^^/(=Fiw, A) 
=^ ^ gq,q'iiuj,A)gs^s'{Ti^,A). 

LO=±A 

The third line follows because (I2.37P is symmetric under the exchange of the summation 
indices (s <^ q,s' q' , iujn T^^^n) (by precisely the same argument used to show the 
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symmetry under the exchange of the first and the last two external indices). Altogether, 
the flow equation for the effective interaction can then be cast in the form 



72(A:'i,fc2; fei, /c2; A) 
9^ ^ ^ 



LJ=itA s' ,q' ,s,q 



-^Gq,q'{iuJ, A)72(g', s';ki,k2; A)gs,s'{-ii^, A)l2{k'i^ k2;s,q;A) 



- Qq,q' {iio, A)72 (/c'l ,q';ki,s; A)Gs,s' {i^, A)72 (4 ,s';h,q; A) 
+ Qq,q' (iuj, A)72(/c2, 9 ; S] k)Qs,s'{i'^, k)-^2{k[, S ] k2,q] A) 



(2.41) 



The special form of the cutoff has completely disappeared from both (|2.4QI) and (|2.4ip . 
ruling numerical problems due to the 'discontinuous' (5-function out from the beginning. 



Finite Temperatures 

At nonzero temperatures one cannot replace the frequency summation in (|2.36p and 
(|2.37|) by an integral. However, using Yliu) — I ~ ^n) we can recast the flow 

equations for A = {71,72} as 

A = [ — 27rry 5(cj - w„)e(|(j| - A)s(e,w), 

J 2n ^ 

with an appropriate choice of B. If we implement the sharp cutoff O as a limit of 
increasingly sharp cutoff functions we cannot apply Morris' lemma because the right 
hand side of the flow equations is now discontinuous. Choosing a smooth cutoff instead 
would though not giving rise to conceptional problems considerably slow down the 
numerics. One way out of the misery that was originally proposed by [Enss 2005b] is 
to take the limits in a different order. If we first replace 5{u — ujn) by a continuous 
5-function 6x{co — cOn) sharply centred around a; = we can implement the sharp 
cutoff and apply Morris' lemma, 

A)B{u;) 



again assuming an appropriate B. If we choose 8x to be a box of height (27rT) ^ and 
width 27rT, we obtain 

^ = ^ E ^('^-)' (2.42) 

which is exactly the zero temperature flow equation, but with the continuous lo = ±A 
replaced by the discrete Matsubara frequency nearest to A, u;„ ±A. 



n 

/duj 
— 27rT^4(a; - Un)5{W 
n 

= rJ^5,(A- |a;„|)SK), 

n 
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Again: Initial Conditions 

We have dropped all convergence factors exp(iu;0~'^) when we initially replaced the 
imaginary-time functional integral for the partition function (|2.ip by its pendant in 
frequency space (|2.5I) . For finite A this exponential factor is irrelevant in the above flow 
equations, but we need it to define their initial conditions at A'°'*'^' = oo properly. This 
can be seen if one integrates (I2.40p from infinity down to some arbitrary large Aq, which 
yields a finite value, 

r-Ao 



1 /" " J- 

— lim / dA 5;e-o+g(iu;,A)72(A;',g';A;,g;A) 



=— lim / dA y Ye^'^o+^jML 

-I j^I^^T.^^'-^- [ A J--^ (2.43) 



q 




In the second line we have replaced the two-particle vertex by its value at A = oo, and 
we have used that Qq^qr{iu>) "^zlSP [Negele&Orland 1988]. With these considerations, 
the initial conditions (I2.33P and (|2.34l) read 

: -V,, , - 

2 ^ 



71 (A:'; k; A'^'^iai ^ oo) = -V^, ^ - ^ v,,^^^^,, 

g (2.44) 

72(/c'i, 4; /ci, fca; A'°"'^' ^ oo) = Vk'^^k'^^k^M^ 

with Vk'^k being an additional one-particle potential not included in the free propaga- 
tor.^^ We can now start the integration of the flow equations at some large arbitrary Aq 
discarding the exponential convergence factor, which is especially important because in 
general this integration has to be carried out numerically. 



Is Everything Weil-Defined? 

Finally, a few general words about potentially ill-defined expressions are in order. If we 
cut out all degrees of freedom at the beginning of the fRG flow, we have to take care 
of all those terms that contain a divsion by the free propag ator g°(A = oo) = 0. One 
such term is the first one in (I2.29jl , but it is cancelled by a contribution from the second 
term if we put in the expansion Q = — Q^^iQ^ + . . .. Otherwise, 9a [^'']^^ appears 
only in the well-defined combination 

^^As explained above, such a potential can be included in Sint on from the beginning and will afltect 
the whole fRG scheme only via the initial conditions l|2.33p and l|2.34p . 
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so that there are no ill-defined expressions in the flow equation (|2.29|) if we choose a 
sharp cutoff. 



^^At least if we implement it as the limit of increasingly sharp broadenend cutoff functions, as argued 
above. 



Chapter 3 

Quantum Dots 



Quantum dots, which is the usual name for mesoscopically confined electrons, are re- 
cently of great experimental interest. This is mainly due to their 'possible' application in 
nano-electron devices and for quantum information processing [Loss&DiVincenzo 1998]. 

The simplest realisation of a quantum dot system is the so-called single-electron tran- 
sistor (SET), which contains a 'droplet' of localised electrons coupled by tunnelling bar- 
riers to a sea of delocalised electrons (the 'leads'). An external voltage can be applied by 
source and drain contacts. If the dot is small enough one would expect its energy level 
spacing to be large, so that only a few levels must be considered, at least if the temper- 
ature is sufficiently low. The simplest model to describe such a situation theoretically 
is the famous single-impurity Anderson model (SIAM). If one considers the simplest 
case of only one level containing spin up and spin down electrons, in the zero tempera- 
ture limit one would intuitively expect a current to flow if the single-particle energy of 
this level Vg crosses the Fermi energy of the leads, the lineshape of the linear-response 
conductance GiVg) only broadened due to the coupling of the confined electrons to the 
latter. This follows indeed if we solve the SIAM assuming that the spin up and down 
electrons do not interact. The lineshape observed in the experiment [GGKSMM 1998] 
is, however, not Lorentzian but rather box-like consistent with the SIAM predictions 
if we take into account a sufficiently large local interaction. This enhancement of the 
conductance is in contrast to the simple picture that the usual effect of the interaction, 
the Coulomb repulsion, in combination with the spatial confinement of the electrons 
should lead to charge quantisation and Coulomb blockade transport properties. It is 
due to the Kondo effect which is generally active below the Kondo temperature if a dot 
with arbitrary many levels and the level spacing being much larger than the level broad- 
ening is occupied by an odd number of electrons. Thus Kondo correlation physics is the 
basis for the use of a single quantum dot as a nano-transistor, since the conductance 
properties can be manipulated by adding or removing single electrons. 

In this chapter we will introduce a Hamiltonian that describes more complex geome- 
tries containing several dots and/or several levels with arbitrary level spacing. We will 
assume two-particle interactions to be present between the localised electrons in the 
dots (the 'interaction region') while we model the source and drain as noninteracting 
semi-infinite tight-binding leads. For application of the fRG, the latter have to be pro- 
jected out in order to get a finite set of fiow equations. Finally we will show how the 
conductance through the dots can be computed in an approach beyond single-particle 
scattering theory which is used to describe transport if no correlations are present. 
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3.1 General Setup 

3.1.1 Experimental Realisation 

A frequently used method to fabricate quantum dot systems (such as the aforemen- 
tioned SET device) is to employ GaAs/AlGaAs heterostructures which contain a two- 
dimensional electron gas that results from the electronic properties of the different 
layers. Adding metallic gates and applying a negative voltage excludes the electrons 
from regions right below these electrodes, rendering it possible to separate a small re- 
gion (the quantum dot) from the rest of the heterostructure by tunnelling harries which 
strength can be tuned by changing the applied voltage. The energy of the dot relative 
to the electron gas can be controlled by an additional gate electrode with potential Vg. 
By adding drain and source contacts one can then measure the linear response conduc- 
tance G as a function of Vg by switching on a small bias voltage. A scanning electron 
microscope (SEM) picture that shows a SET used to measure the G{Vg) dependency 
at very low temperatures exhibiting the typical Kondo box-like lineshape is shown in 
Fig. 13.11 Technical details on the fabrication of the quantum dot system can be found 
in [GSMAMK 1998]. 

At T = 0, which is a meaningful limit since in experiments the temperature can be 





FIG. 3.1: SEM picture of a GaAs/AlGaAs 
heterostructure and metallic gate electrodes. 
The electrode on the right and the upper 
and lower ones on the left define the quan- 
tum dot and control the strength of the tun- 
nelling barriers to the rest of the electron 
gas (the 'leads') by varying the magnitude 
of the applied voltage. The middle electrode 
on the left can be used to change the energy 
of the dot. Source and drain contacts are 
not shown. (Reprinted by permission from 
Macmillan Publishers Ltd: Nature 391, 156, 
copyright 1998.) 



FIG. 3.2: SEM image of a two-path inter- 
ferometer formed on a GaAs/AlGaAs het- 
erostructure. Electrons can travel from 
the emitter (E) to the collector (C) either 
through the reference arm (right path) or 
through the quantum dot (QD, left path) 
which are separated by an island. The 
ground regions (B) are based and serve as 
draining reservoirs, ensuring that interfer- 
ence is only between the two path. The en- 
ergy of the QD is controlled via a plunger 
gate voltage. An additional quantum point 
contact (QPC) is capacitively coupled to the 
QD, allowing for measurement of its occu- 
pancy. All other electrodes are used for 
fine-tuning. (Reprinted by permission from 
Macmillan Publishers Ltd: Nature 436, 529, 
copyright 2005.) 
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tuned to be the smallest energy scale in the system (for example about one percent of 
the level-lead hybridisation strength in [AHZMU 2005]), the conductance G is connected 
to the transmission probability T by a constant of proportionality. The transmission 
probability itself is the absolute square of the transmission amplitude T = |ip that is 
related to matrix elements of the propagator by ordinary scattering theory. We will for 
the moment stick to this single-particle picture. Further comments on how to calculate 
the conductance of an interacting system and the connection to the ordinary scattering 
theory approach will be given in Sec. 13. 2.21 The transmission phase is now defined as 
the argument of the complex number t, which is nothing else but the phase change of 
the wavefunction of the electron passing through the system. Thus to fully describe 
transport both the transmission probability and phase have to be determined. Theoret- 
ically the computation of the latter poses no problem as soon as the exact propagator 
of the system is known. Experimentally it can be extracted by placing the quantum 
dot in one arm of a double-slit interferometer and piercing the whole system by a mag- 
netic flux The transmission probability Tq through the whole system is then given 
by Tq = |t -|- e^^'^tsp, where tg refers to the amplitude through the second path of 
the interferometer and A93 is an additional phase difference induced by the flux. The 
interferometer has to be sufficiently open to avoid multi-path interference. Assuming 
fully coherent transport, as should be expected at sufficiently low temperatures, the 
interference term in this expression reads \t\\ts\ cos(A(^ -|- arg(tcj) + arg(t)). The trans- 
mission probability is therefore expected to oscillate as a function of the magnetic field 
strength, as it is indeed observed in the experiment [AHZMU 2005]. If we assume the 
transmission through the second arm tg to be constant (in particular independent of the 
energy of the quantum dot), a change in the phase of t will lead to a similar change in 
phase of the oscillations, which therefore allows for a direct measurement of the former. 
An SEM picture of an experimental setup that establishes an interferometer with a 
quantum dot placed in one arm is depicted in Fig. 13.21 



3.1.2 Model Hamiltonian 

We will now specify a certain quantum-mechanical model to describe interacting quan- 
tum dots coupled to noninteracting semi-infinite leads. We assume that our general 
Hamiltonian consists of three parts, namely 

= -f^lead + -f^dot + -f^lead-dot • (3.1) 

Here -f/^iead describes the leads, i^dot the interaction region of the dots and i?iead-dot the 
coupling between the two. For simplicity, we assume the two leads to be equal and 
model them by a tight-binding approach, 

oo 

-f^lead = ~ 1^''" {^m,(T,s^rn+l,rj,s + Cm,cT,sC'rn+l,a,s^ + l^sC[n,(T,s'^rn,cT,s ; (3-2) 

s=L,R a m=0 

with Sm,u,s being the annihilation (creation) operator for an electron with spin direction 
o" = t, i localised on lattice site m of the left (s = L) or right (s = i?) lead, r denotes 
the hopping amplitude between two nearest-neighbour sites m and m -|- 1, and [is the 
chemical potential of the left or the right lead, respectively. For the time being, we will 
always assume our model to include spin degrees of freedom, since the Hamiltonian to 
describe spin-polarised situations just follows by dropping all spin indices in (|3.ip . 
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The Hamiltonian that describes the quantum dots is made up of three terms, 

-f^dot = ^^dot + ^dot + -^doV (3-3) 
The first term denotes the on-site energies of the different dot levels, 

^dot = EE^'.-<A-' (3.4) 

with dfl being the dot electron annihilation (creation) operators and e;^o- the one- 
particle dispersion^ of the states from the interaction region. In general, we will choose 
this dispersion to include a constant position (that can be different) for each level as 
well as a variable gate voltage Vg, 

ei,. = el^ + Vg. (3.5) 

Later on, we will mainly focus on studying the properties of the system described by 
(|3.1I) as a function of Vg. If we apply a magnetic field i3, the spin dependence of e^^ 
reads 

el = 6? - 4' (3-6) 

with the convention that a = ] corresponds to a = +1. The second term in (|3.3I1 
introduces a hopping between the different dots. 



<T = - E E [hi'dldi',<^ + H.c] , (3.7) 



o- l>l' 

while the last one accounts for the two-particle interactions, 

^^To^t = ^ E E Ki'' {<A,^ - I) {dlydi,^^ - . (3.8) 

We choose U to be symmetric, C/^'^^'f = U", f , and of course we set U^f = 0. The 
additional shift of the single-particle energies was chosen such that = corresponds 
to the particle-hole symmetric point. As suggested by our notation, we assume this 
energy shift not to be included in the free propagator. It will then affect the fRG 
scheme via the initial condition (I2.44p . which now reads 

■~fi{l' ,a';l,a;A = oo) = - E ^i,i'^cr,cr' - ^ E '^i''^' M,i"cr" 

' l",a" 
1",(t" l",a" 

We have used that in order to obtain (|3.8I1 . we have to define 



^Iii fact, this name is slightly misleading because H^°^ is of course a single-particle term as well 
which together with -ffdot would determine the one-particle dispersion ej- ^ . 
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FIG. 3.3: General setup described by the Hamiltonian l|3.ip . The two noninteracting leads 
(green) are modeled by a tight-binding approach with a hopping matrix element r, and they are 
connected to the interaction region (red) via tunnelling harries tj~^'^. The interaction region 
of the dots itself comprises levels with arbitrary onsite energies e/.cr, two-particle interactions 
U[l7 , and inter-dot hoppings ti^ii. The latter induce a spatial geometry into the system, and 
the whole Hamiltonian is qualified to model arbitrary many dots each with arbitrary many 
levels. 



Up to now we were quite sloppy in using the terms 'dot' and 'level' when referring 
to the region in between the noninteracting leads. A more precise notion would be 
the following. If we want to use (13.31) to model several spatially separated quantum 
dots each containing several energy levels, it is meaningful to introduce hopping matrix 
elements f/^// only between those /'s that belong to different dots. The choice of these 
hoppings therefore induces the spatial geometry of our system in the Hamiltonian. For 
each level there should be a local interaction Uj^f (at least if we want to include spin 
degrees of freedom), and we can introduce both inter-level and inter-dot interactions 
Upi^^ and U"l7 ■ Having all this in mind, we will continue the aforementioned sloppiness 
and speak of 'dots' when referring to the interacting region to keep the notation short. 

The last term in (13. 1|) that needs to be discussed is the coupling between our quantum 
dots and the noninteracting leads. It is modeled by 

^lead-dot = - E E [t!cl.,sdl,. + H.C.) , (3.9) 

s=L,R I 

where only the hoppings tf do not vanish if the corresponding dot I is connected to the 
lead s. The general setup is again depicted in Fig. 13.31 



3.2 Theoretical Approach 

In this section we will show how the fRG flow equations derived in the previous chapter 
are applied to our dot system. We will map the infinite system described by the general 
Hamiltonian (|3.ip to a finite one by projecting out the noninteracting leads. Finally 
we will show how quantities that might be measurable in the experiment are calculated 
in an interacting system. For our main quantity of interest, the conductance, we will 
derive a generalisation of the Landauer-Biittiker formula, the latter describing transport 
through a noninteracting system. 
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3.2.1 Application of the fRG 
General Considerations 

The starting point in order to treat ttie interacting many-particle Hamiltonian (13.11) 
are the coupled equations (I2.40p and (I2.4ip for the flow of the self-energy and the two- 
particle vertex evaluated at zero frequency. But since our system is infinitely large due 
to the semi-infinite leads, solving these equations would mean that one would again 
have to tackle infinitely many coupled differential equations. Fortunately this set can 
be boiled down to one of size A^^ -|- A^^ (not accounting for symmetries), where A^ is the 
number of degrees of freedom within the dot region, by the following considerations. 
Since the interaction region is finite at the beginning, it will remain finite during the 
fRG flow. In particular, no additional single-particle or interaction terms with indices 
outside the dot region can be generated by (I2.40p and (|2.41|) . since all external indices 
on the left hand side of these equations are connected to a two-particle vertex on the 
right hand side, which by deflnition vanishes for A = A'°'*'^' if one index is chosen 
from the leads. ^ Furthermore it follows from the usual expansion (with the implicit 
understanding of G referring to the interaction region only). 



(^0)^1 + 71 

that the free propagator entering the flow equations only needs to be evaluated at dot 
indices I, and hence we can replace by its projection on the dot region, PQ^P. The 
computation of the latter will turn out to be very easy by applying a standard projection 
method presented below. Finally, ^(A) can be calculated as the inverse of an A^ x A^ - 
matrix 

= (Pg°P)-.+7.(A) - 

It is important to note that the projection technique that leads to (I3.10|) is an exact 
procedure. 

Projecting Out the Leads 

The propagator of a noninteracting system'^ described by the Hamiltonian H reads 

1 



z-H 



In our case, H is obtained from the general Hamilton operator (13. 1|) by setting the 
interaction U to zero (and replacing the many-particle Hamiltonian by its pendant in 
the standard single-particle space). If we now define operators P and Q that split the 



^Of course this does not mean that the interaction in the dot region does not affect the leads. 
If, for example, we wanted to calculate the full propagator Qm.m' with lead indices m and m', there 
would always be a contribution from diagrams (if for the moment we think of the propagator exactly 
expanded in an infinite perturbation series) containing dot indices I since both are coupled via l|3.9p . 
It only means that the single-particle energies in the leads are not renormalized by the flow. 

^ Later on, we will also need this projection technique when the interaction in the dot region is 
present. Generalising the results derived in this section will, however, turn out to be quite simple, so 
that we will stick to the noninteracting case for the time being. 
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Hilbert space by projecting on the dots (P) or the leads (Q), we can rewrite Q as 

g^^^ (Spp{z) Qpq{z] 
\GQpiz) Gqq{z\ 

The projection on the dot's region Qpp := PQP can be computed by considering 

{z-H)Q{z) = 1 

{z- H){P + Q)g{z)P = P 

{zP - PHP)Pg{z)P - PHQQg{z)P = P 

{z-H)g{z) = 1 

^ {z - H){Q + P)g{z)P = P 

(zQ - QHQ)Qg{z)P - QHPPg{z)P = 0, 

where we have used that P and Q are projectors, that is P^ = P, = Q and P+Q = 1. 
By eliminating gqp = QgP from these two equations, we obtain 

gpp{z) = ^ 

zP - Hpp - Hpqjq:^^Hqp ^^^^^ 



with the obvious definition Hpp := PHP, . . ., and likewise for the other components: 

1 



Gpq{z) = - gpp{z)HpQ 



zQ-Hc 



Qqp{z) = - ^ ^„ HQpgpp{z) (3.12) 

zQ - Hqq 

GQQiz) =-7^ E7 — + E7 — HQpgpp{z)HpQ-^ 



zQ - Hqq zQ - Hqq zQ - Hqq 

Since only i?iead-dot will contribute when calculating HpQ and Hqp, we again write 
down its single-particle version i?iead-dot to avoid notational confusion: 

^lead-dot = "5^ 5Z [*?|0,cT,s)(^,a| + (tf)*|/,a)(0,cT,s|]. (3.13) 

o- s=L,R I 

The states \0,a,s) and |/,(t) denote the wavefunction with spin a localised at the last 
site (m = 0) of the left (s = L) or right (s = R) lead and of the dot level /, respectively. 
We can then compute the effective Hamiltonian Hqq that determines the propagator g 
projected on the dots, 

Hes = Hpp + HpQ— — — Hqp 

zQ - Hqq 

= Hpp + Y, E Y.[^'^tT4\l,a){<d,a,s\{zQ-HQQ)-'%a',s'){l\a'\ 

a,a' s,s'=L,R 1,1' 



HPP + Y. E E V'^tT^t'V^ ^) ^' ^K^Q - Hqq)-^% cj, s) {l',a\ 

a s=L,R 1,1' " ' 



(3.14) 
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The last line follows because the left and the right lead do not couple directly. Except 
for the calculation of the function g^iz) we have now succeeded in deriving a finite set 
of flow equations that describe the effects of the interaction in between the quantum 
dots. 

As we will need them later on, we also calculate the following matrix elements: 

gr.oSz): = {l,(y\GpQ{z)\0,a,s) 

= -E E Y.(^,<^\Gppm','^')it!'r{0,'^',s'\{zQ - HQQ)-^\0,a,s) 

a' s'=L,R I' 

I' 

gi.l{z): = {0,a,s\gQp{z)\l,a) 

= -E E Y.^O,a,s\{zQ - HQQr'\0,a',s')tf;{l',a'\gppiz)\l,a) 

a' s'=L,R I' 

= -E^?;«(^)' 
I' 

(3.15) 

for an arbitrary choice of s = L, R, and, 
gS^.o,{z): = {0,a,R\gQQ{z)\0,a,L) 

(0, a, R\izQ - Hqq)-'\0, a', s')4{l', a'\gpp{z)\l" , a") 
itt::r{0,a" ,s"\{zQ - HQQ)-'\0,a,L) 

l',i" 

(3.16) 

Here / is always assumed to be a single-particle index from the dot region. 
Calculation of g'^{z) 

Finally we have to calculate the propagator g"{z) of a semi-infinite noninteracting lead 
described by (|3.2p evaluated at the last lattice site. This is most easily achieved by a 
symmetry argument. Since both leads are assumed to be equally modeled by a tight- 
binding approach, their propagator g'^{z) will be of the same structure, so that without 
loss of generality we can consider the left semi-infinite lead with chemical potential 
ranging from — oo up to some arbitrary site m. The desired function gL{z) is then 
given by giiz) = (m|(z — K)~^\m), where K denotes the single-particle version of 
the standard tight-binding Hamiltonian (|3.2p . and we suppress the index a from now 
on since the leads are symmetric in the spin degree of freedom. Because the lead is 
semi-infinite and homogeneous, adding one additional site at the right end with the 
same hopping amplitude r will not change the Green function at the last site, that is 



= E E E 

o-'.o-" s'.s"=L.Rl'.l" 
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{m + l\{z — K) + 1) = gL{z) = giiz)- If we now apply the projection (|3.1ip with 
P projecting on the site m + 1 and Q on the rest of the lead, we obtain 

[gUzT^ = {m + l\iz-K)\m + l) 

= {m+l\{z - K)\m + 1) - {m + l\K\m){m\{zQ - QKQ^'^lm) {m\K\m + 1) 

= z + fiL- r^giiz) = [gL{z)]~^. 

In the first line we could interchange the order of calculating the 'inverse' and evaluating 
the scalar product since in this case P is projecting on a one-dimensional space. The 
complex roots of the resulting quadratic equation, 



T^gliz) - (z + fiL)9L{z) + I = 0, 



are given by 



+ i^ J^r'^ -{z + J^ ] if Im(z) > 
127^[^ + /"l + V4t2-(z + Ml)2] iflm(z)<0. 

The sign was determined such that the imaginary part of giiz) has a branch cut at the 
real axis and that lim^^-too = holds. The corresponding density of states at 

the last lattice site reads 

- Imgiiuj + tO) = - — 



PleadA^) = -- ^^9l{oJ + iO) = 2^4r2 - {uj + . (3.18) 



Later on, we will always take /9iead,L('^) to be energy-independent, i.e. perform the 
so-called wide-band limit. 



Generalisation to the Interacting Problem 

Up to now, we have only developed the projection method for a completely noninter- 
acting problem (arising from (13. ip by setting Uil" =0). In order to set up a finite set 
of fRG flow equations, this is all that we need, since these equations include only the 
free propagator evaluated at dot indices. To calculate the conductance of an interacting 
system we will, however, also need to project out the leads when the interaction terms 
are present, that is in the full propagator. 

Fortunately, it is easy to see that (|3.15|) and (|3.16|) also hold in the interacting case. 
Therefore we expand the full propagator as usual, 

where we have suppressed the frequency dependence. The summation extends only over 
the interaction region, since all self-energy diagrams vanish if one index is taken from 
the leads. If we now plug in the results from projecting the free propagator, we recover 
the noninteracting expression, 

^/To. = - E(^/-/')°(*^')*5^ +EE^Mi7i(^i,^;^2,cT)(gr,^,)°(i?')*5s 

(3.19) 

I' 
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likewise for Gq^.i, and 

l',l" I" hh 



I'V' 



(^?;/")° - ^Qi'-hli{h,(y\l2,(y){Qi^.i„f 



(3.20) 



^Qi'-i"{ti")*^v9LgR- 



V,l" 

3.2.2 The Conductance for an Interacting System 

Transport through a noninteracting one-dimensional system can be described by the 
Landauer-Biittiker formalism [Bruus&Flensberg 2004], but since we are dealing with a 
system of interacting fermions, it cannot be applied here. One could now argue that 
since we have set up an approximation scheme where the self-energy remains frequency- 
independent during the fRG flow, we map our interacting system to a noninteracting 
one, so that at the end all properties of the system (such as the conductance) can 
be calculated as in the well-known noninteracting model. For reasons of consistency 
it is, however, better to derive expressions for the desired quantities within the full 
interacting problem and then to verify that in our case these expressions correspond 
to the noninteracting ones (with renormalized parameters) without applying further 
approximations. Basically, this has been done before [Enss 2005a], but here we will 
present a slightly more general approach that allows for arbitrary couplings of the dots 
region with the noninteracting leads. 

Transport in Linear Response - Kubo Formula 

If we apply a voltage V = fiR — fii between the 'ends' of our two leads, we expect a 
current to flow. The linear conductance can in general be computed as [Oguri 2001, 
Bruus&Flensberg 2004] 

C=4lin/-»-^'°', (3.21) 

where the retarded current-current correlation function in frequency space K{z) is given 
by the analytical continuation of 

.13 



K{iuj) = [ dTe'^^{TjR{T)JLm- (3-22) 
Jo 



Here w is an even (bosonic) Matsubara frequency. Jl and Jr are the current operators 
at the left and right ends of the system. They are defined straight-forwardly as 

Jl ■■ = -Ml = ih[J\fL, H] = ih[ML, i^lead-dot] 



cr I 

(3.23) 

Jr ■■=Mr = -ihlAfn, H] = -ih{J\fR, -fflead-dot] 
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with Ml,r being the particle number operator of the left or right lead, respectively. 
Calculation of 

We will begin the calculation of the correlation function (I3.22p by expanding it into 
an infinite perturbation series. Due to the linked-cluster theorem, only connected di- 
agrams will contribute. Importantly, it is not necessary that all four external legs are 
connected to an overall connected diagram. Put differently: we can split up the four- 
point function into two parts, the first one {K'^) comprising all diagrams that consist of 
two unconnected parts each with two external legs, and the second one (A'"^) describing 
the overall connected diagrams (see also [Negele&Orland 1988, Eqn. (2.158b)]). The 
contribution of the former to (|3.22l) reads 



^'(^^)^=EE/ dTe'-^{- tX ('^rf/',.'(0)c^^(r))(Td,,,(r)c;,^,_^(0)) 

+ (tf)*t^ {Tdi,y{0)dl{T)){Tco,.,n{r)cl^,Lm 
- (tf)*(t^)* (Tco,.^L(0)<,(r))(rco,.,H,(r)4,,,(0))}. 

Those terms where one creation is paired with one annihilation operator at equal times 
vanish. Since the Hamiltonian is time translational invariant, they do not depend on 
time at all, and the integral dr exp{iuJT) with oj being an even Matsubara frequency 
is zero. If we now plug in the Fourier expansion of the Green functions and carry out the 
T-integration (which yields the inverse temperature times a delta function) as well as 
one spin summation (which is trivial because we assume spin conservation), we obtain 

Next, we express the correlation function in terms of Green functions with indices only 
in the interaction region by applying the projection technique for the full propagator 

K^{iLo) = T ^[-tft^(tf)*(t^,)*5fl(ie)9L(^u; + ie) G'{,.^^{ie)gi^,{iu: + ie) 

^'"'"''''^ " +tf{tf,rt^,{t^,rgn{^e)gL{^e) Gl.^,,{te)Gt^,{iu + ze) 

+ {tfTtvtk{twT9R{i^ + i^)9L{iu^ + ie)Ql^i{ie)gi.k'{i^ + ie) 
- {tf'ntlirtj^t^gniiu + ie)gL{ie) QlAie)Ql,,,(iu^ + ie) 



l,l',k,k';a ie 

(3.24) 
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(3.25) 



where we have defined 

fL{ie,ie + iuj) := giiie + ito) - giiie) 
fniie, ie + iuj) := gniie) - gniie + iuj) 
f{i€, ie + iuj) := fiiie, ie + iuj)fR{ie, ie + iuj). 

In order to carry out the ze-summation, we rewrite (|3.24|) complex contour 

integral with an appropriate closed integration path C by substituting T Fii^) — > 
— {2TTi)~^J^^dzF{z)g{z), where the simple poles at z = ie of g{z) should not coincide 
with those of F{z).'^ Since e is an odd Matsubara frequency in contrast to uj, the Fermi 
function fullfills this condition. Hence we choose g{e) = l/(exp(/3e) + 1) from now on.^ 
In our case, F{z) has a branch cut at the real axis Im(z) = as well as at the line 
lm{z) = —uj (because the propagators are nonanalytic at z = iui), so that the contour 
C consists of three parts. The Fermi function falls off exponentially for large positive 
arguments while it tends to one for large negative arguments. But since each propagator 
falls off linearly, we can deform the integration path into four lines slightly above and 
below the axes of nonanalyticity. This yields 

K'^{iu:)=Y, tXitkYi^r Jde^{ f{e + iO, e + zu;)gr,.^,{e + iO)g^^,,{e + iu,) 

-f{e-iO,e + zu:)gr,,,{e-zO)gr^,,{e + zu;) 
+/(e -iu;,e + 20)g,T.fc(e - iuj)gr.f,,ie + iO) 

-f{e -iLO,e- iO)g?,k{e - iuj)g^^k'(.^ " ^0) }. 

(3.26) 

The first and fourth term are of 0(a;^), and hence they will vanish in the limit uj ^ 
in (I3.2ip . The former can be seen if we consider 

1 



gs{e±iS) = ^^e±i6 + i\J ^t'^ - (e ± i5 + ^)2 



1 

2^2 



with [i being the chemical potential of both the left and right lead, which is now equal 
since we are calculating the conductance in linear response. Thus, for small uo 

/L,R(e±iu;,e±iO) = 0(a;) 

i 
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fL{e 
R 



fL{e 
R 



iO, e + iuj) 
ico, e + iO) 



=F- 



=F- 



i 



(e + fif + Oiio) 



V4r2-(e + ^)2 + 0(c^). 



Performing the analytical continuation iuJ 
second term of (|3.26p . we obtain 

g{e + uj 



l,l',k,k';a 



de 



a; + iO and substituting e ^ e + u; in the 
9{e) 



2TTi 



{/(e -i^,e + LO + iO)g,T.fc(e - iO)a,V(e + a; + iO)} + 0{uj'' 



(3.27) 



*This is the so-called Possion summation formula, for a complete review see e.g. [Mahan 2000]. 
^Be careful not to confuse this with the propagator gs{z). 
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The conductance that follows from this part of the correlation function then takes the 
simple form 

= [de g'{e)f{e -iO,e + iO) ^ iftz^(tf - iO)g?,k'{e + iO) 



-"^Y. [ 5'(e)|27rpiead(e)5;tf(t^)*a,V(e + iO) 

<T LI' 



(3.28) 



where we have used that g[it{z) = Qf,.i{z) due to time-reversal symmetry of the Hamil- 
tonian,^ and Q{ijJ + id) = {Qiuj — i6))* , and we have reintroduced the density of states 
at the last site pieadi^)- 



Vertex Corrections 



Before further commenting on (|3.28ll . we have to calculate the second contribution to 
()3.21|) which we will call K^. It arises from those diagrams where all four external 
legs are connected to an overall connected part, i.e. the connected two-particle Green 
functions, which can be related to the two-particle vertex functions via (|2.22|) . This 
will prove meaningful since we are directly calculating this quantity within our fRG 
approach. If we then carry out the r integration in (|3.22l) (which kills one frequency 
summation) and use that the two-particle vertex is frequency-conserving (which kills 
another frequency summation), we obtain 



E E E h'<;0.(^e + ^^)Or'',l,i^^) - {tf'TQCA^e + iu^)Qi^,,^(ie) 

1,1' li,...,l4, ie,ie' 

(3.29) 

where the li summations run only over the interaction region (since the two-particle 
vertex vanishes outside), and we have defined 

^ri}2;«3,u(^'^l' ^"^2; ^^3, iui4) := 72 {{h,cr', iuJi}, {^2, o", 1(^2}; {k, o"', ^"^3}, {^4, cr, ^^4}) . 



^Imagine the exact propagator Qi-ii expanded in an infinite perturbation series. If we now in each 
diagram interchange the direction of every internal line (which is possible if we assume that the single- 
particle dispersion is symmetric) and shift the time by r ^ — r (which is possible because the Hamil- 
tonian was assumed to be time-reversal invariant), we will end up with a diagram of precisely the same 
structure that would appear in the expansion of Qi'-i. 



40 CHAPTER 3. QUANTUM DOTS 



Applying the projection technique (|3.19l I3.20p yields 

a-,a' 1,1' ,k,k' li,...,l4 ie.ie' 

= ^ E E E E ihtkYVC^kiie + ioj)gr,[i^iie)fLiie, ie + iu;)Alli^iie; iu) 

a' I' ,k li,l2 it 

(3.30) 

with the definition 

(3.31) 

Up to now, we have not applied any approximations except for the fact that we are 
describing the transport through our interacting system in linear response. This means 
that ()3.28|) and (|3.30ll give an exact expression for the conductance, but in general it 
is impossible to evaluate at least the second one exactly even if we knew the exact 
propagator. In contrast, in the context of our fRG approximation scheme this becomes 
very simple which can be seen by an argument following [Enss 2005a]. Since we have 
chosen a truncation scheme that keeps the two-particle vertex frequency-independent, 
it is easy to apply Poisson's summation formula again in complete analogy to (|3.26|) to 
write down 

i,fc/;Z3,Z4 

gt.,^{e' + iu;)Ql,A^' + iO)fR{e' + iO, e' + ioo) 
-gr,,^{e' + iu;)ar,;fc'(e' - iO)/^(e' - zO, e' + iu) 
+gr.i^ie' + iO)g^^.,,{e' - iu;)fR{e' - iu;, e' + iO) 

-o?,i,{e' - io)grM'' - ^^)fR(^' - ^ - ^0) }• 

As before, in the first and the fourth term the function //j has frequency arguments at 
the same side of the branch cut and hence these terms are by an order of uj smaller than 



3.2. THEORETICAL APPROACH 41 



de' ( 
2^1 



the other two. Performing the analytical continuation we can thus write 

a 

[g{e' + u;) - 5(e')]^M3 i^' + ^ + i^)Qu A^' " iO)/R(e' - iO, e' + w + iO)} 
= 0M. 

This expression is independent of ie, hence we can perform the frequency summation 
in (|3.3QI) as well. Since the structure of this equation is precisely the same, it will also 
be of 0{uj), so that for the vertex corrections to the conductance we obtain 

= ^ lim ^ = 0. (3.32) 

Hence (|3.28|) is the expression for the conductance consistent with our approximation 
scheme. In the zero temperature limit, this result is much more general. Under quite 
weak assumptions for the two-particle vertex it is possible to show that the vertex 
corrections exactly vanish [Oguri 2001]. 



Special Cases 



We will now discuss the form of (I3.28P in some special cases. In this thesis, we will 
mainly focus on the zero temperature limit, in which minus the derivative of the Fermi 
function becomes a (5-function, and the conductance reads 



(3.33) 



If only two dots from the interaction region (1 and N) are connected to the left or right 
lead, respectively, we recover the familiar result 



The general form of (|3.28p naturally allows for defining the partial conductance of the 
spin up and down electrons, 



G — G^+Gi 



1 



27rpiead(0) j;tf(t^)*aj,,(0) + 2^piead(0)J]tf(t^)*a^,,(0) 



LI' 



LI' 



The name 'partial conductance' is justified because the so-defined expression for G^ 
would follow if we calculate the current-current response function for the spin up and 
down electrons separately by replacing J — > J^. := —Ma in (|3.22l) . 

There is one case where one can obtain a very simple expression for the conductance 
valid at arbitrary temperature, and we will note it since it is frequently used in the 
literature to describe transport even at T = where (I3.28P is exact. Namely, if the 
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interaction region comprises only one single dot (with interacting spin up and down 
electrons), the conductance can be expressed as [Meir&Wingreen 1992] 

2 2 4^2 ^2 [■ 

G = -^TTt7rE / 5'MpM/5.(u^), (3.34) 

with ti and tR being the couplings to the left and right lead, g{oj) the Fermi function, 
and p{ijj) and Pa{^^) the density of states at the last site of the leads and at the dot, 
respectively. 



Connection to Scattering Theory 

Everybody who is a bit confused by the form of the conductance if an arbitrary number 
of dots is connected to the leads or by its derivation should consider the following. In 
the noninteracting case (where (|3.28p is exact) one usually defines the zero temperature 
conductance using ordinary single-particle scattering theory as a factor of e'^/h times 
the absolute square of the transmission amplitude.^ The latter is most easily derived 
by looking at an eigenstate of the isolated left-lead Hamiltonian, 



{m\k, L) 



- sin [k{m — M)] for m from the left lead 
otherwise. 



M is implicitly assumed to be chosen such that this linear combination of left- and 
right-moving waves vanishes at an imaginary site added to the right end of the lead. 
Since everything is assumed to be diagonal in spin space anyway, we suppress the a 
index. Scattering states are the defined as usual, 

\k,L,±) := lim{±ir])g{ek ± iri)\k, L) . 

The free propagator Q of the system fullfills the Dyson equation, 

G = G + GHcoupG = G + GHcoupG, 

with G being the propagator of the full (noninteracting) system excluding the connection 
-f^coup between the leads and the dot region, G{z) = {z — H + -ffcoup)"^- Thus, for any 
m from the right lead we can write 



ijn\k, L, +) ={m 



^0 ^/ ^rr 

+ G{ek + tO)Hcoup- 



efc - + Feoup + ^0 --'^ ' ^""^efc-/7 + i/coup + ^0 

{m\k,L) +{m\G{ek + iQ)Hcoup\k, L) 

=0 

{m\G{ek + iO)HcoupG{ek + iO)Hcoap\k, L), 



\k,L) 



^One should note, however, that although it frequently appears in the literature (see e.g. [Datta 
1995]), there is a lot of sloppiness buried in such a definition. An equally intuitive but much more 
stringent approach to derive the Landauer-Biittiker formula starts out from the time evolution of the 
current operator, limt-too{4'{t)\j\4>{t)) (implicitly weighted by the Fermi function), with \4>{0)) being the 
ground state of the isolated leads (see [Schonhammer 2005]). The subsequent computation, however, 
is crucially based on the single-particle nature of the Hamiltonian governing the system. 
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which by plugging in the single-particle version of the level-lead coupling Hamiltonian 



Hcoup = - Yl 5^[tf|0,s)ai+H.c.], 

s=L,R I 

simplifies to 

{m\k, L, +) = J2{tfytjf{m\g{ek + iO)|0, R){l\g{ek + mi'){0, L\k, L) 



2-n 



2sin(A:) 



Y.^tf,rtf{l\G{e, + tm') 



1,1' 



Since the matrix element of Q is just an outgoing wave, (m|^|0, R) = t ^ exp{ikm — iip), 
we can now read off the transmission probability r(e) as 



Tie) 



2sm{k) 



LV 



(3.35) 



If we furthermore introduce the density of states at the last site of the lead, 



Plead (e) 



1 



27rT2 



V47 



1 



2-KT 



\/l — cos^{k) 



sin(A;) 



TTT 



we precisely recover the expression (|3.33l) derived above, showing that in the noninter- 
acting case both approaches are equivalent. At nonzero temperatures, the distribution 
of the lead electrons is governed by the Fermi function, such that the obvious T > 
generalisation of the noninteracting scattering theory definition of the conductance is 
to introduce another weighted energy integration, which is again consistent with (I3.28p . 

Our fRG approximation scheme maps the general interacting problem to an effective 
noninteracting one, and it turned out that we can compute the conductance using the 
[/ = expressions with renormalized parameters. However, since in the interacting 
case we have to define this quantity in a different way (as we cannot use single-particle 
scattering theory), this is not obvious. 

From the scattering theory approach it is also immediately clear that since the trans- 
mission probability is bounded by one, the maximum conductance is given by e^//i per 
channel, that is 2e^ /h if take into account the spin, and /h in spin-polarised situations. 



Particle Numbers 

An important quantity that is also accessible in the experiment is the occupation number 
of each dot in the interaction region. For the dot with single-particle index / it is given 

by 

with 6 going to zero. For T ^ the sum can be written as an integral, 

«) = ^ /'^^ e^'^^iiu:), (3.36) 
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which is easy to tackle numerically, since the contribution from large uj can be computed 
similarly to (|2.43p . At nonzero temperatures we use Poisson's summation formula to 
obtain 

iu)n 

dz e^'Qli{z)g{z) 

(3.37) 

which is easy to calculate numerically as well. 




Chapter 4 

Numerical Results 



In this chapter we will present results from the application of the fRG to dot systems 
with a few levels for various geometries. Unless stated otherwise, we will always use the 
truncation scheme that includes the flow of the (frequency-independent) two-particle 
vertex. Despite that fact this approximation can strictly be justified only in the limit of 
small two-particle interactions, we will also apply it to fairly large U. Allover this chap- 
ter, we will only sporadically comment on what 'fairly large' means and on the question 
which parameter regions are for sure out of reach within our approach. We will focus 
here on describing the physics that arises due to the presence of the interaction with the 
implicit understanding that the fRG produces reliable results for all situations shown. 
In the next chapter, by comparison to other data either obtained from a (complicated) 
exact solution or numerical methods known to be very precise we will verify that this 
is indeed in the case. We will furthermore specify in more detail which interaction 
strength can no longer be tackled by our simple fRG truncation scheme. 

Since our main interest is in transport properties of dot systems, we will compute the 
linear-response conductance G as well as the transmission phase a (at zero temperature, 
this is the argument of the transmission probability) as a function of the gate voltage Vg 
that shifts the single-particle energies of each dot. For situations with spin degeneracy, in 
particular in absence of magnetic fields, the transmission phase of the spin up and down 
electrons is equal and will be denoted as a := = a^. Frequently, we will also calculate 
the occupancy of each level since this often facilitates the physical interpretation of the 
results. Unless stated otherwise (that is everywhere outside the section called 'finite 
temperatures'), we will focus on the low-energy physics in the T = limit. 

Every system under consideration exhibits three typical energy scales that determine 
the gate voltage dependence of the conductance, namely a single-particle energy A 
(which might be a level detuning between parallel dots or a nearest-neighbour hopping 
between sites of a chain), the two-particle interaction [/, and the hybridisation strength 
with the leads F. To systematically study their influence on the physics that governs 
G{Vg), we will pursue the following course of action which is guided by the limitations 
imposed by the perturbative nature of the fRG. First, we will discuss the noninteracting 
limit [/ = 0, varying the single-particle spacing from A <C F to A ^> F. Next, we will 
turn on an interaction such that the physics is significantly influenced (which will mostly 
turn out to be [/ F). If possible, we will finally increase U such that all important 
limits can be captured ({A <^r,U ^ F}, . . .), allowing us to determine which quantity 
governs the physics in each case. This will, however, only turn out to be possible in the 
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spinless two-level case and for the SIAM. Hence it is meaningful to focus on the effects of 
gradually turned on small interactions rather than concentrating on the extreme limits 
where one quantity is much larger (smaller) than the others. 

It is very important to point out that everything shown in this chapter represents 
the generic behaviour of each particular system under consideration, and, even more, 
we will refrain from discussing special non-generic situations. In order to obtain a 
complete picture of what should be associated with the former and what with the 
latter, we have taken advantage of the small computational power required by fRG to 
scan wide regions of the parameter space. The consequences showing up if one carelessly 
sticks to situations of too high symmetry, which are a formidable candidate for yielding 
non-generic behaviour, can be quite drastic. In [SOG 2002], a noninteracting spinless 
parallel double dot was studied, and the authors focused on completely symmetric 
hybridisations. The paper aimed at giving a clue towards an explanation of the evolution 
of the transmission phase observed in quantum dot experiments (an issue that we will 
also tackle later on). Though giving many important insights, the results would have 
been more enlightening if generic hybridisations had been considered, as later on stated 
by the authors themselves [Golosov&Gefen 2006]. 

From now on, we set the chemical potential of the leads to zero. Furthermore, we 
perform the wide-band limit, which is achieved by substituting r — > r/r for the hopping 
matrix element of the leads and tf — > y^tf for all hoppings from the last site of the 
latter into the dot region and taking the limit 77 — > oo. The propagator 5^ (^a;) then 
reads 

-■ir"^ Im(u;) > 
+iT^'^ lm{uj) < 0, 

and the density of states piead becomes energy-independent, which justifies the name 
wide-band limit. It is important to note that we do not perform this approximation 
for computational reasons, but only because it usually appears in the literature. It is 
motivated by the fact that the details of the leads should not influence the transport 
properties of the system dramatically. For the system of two parallel spin-polarised dots 
(Sec. I4.1.2P we have verified that this is indeed the case. Since the density of states is 
independent of energy, the same will hold for the hybridisations 

Tf := vrltf iVlead. (4.2) 
Allover this chapter, the unit of energy is chosen to be^ 

^■■=T.^r.= E rf, (4.3) 

/ s=L,R 

and we introduce the shorthand notation F = {F^/F F^/F F^/F . . .}. 



4.1 Spin-Polarised Dots 

In this section we will apply the fRG to quantum dot systems where the spin degree 
of freedom is ignored. As mentioned above, the Hamiltonian to describe spin-polarised 

^Choosing F as the unit of energy might be questionable, especially if the system under consideration 
comprises a large number of dots A*' where one would expect the physics to be governed by the ratio 
U/Ftyp and A/Ftyp (with Ftyp :— ^/N) rather than by U/F and A/F. Practically it turns out, however, 
that for all situations considered here F is indeed the most suitable choice of the unit of energy. 
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systems is obtained by dropping all spin indices in the general Hamiltonian (I3.1|l . We 
will focus on systems of up to six parallel dots each containing one level (or one dot 
containing up to six levels, or mixtures of both situations; we will keep the aforemen- 
tioned linguistic sloppiness and do not distinguish between 'dots' and 'levels', ignoring 
the spatial geometry that is induced by the choice of the hoppings t;,/')- The term 
'parallel' means that every dot is connected to the leads by nonzero tunnelling barriers 



On the one hand, the physical importance of a model neglecting the spin degree of 
freedom arises as it should reproduce the results of a model containing spin if a large 
magnetic field is applied. In the next section, we will exemplarily demonstrate that this 
is indeed true. On the other hand if it occurs in the experiment that the spin degree of 
freedom does not seem to play a role (signalised by the absence of Kondo physics), it 
might be justified to try to explain these experiments using a spin-polarised model. 

This section is organised as follows. First, we will as an introduction describe transport 
through a single impurity in a homogeneous chain, which is by definition a noninter- 
acting model and therefore exactly solvable. It will serve for a better understanding of 
transport through more than one dot if the singe-particle level spacing is assumed to 
be large so that the different levels do not overlap. This picture of transport occurring 
through each level individually will remain valid in presence of a two-particle interac- 
tion, only that the separation of the transmission resonances is enlarged due to Coulomb 
repulsion. In contrast, the effect of the interaction will be much more dramatic if the 
spacing between the levels decreases so that they overlap significantly. In a noninter- 
acting picture, this would imply that they simultaneously contribute to the transport 
so that no well-separated resonances are to be expected. In presence of interactions, 
we will observe such peaks of good separation nevertheless due to Coulomb repulsion 
('Coulomb blockade peaks'). Furthermore, we will frequently find the GiVg) curve to 
exhibits additional correlation induced resonances (CIRs) if the integration between the 
electrons exceeds a certain critical value depending on the dot parameters. These CIRs 
were first predicted for a two-level dot by [Meden&Marquardt 2006] using the fRG. 

As stated in the introduction, our main focus within this section will be on the number 
of resonance peaks and transmission zeros (and corresponding jumps in the phase a) 
in the conductance G{Vg) as a function of the dot's single-particle level spacing. If 
an experimentally realised quantum dot is successively filled with electrons, one would 
expect the latter to decrease for those electrons that occupy the states of highest energy. 
Taking this as an input from the experiment, our results might serve as an explanation of 
the puzzling behaviour of the phase evolution a{Vg) experimentally observed in quantum 
dots [YHMS 1995,SBHMUS 1997,AHZMU 2005]. 

4.1.1 A Single Impurity 

For a single dot coupled to tight-binding leads the Hamiltonian (|3.ip reads 



t: 



s=L,R 




The propagator projected on the dot (13. lip is then given by 



g{iLo) 



1 



(4.4) 



iu — Vg + i sgn(a;)r 
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and the linear-response conductance reads 



G{Vg) 



h + r2 ■ 



(4.5) 



0. 



2 r r 

This is a Lorentzian of full width 2T and height 4^ ^^^^^^ centred around Vg 
In the left-right symmetric case (Tl = ^r) the conductance reaches the unitary limit, 
and the more asymmetric the hybridisations are chosen, the more it is suppressed. The 
transmission phase a, which can be computed as 



a{Vg) 



arctan ( — 

V a 



(4.6) 



changes by vr over the resonance. In fact, for a single dot at T = it is no independent 
quantity in the sense that it is directly connected to the conductance by a generalised 
Friedel sum rule [Hewson 1993]. The same holds for the level occupancy, which equals 
the transmission phase in units of tt. The conductance and the transmission phase are 
depicted in Fig. 14.11 



4.1.2 Parallel Double Dots 
Application of the fRG 

We will now turn to a more complex system of two parallel spin-polarised dots. The 
free propagator with the leads projected out can be cast in the form 



lUJ 



^9 + T + ^ sgn(u;) Ta i sgn(u;) 



\ 



Vn 



+ i sgn(cj) 



/ 

(4.7) 

where we have introduced the level spacing A, and a vanishing direct hopping between 
both levels was assumed, since it can always be tuned away by a basis transformation of 
the dot states (leading to an additional level detuning). The notation used to derive the 
general expressions for the effective projected Hamiltonian (13. lip and the conductance 
(|3.28|) allows for all hopping matrix elements being arbitrary complex numbers. This is 
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necessary if we for example want to describe the effect of a magnetic flux piercing the 
ring geometry. Here we will concentrate on real hoppings, in which case the off-diagonal 
element of the free propagator can be expressed by the square roots of the hybridisations 
and the relative sign s of all hopping matrix elements, s = sgn(t^t^t^t^). The fact 
that it is only the relative sign that enters the calculation of all quantities of interest 
is easily seen. Assuming an arbitrary choice for sf := sgn{tf) and flipping <-> 
does not change the projected part of the free propagator (13. lip or the expression for 
the conductance (|3.33|) . The interchange of with changes the overall sign of the 
off-diagonal elements of Q^, which implies that also the sign of the off-diagonal elements 
of the full propagator is flipped, while that of the diagonal terms is unaffected. This 
becomes clear if we expand the full propagator in an infinite perturbation series. Here 
this series only contains one single type of interaction vertices {va,b,A,b), and hence 
each term in the expansion of Qi-i (Gi-i) has to contain an even (odd) number of factors 
(Qij)^- However, in the expression of the conductance every off-diagonal term acquires 
an additional minus sign as well, and all that enters the computation of the average level 
occupancies are the diagonal parts of the propagator. Hence all quantities of interest 
depend only on the relative sign of all couplings as claimed. By writing down (14. 7|) we 
have implicitly assumed that this relative sign is realised such that tf" > 0, > 0, and 
sgn(if ) = s. 

The fRG flow equations (|2.40p for the effective level position V/=a,b(A) := 
and inter-level hopping t{A) := ■ji{A; B) read 



dAViiA) 
dAt{A) 



LU = ±A 



mi 

2tt 



E ^a;b(«w;A), 



(4.8) 



=±A 



with / denoting the complement, A = B, and the flow of the effective interaction 
[/(A) := 72 (A, B; A, B) (CTIl is given by 



dKU{k) 



2tt 



[ - ^AiAli"^; A)^B;B(-ia;; A) -h^A;B(«w;A)^B;A(-iu;;A) 

+ GA;B{i(^] A)^B;A(2W; A) - ^yl;A(«W; A)g B;B{ii^] A) 

(4.9) 



a;=±A 



We have chosen the matrix Q as 

-1 



e(iw;A) 



iuo - Va{A) + i sgn(u;) Ta 
t{A) + i sgn(u; 



t{A) + i sgn(LL' 

iio - Vb(A) -I- i sgn(u;) Tb 



so that the initial conditions (|2.44l) read Va{A —>■ oo) = — ^, yB(A — > oo) = Vg + ^, 
t{A oo) = 0, and U{A — > oo) = U, with U being the interaction between the electrons 
of both levels. Solving these differential equations provides us with the full propagator 
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G = G{A = 0) of the system and we can compute the zero temperature conductance 
using (|3.28l) . which here takes the form 



e2 



(4.10) 

Again, we have assumed that the relative sign of all couplings is realised in a way that 
tf > 0, > 0, and sgn(t^) = s. 

The Noninteracting Case 

Since the double dot system is already a quite complex geometry (in the sense that it 
has a fairly large number of parameters to be varied), we will first study it without 
interaction. Replacing Q by the free propagator (|4.7|) in (|4.10p then yields 



GiVg)=4-- , T2 ? T2' (4-11) 



+ 



with the obvious definition €a = Vg — ^ and cb = Vg + ^ ■ As mentioned above, we will 
always focus on studying the generic behaviour of the conductance and transmission 
phase as a function of the gate voltage. The term generic refers to that fact that we 
would expect the qualitative features and dependence on the dot parameters U, A and 
s of both curves to be independent of the choice of the hybridisations. This is indeed 
true for all except for a three-dimensional manifold in the four- dimensional space 
of all couplings. In particular, it will prove that a left-right and A-B asymmetry only 
influences the overall height of the conductance curve, the set-in of the crossover regime 
at scale Across, and the correlation induced structures which vanish for A > Acir (see 
below). An indication for this is already given by the noninteracting case if one analyses 
the exact expression (|4.11|) . For the interacting case, we have verified numerically by 
scanning a large part of the parameter space that the results presented in this section 
are indeed the generic ones. Many parameter sets {T^} included in the manifold of non- 
generic hybridisations are not of physical relevance anyway. For example, it shows that 
the qualitative behaviour of GiVg) for degenerate levels and s = -|- is not generic in the 
left-right symmetric case. For a more detailed discussion of the manifold of non-generic 
hybridisations see [Meden&Marquardt 2006]. 

For degenerate levels the curve described by (|4.1ip is a Lorentzian of full width 2r 
for large gate voltages that has a dip and a transmission zero located at = 0. For 
a large A-B asymmetry (that is for Ta/^b ^ 1, or vice versa), this can be viewed as 
a Fano anti-resonance resulting from destructive interference between the two different 
path through the system. The transmission phase changes by vr over each 'peak' and 
shows a phase jump corresponding to the transmission zero in between. Increasing Vg 
both levels are heightened in energy and therefore successively depleted. 

For s = + increasing A shifts the transmission peaks to larger Vg, the transmission 
zero always remaining in between. For large A ^> F their separation is given by A, 
while the width and height of each peak is determined by the individual hybridisations 
F^'^ and F^'^. This is more or less obvious since in the limit of large A, the domi- 
nant contributions to the conductance (I4.10p will come from the diagonal parts of the 
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FIG. 4.2: Gate voltage dependence of the conductance G (black) and transmission phase a (red) 
of a noninteracting two-level dot with generic level-lead hybridisations F = {0.1 0.3 0.4 0.2} for 
one small and one large level spacing. Level occupancies for dot A (green) and dot B (blue) are 
shown as well. 



propagator which have support on Fg-regions separated by A. The transmission phase 
changes by n over each peak (and jumps in between), and again the form of a{Vg) is 
the same as in the single-level case. Both dots are depleted individually over their asso- 
ciated resonances, which finally manifests the intuitive picture that for nondegenerate 
levels with level spacing larger than the level broadening transport is observed each time 
the energy of one level crosses the Fermi energy of the leads. For s = — the crossover 
regime from A ^ F to A ^ F is very different. One of the peaks present if the levels 
are degenerate splits up while the other becomes vanishingly small. For large A we 
recover the resonances corresponding to the individual levels, but now the transmission 
zero lies outside next to the (not observable) third resonance. Fig. 14.21 again shows the 
described behaviour for the two regimes A <C F and A S> F for s = -|- and s = —. 

Interacting Case, Nearly Degenerate Levels 

If we add an interaction U of order F between the electrons on both levels, the physics 
changes. For nearly degenerate levels A <C F increasing U enlarges the separation of 
the two peaks (the wings of the Lorentzian with the dip at Vg ^ 0) and changes their 
lineshape. The height of the peaks increases as well as their separation which is given 
by U and hence comparable to their width which is of order F (for [/ w F < Uc). At a 
certain critical strength of the interaction (depending on the level-lead hybridisations) 
C/c({rf }i s) = 0(T) each peak splits up into two, the height of all four resonances being 
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FIG. 4.3: Gate voltage dependence of the conductance G of a two-level dot with generic level- 
lead hybridisations T = {0.1 0.3 0.4 0.2} at different local interactions. If the latter exceeds 
a critical strength depending on the dot parameters, additional sharp resonances split of the 
Coulomb blockade peaks located at = ztU/2. Note that the curves for the largest U around 
= are only shown in the insets. 



equal and now independent of U. For large U = 0(r), the outer peaks are located 
approximately at Vg ^ /2, and their width is again of order T and not of order U 
or Ti. They are frequently interpreted in a usual Coulomb blockade picture. The inner 
resonances result from a subtle many-particle effect (justifying the name correlation 
induced resonances, CIRs) and cannot be explained in a simple way. Only recently they 
were predicted by [Meden&Marquardt 2006] using the same fRG approach we employ 
here. Increasing U , the CIRs become very sharp and are located exponentially close to 
Vg = 0. All this is shown in Fig. I43l 

The critical Uc that determines the interaction strength at which the additional res- 
onances show up depends on the hybridisations Tf as well as on the relative sign of the 
level-lead couplings s. If we fix F as the unit of energy, we have to analyse a three- 
dimensional parameters space for each s = -|-, — in order to obtain a complete picture 
of Uc that might serve to better understand the nature of these structures. Since this is 
optically impossible, we stick to two-dimensional manifolds, and we start with s = —. 
One important choice for such a manifold is certainly that where the total coupling 
strength of dots A and B (governed by F^/F) as well as the left-right asymmetry as- 
sumed to be equal for both levels, F^/F^ = F|^/F^ = Tr/Tl, is varied (Fig. [13 left 
panel). We find that Uc is always roughly of order F, and it decreases for decreasing 
symmetry [A-B or left-right) of the system. Another important parameter manifold 
arises for fixed Ta/Tb but arbitrary F^F^ and FyFf (Fig. ri 

ght panel). If we 

assume dot B being stronger coupled than dot A, we find that for constant F^/F^, Uc 
decreases the more Fg/F^ deviates from one. Furthermore, Uc becomes larger (smaller) 
if for fixed Fg/F^ < 1 (Fg/F^ > 1) F^/F^ is increased. In the s = -h case, an equal 
left-right asymmetry does not represent generic parameters; computing Uc as a function 
of Ta/^b and F^/F^ = 2F|j/F^ yields the same picture as in the s = — case with 
fixed asymmetry. We will refrain from a more detailed discussion here. 

The transmission phase a changes by vr over each Coulomb blockade peak and jumps 
by vr at the transmission zero which remains close to = for arbitrary strength of 
the interaction. If present, the phase changes steeply when crossing the CIRs (see the 
upper left panel of Fig. I4.5p . but basically it exhibits an S-like lineshape that becomes 
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FIG. 4.4: The critical Uc that determines the appearance of the CIRs in the conductance of a 
double dot with degenerate levels and s = — for fixed left-right asymmetry F^/F^ = T^/F^ = 
^rI^l (left panel) and constant A-B asymmetry F^/F = 0.25 (right panel). 

more pronounced the larger the interaction is chosen (compare the upper left panel 
of Fig. 14.61 with U/T = 10.0 in Fig. 14. 8|) . The level occupancies experience a more 
dramatic change if the interaction is switched on. In the limit of large U, they no longer 
depend monotonically on the level energy. Coming from large negative Vg, the dot that 
is more strongly coupled gets depleted over the first Coulomb blockade peak, while the 
occupation of the other stays close to one. At Vg w 0, the former is filled again while 
now the latter gets depleted, and a transmission zero instead of a resonance is observed. 
This inversion of the level populations if followed by another depletion of the more 
strongly coupled level over the second Coulomb blockade peak. In the limit of large 
Vg, both levels get depleted. It is important to note that the population inversion is 
not directly related to the appearance of the CIRs, since the strength of the interaction 
where the former sets in does not coincide with the critical Uc (see for example Figs. 
14.51 &: 1461) . Furthermore, one should be aware that the total occupancy still depends 
monotonically on the gate voltage and changes by one over each resonance. 

Interacting Case, Large Level Spacing 

If we increase the level spacing A, the CIRs gradually vanish. The largest A =: Acir 
where they are still observable is a function of the other parameters of the dot, but 
roughly it is given by A r/50. Applying a larger level detuning leaves only the 
Coulomb blockade peaks at position Vg ^ A ± U/2 which have width of order T and 
the transmission zero in between. The subsequent evolution from A < F to A S> F 
(the crossover regime) is completely similar to the noninteracting case. For s = + the 
separation of the peaks increases, but the conductance always vanishes at some point 
close to = 0. For s = — one of the peaks splits up (the one that one would associate 
with the more strongly coupled level^), while the other becomes vanishingly small. 

^In that fashion that one would think of the level with higher energy to get depleted first; however, 
the universal peak width of F shows that such an interpretation is incorrect. In fact, it would be best 
just to say that if dot A (B) is the more strongly coupled one, it is the right (left) peak that splits. 
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FIG. 4.5: Gate voltage dependence of the conductance G (black) and transmission phase a 
(red) of a two-level dot with generic level-lead hybridisations T = {0.1 0.3 0.4 0.2} for four 
different level spacings A. Level occupancies for dot A (green) and dot B (blue) are shown as 
well. For nearly degenerate levels, the phase changes rapidly because of the remnants of the 
CIRs. The effect of the population inversion is only weakly developed since the interaction is 
still too small. For large level spacings, the conductance is governed by transport through the 
individual levels. 



The scale Across that separates the A <C T regime where one observes two resonances 
with a transmission zero in between from the crossover regime where one of the peaks 
significantly decreases while the other begins to split depends on the parameters of the 
dot (and of course on its precise definition), but a very rough estimate should be Across = 
r/10. Surprisingly, Across is always of order T no matter how U is chosen (Fig. 14. 8| ). 
If the CIRs are present, it might not be possible to distinguish their vanishing and the 
splitting up of one peak if for certain dot parameters it happens that A^r ~ Across, 
as it is likely to happen if these structures are well-pronounced, that is for sufficiently 
large U. In both the s = + and the s = — case in the limit A ^> F one ends up with 
two peaks either with a transmission zero and associated phase jump of vr in between 
(s = -|-), or the conductance stays finite and the phase evolves continuously (s = —). 
As in the noninteracting case, this situation corresponds to separate transport through 
each of the original levels, but in contrast to U = 0, the separation of both peaks is now 
given by U + A. 

The latter follows from the intuitive picture that each time one dot gets depleted, all 
other levels are lowered in energy by U due to the lack of the Coulomb repulsion. A 
more precise understanding can be gained if one diagonalises the interacting isolated dot 
system decoupled from the leads, which is meaningful if the energies of the former {U 
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FIG. 4.6: The same as Fig. 14. 5^ but for s = — . Here, the interaction does not exceed the critical 
Uc so that the CIRs are not observed for nearly degenerate levels, while the population inversion 
is already visible. In the subsequent evolution, one peak splits up and the transmission zero 
moves outwards. 



and A) are much larger than the energy scale that is determined by the coupling to the 
latter. In the occupation number basis {|0, 0), |1, 0), |0, 1), |1, 1)} the dot Hamiltonian 
reads 



H, 



dot 



"9 2 

2 

V / 



One would now expect transport to be possible when the states with zero and one or 
those with one and two electrons have the same energy. At = — + A)/2 this holds 
for |0, 0) and the single-particle state with lowest energy, |1, 0), and at = (A + U)/2 
the latter and the two-particle state |1,1) are degenerate. These are approximately 
the positions of the resonance peaks for large level spacing. The same argumentation 
holds if the dot system comprises arbitrary many levels each with levels spacing A 
and two-particle interaction U but no direct hopping. For arbitrary Vg, the n + 1- 
particle state of lowest energy will then exceed the energy of that with n particles by 
Vg + nA + nil, and likewise the latter exceeds the energy of the n — 1-particle state by 
Vg + {n — 1)A + {n — 1)U . Hence the differences differ hy A + U as claimed. 

A final overview of the conductance and transmission phase of the double dot geometry 
that particularly accounts for the different energy scales within the system (though 
containing nothing that has not been mentioned before) is shown in Fig. 14.81 
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Effective Level Interpretation 

To gain more insight into the interacting case it is often useful to exploit that the self- 
energy computed by our fRG scheme is frequency-independent. Hence the effect of the 
electron correlations is buried in the renormalized level positions and hoppings (which 
are not present in the free propagator, but in general they are generated by the flow) 
as a function of Vg . Equivalently, one can consider the eigenvalues and eigenvectors 
V"^ of the dot system with these effective parameters decoupled from the leads. The 
eigenvectors {(Ai, A2), {fj-i, ^^2)} determine the hybridisations of new levels with energies 
^iiVg) to the leads by expressing the coupling part of the Hamiltonian in terms of the 
new eigenstates of the effective isolated dot Hamiltonian and reading of the hopping 
matrix elements. More precisely, the old dot operators dA,B can be written as 

dy = Ai4t)+^4t) 



hence the new couplings read 



^2 — + t^2tB^ 

and the hybridisations follow as 7^ = vr|t|p/3iead- Transport through this new system 
is then equivalent to transport to the original one, which is evident because the Green 
function ^1^1^ is of course independent of the choice of a basis in the dot region, and 
in an (effective) noninteracting problem it is irrelevant which matrix element of the 
propagator we use to calculate the conductance, so that Gi^-i^ determines the latter 
completely. If we have an arbitrary number of parallel dots, the new couplings follow 
similarly, 

N 

One should note that they are normalized because of the orthogonality of the V^, 

i,s i,s j,j' s,j 

In the case of nearly degenerate levels (Fig. 14.71 left panel) and regardless of the 
relative sign s, we find that for Vg — > =Foo the new effective levels (the eigenvalues 
of the isolated dot system at the end of the flow) are given hy tOi = Vg ± U/2. This is 
plausible since in the limit of large gate voltages, the latter is the dominant energy so 
that for — > 00 we would indeed expect tOi = Vg — U/2. If on the other hand both 
dots are filled {Vg — > —00) their energy is shifted by the constant ?7, which explains 
the observed behaviour Ui = Vg + U/2. In contrast, in the region around = 0, the 
interaction (assumed to be sufficiently large) has a more dramatic effect, in particular 
the eigenergies now no longer depend monotonically on the bare energy of the original 
levels. The more strongly coupled level (referring to the new couplings determined from 
the eigenvectors) crosses the Fermi energy of the leads at the first Coulomb blockade 
peak, but close to = both levels cross (and cross /i) so that a transmission zero 
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FIG. 4.7: Conductance G (black) of a two-level dot with generic level-lead hybridisations 
r = {0.1 0.3 0.4 0.2} for two different level spacings. For almost degenerate levels, the level 
occupancies (dot A: green, dot B: blue) depend non-monotonically on the gate voltage Vg. The 
same holds for the effective eigenenergies Wi (more weakly coupled dot: red, more strongly 
coupled dot: orange). Note that it is always the more strongly coupled effective level (original 
dot) that crosses the chemical potential (gets depleted) at each Coulomb peak. For large level 
spacings, LUi crosses /i linearly at each peak. The new hybridisations almost coincide with the 
original ones and are therefore not shown (this does not longer hold for A > 2 and A <C F 
where correlations enhance the asymmetry between the dots). 



(a type of Fano-antiresonance) instead of a peak is observed. At the second Coulomb 
resonance it is then again the more strongly coupled level that crosses the chemical 
potential. This is closely related to the aforementioned population inversion, since it 
turns out that the renormalized level positions are the dominant part of the effective 
Hamiltonian decoupled from the leads that one diagonalises in order to obtain the 
Wj, implying that the hybridisations barely change. The whole scenario of the most 
strongly coupled level crossing the chemical potential again and again at each peak is in 
close analogy to the one observed by [Silvestrov&Imry 2000] (SI) who studied chaotic 
quantum dots containing a few hundred levels. One should note, however, that the 
situation described here is more general since it does not require Si's assumption of 
strong asymmetry in the hybridisations. In particular, the behaviour of the level with 
largest 7^ crossing ^ at each resonance is observed generically even if not one of the 
original Tj is much larger than the other. 

For large level spacings (Fig. 14.71 right panel) it shows that the renormalized level 
positions are not only dominant (as for nearly degenerate levels) but at least two orders 
of magnitude larger than the hoppings generated by the flow, and therefore the former 
are the eigenenergies of the effective dot system with unchanged and s. At each 
Coulomb blockade peak, one of them (say wa) crosses the chemical potential. For all 
gate voltages except those close to the other Coulomb peak uja depends linearly on 
the gate voltage {uja = Vg + c). This explains why the form of the resonances and 
the behaviour of the transmission phase are precisely the same as in the noninteracting 
case. In particular, the height and width of the peaks are determined by the parameters 
of the corresponding level through which transport takes place, and a phase jump 
is observed in between for s = — . At the other Coulomb peak (when u!b crosses 
the chemical potential and the other dot is depleted), the energy oja is shifted by U 
due to the interaction of both electrons, so that in the limit Vg ^00 we recover 
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FIG. 4.8: Systematic account for the three different energy scales (nearest-neighbour interac- 
tion [/, level spacing A, and hybridisation strength F) that govern the conductance G (black) 
and transmission phase a (red) of a spinless two-level dot. The relative sign of the level-lead 
couplings is s = — for the solid and s = + for the dashed curves, and the hybridisations read 
F = {0.1 0.3 0.4 0.2}. Note that the depicted behaviour is the generic one and in particular 
qualitatively independent of the actual choice of F = {. . .}. The most fundamental points of this 
figure are the following. For nearly degenerate levels, we observe two transmission resonances 
(Coulomb blockade peaks) with width of order F separated by U . They are accompanied by 
additional correlation induced resonances for appropriate choices of U and A ^ F. The trans- 
mission phase shows an S-like change of tt over each peak and jumps by tt at the transmission 
zero located close to = 0. Increasing the level spacing, the behaviour of the curves begins to 
depend on s. In the Hmit A ^ F, we find two resonances of width F; separated by -I- A. a 
drops by tt over each of them and jumps by tt at the zero which is still in between for s = +, 
while evolving continuously for s = — . 
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uji = Vg±{U + A)/2. 

A word of warning about the interpretation of the system in an effective single-particle 
picture is in order. Since at the end of the fRG flow we are left with a frequency- 
independent self-energy it is of course possible to describe transport from such a one- 
particle point of view, or, even simpler, to diagonalise the effective Hamiltonian to 
obtain new level energies Ui with new hybridisations as a function of the gate voltage. 
However, it is important to note that in the first instance everything that can be ex- 
tracted from this picture is that if we had a noninteraction dot system with energies 
and corresponding level-lead couplings t|, then transport through such a system would 
coincide with transport through the effective interacting one and therefore allow for a 
more or less easy interpretation.^ But what we a priori do not know is why the in- 
teraction renormalizes the levels the way it does. The effective picture may therefore 
in some cases facilitate the interpretation of the results (as in the case of large level 
spacing, where the effective energies depend linearly on Vg and are only shifted by U 
if one level gets filled, whereas the new hybridisations coincide with the original ones 
and are in particular independent of Vg) while it proves to be of limited validity in oth- 
ers. This is most strikingly shown at nearly degenerate levels where the conductance 
exhibits Coulomb blockade peaks of width F. But since the hybridisations of the new 
levels basically coincide with the original ones (because the hopping generated by the 
interaction is small compared to renormalized level energy), it is impossible to interpret 
them as resonances arising from transport through separate levels. To put it the other 
way round, those peaks have to be caused by several overlapping effective levels, and it 
is of course far from evident why the interaction produces such a situation. 

4.1.3 Parallel Triple Dots 
General Considerations 

Next, we consider the parallel geometry that consists of three dots. The free propagator 
of this system can in complete analogy to the two-level case be cast in the form 

/iu)-Vg + /\ \ 
[G^iuj)] = iu;-Vg - ^eff + Hpp, (4.13) 

V iu:-Vg-Aj 

with the last term being the contribution from the projected leads. 



-ffeff — Hpp 
i sgn(u;) 




As usual, is a variable gate voltage that simultaneously shifts the on-site energies of 
the dots, A denotes a level spacing identical for all levels, and si = sgn(t^t^t^t^) and 

^One should note, however, that even for (/ = and large A no simple explanation (besides the one 
buried in the exact expression for the conductance) of why s governs the occurrence of a phase jump 
between the transmission resonances exists. 
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S2 = sgn{t^t^t^t^) the relative signs of the couplings, and the energy-independent 
hybridisations with the leads in the wide-band limit. As before, the relative signs are 
assumed to be realised such that = si, = S2, and > otherwise. Furthermore, 
we introduce the shorthand notation s = {S1S2}. The general expression for the zero 



temperature conductance (|3.33|) then reads 

=2 



G 



e 



2vrpiead(0) Vs/,/,Jrfr^^gi;Z,(0) , (4.14) 



with 



51 for I = B and arbitrary V 

52 for Z = C and arbitrary I' 
1 otherwise. 



The full propagator in presence of interactions between all electrons, 

U := VA,B,A,B = VA,C,A,C = Vb,C,B,C, 

is obtained by solving the flow equations for the self-energy and the effective interaction 
()2.4Ql I2.4ip . Alltogether, they comprise twelve equations, six for the effective on-site 
energies V;(A) and inter-dot hoppings ti^ii{K) (which are in general of course generated 
by the flow), and six for all the independent components of the two-particle vertex, 
72(i, j;A;,^; A) with {ijkl) = {ABAB,ABAC,ABBC,ACAC,ACBC,BCBC}. 

It is important to note that for otherwise arbitrary parameters relaxing the assumption 
of equal interactions between all electrons and introducing asymmetric level spacings 
only leads to quantitative changes but does not influence our results qualitatively. 

In contrast to the previously studied double dot geometry, we will refrain from first 
describing the noninteracting case here. This is because for small level spacings A <^ 
r the U = case is qualitatively very different from that with U > and hence 
no additional insight into the interacting problem can be gained from considering the 
noninteracting case. Therefore we will delay the discussion of the latter until the end 
of this section. As before, we will focus on the generic behaviour of our system. In 
particular, we will ignore all choices of the hybridisations (and call them non-generic) 
where the conductance and transmission phase as a function of V^, A, U and s deviate 
from the behaviour observed for general Tf if the dimension of the manifold of these 
hybridisations is smaller than six. By computing approximately 6500 data sets in the 
{U,A,s,r) - parameter space (~ 5 • 10^ data points G{Vg)) we have verified that the 
behaviour presented in this section is indeed the generic one. 

Interacting Case, Nearly Degenerate Levels 

For nearly degenerate levels, nonzero interactions of order F and arbitrary {F, s} we 
find three resonances separated by U in the conductance if the gate voltage is varied. 
They are of almost equal height and width, the latter always being of the of the order 
of F rather than of F; or U. The transmission phase a changes by vr over each of the 
peaks and jumps by vr at the transmission zeros located at ±U/2 in between. 

The average level occupancies show a complicated behaviour. As for the double dot 
geometry, they depend nonmonotonically on the gate voltage and show population 
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FIG. 4.9: Conductance G (black), transmission phase a (red), and average level occupancies 
(dot A: green, dot B: blue, dot C: orange) as a function of the gate voltage for parallel triple 
dots with nearly degenerate levels and different hybridisations F = {0.06 0.14 0.07 0.03 0.3 0.4} 
(upper left and lower right), F = {0.06 0.14 0.3 0.4 0.07 0.03} (upper right), F = 
{0.2 0.3 0.25 0.15 0.07 0.03} (lower left). 



inversion, but here it is not always the most strongly coupled dot that is depleted each 
time Vg is on-resonance. All this is shown in Fig. 14.91 

Another similarity with the case of two dots is the gradual appearance of additional 
structures in the conductance if the interaction is increased (Fig. 14.10]) . First, above a 
certain critical strength Uc depending on the hybridisations Tf and the relative signs 
s, the peaks located at Vg ^ ±U split up into two. For U > Uc the outer of the new 
peaks are still located at Vg w ±U, while the inner ones become more and more sharply 
centred at the transmission zeros at Vg ^ such that one would expect only the 

outer peaks to remain observable in the limit of large U . Unfortunately, we are not 
able to tackle interactions large enough to confirm this within our fRG approximation 
scheme (see next chapter). The additional resonances vanish if one increases the level 
spacing above a scale Acir <C F. This whole scenario is identical to the appearance of 
the correlation induced resonances for the double dot geometry. 

Furthermore, the third peak centred at shows very sharp structures that 

become further pronounced if the interaction is increased. Unfortunately, the precise 
shape of these 'sharp structures' is not universal, but most frequently they take the 
form of an additional peak (sometimes accompanied by dips on its left and right side) 
of width much smaller than F (Figs. 14.91 fc [4?T0]) . In general, the structures vanish for 
^CIR ^ r. If (and only if) we choose si = S2 = — and F^ < F^ -|- Fc, the dips in the 
conductance next to the additional resonance may even be transmission zeros for some 
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FIG. 4.10: Appearance of additional correlation induced structures in the conductance of par- 
allel triple dots with T = {0.03 0.07 0.04 0.06 0.2 0.6}. 



special choices of the hybridisations (Fig. 14.111 upper left panel). In particular, this is 
always the case for left-right symmetric Tf . Turning on slight asymmetries (which is the 
generic case) the zeros gradually disappear (Fig. 14.111 other panels). Surprisingly, they 
do not vanish if the level spacing is increased but always appear above a scale Adpl 
which depends on the dot parameters, in particular on left-right asymmetry (Fig. 14.12]) . 
For symmetric hybridisations, we have Adpl = 0, and this scale becomes larger and 
in particular no longer much smaller than F the more asymmetric the F^'^ are chosen. 
If it happens that Adpl and the level spacing where the crossover regime (see below) 
sets in become comparable, the additional zeros will never appear in between the limits 
A <C F and A > F. 

As for the double dot geometry, the correlation induced structures are accompanied by 
a steep change of the transmission phase. If for si = S2 = — the additional transmission 
zeros are present, a jumps by vr at each of them, so that we observe double phase lapses 
(DPLs) next to the sharp resonance. 

Due to the very sharp structures close to = 0, one word about the reliability 
of the fRG results is in order. At first sight they might be interpreted as an artifact 
of our approximation, especially since they become more pronounced the larger the 
interaction is chosen. However, there are strong indications that these structures are 
indeed believable features of our interacting quantum dot model. First, they appear 
consistently within the different fRG truncation schemes, and second, they are confirmed 
in some special cases that can be tackled by an NRG calculation believed to give very 
precise results. All this will be explained in much more detail in the next chapter. Here, 
we should only keep in mind that everything presented should be taken to be reliable, 
as promised in the introduction. 



From Small to Large Level Spacings: the Crossover Regime 

As mentioned above, the qualitative behaviour of the conductance as a function of Vg is 
independent of the choice of the hybridisations (as long as they are generic) for nearly 
degenerate levels. In particular, similar to the double dot case the only influence of the 
left-right or the A-B-C asymmetry is on the overall height of the conductance curves, 
on the scales Acir, Acir, and Adpl related to the correlation induced structures, and 
on the actual shape of the sharp feature close to = (Figs. 149) &: 14.1 1|) . This does no 
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FIG. 4.11: Illustration of the dependence of Adpl on left-right asymmetry Tf^/Tf. For level 
spacings within Adpl < A < Across, correlation induced double phase lapses (DPLs) are 

observed in the transmission phase a (red) of triple dots with s = { }. The couplings read F = 

{0.15 0.15 0.2 0.2 0.15 0.15}, F = {0.1 0.2 0.1 0.3 0.12 0.18}, F = {0.05 0.25 0.05 0.35 0.05 0.25}, 
F = {0.1 0.2 0.05 0.35 0.2 0.1} (from upper left to lower right). The conductance G (black) is 
shown as well. 
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FIG. 4.12: The same as Fig.lHIl but for F = {0.16 0.24 0.13 0.07 0.24 0.16} and different level 
spacings. The DPLs appear for Adpl < A < Across- 
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FIG. 4.13: Conductance G (black) and transmission phase a (red) as a function of the gate volt- 
age for parallel triple dots with generic level-lead hybridisations F = {0.1 0.2 0.06 0.14 0.4 0.1}. 



longer hold (see below) when the level spacing exceeds Across (r,s) and the dot merges 
into the crossover regime (as we will call the evolution from A <C F to A » F), which 
should be clear because similar to the double dot we would expect a{Vg) to be strongly 
dependend on s for large level detunings. 

In order to search for the different parameter regions that lead to distinct behaviour 
in the crossover regime, first, one has to convince oneself that there are only three 

independent choices of the two relative signs si and S2, namely {S1S2} = {++, H — , }. 

This can be seen by considering the particle-hole transformation di — > h\ under which 
a Hamiltonian of the type 



H=Y. Mil + "Y. {4i< - 5) (4*' - 5) 



{N -l)U 



d\di + U dldid\,di' + const. 



i>i' 



goes over to 



{N-l)U 



J (1 - hlhi) + U ^(1 - h}hi){l - hlhv) + const. 
^ i>i' 



(-^ - 2^^^ ) h\hi + uY,h\hih\,hv + const., 
I ^ ^ i>i' 



with being the number of levels in the system, 
that the choice si = — , S2 = + is equivalent to si = 



Hence, this transformation implies 
-|-, S2 = — if we additionally replace 



Vg — > —Vg- As argued above, the crossover regime will strongly depend on which of the 
three independent possibilities is realised. For s = {++}, it turns out to be independent 

of the size of the hybridisations, while for s = {-\ — , } it is the strength of the dot A 

which couplings differ in sign from the other two, sgn(f^t^) 7^ sgn(t^t^) = sgn(t^t^), 
that determines the behaviour of the conductance in between A <C F and A F. In 
particular, there are two distinct parameter regions in each case, namely F^ < F^ ~^'^c 
(regime W), and T ^>V j^ + T^ (regime S). 



The most simple case is that with s = {++}. It is depicted in Fig. 14.131 Starting 
from nearly degenerate levels, the correlations induced resonances close to = ±?7/2 
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FIG. 4.14: Gate voltage dependence of the conductance G (black) and transmis- 
sion phase a (red) for parallel triple dots with generic level-lead hybridisations F = 
{0.06 0.14 0.07 0.03 0.3 0.4}, illustrating the crossover from A <C F to A > F in regime 
S. 
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FIG. 4.15: The same as Fig. 14.141 but for generic hybridisations from regime W, F = 
{0.06 0.14 0.25 0.35 0.13 0.07}. 
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and the sharp feature at Vg ^ vanish (if they are present at all) for A Acir 
and A « Acir, respectively. The separation of the three remaining resonances (the 
Coulomb blockade peaks) increases, their lineshape gradually changes and for A » F 
we finally observe three Lorentzian resonances separated by C/ + A. They correspond 
to transport through the individual levels, and accordingly the height and width of the 
peak over which dot / gets depleted is determined by the hybridisations Tf" and Vf. As 
for A ^ r, the transmission phase changes by vr with the noninteracting single-level 
lineshape over each resonance and jumps by vr at the two transmission zeros in between. 

For s = {H — }, the transition from small to large level spacings is more complicated. 
As before, the additional features observed in the three-peak structure gradually dis- 
appear if the level spacing is increased. As mentioned, the subsequent evolution in the 
crossover regime depends on the relative coupling strength of the different dots. In- 
creasing A, one of the Coulomb peaks splits up (the left one 'corresponding' to dot A in 
regime S; the central resonance in regime W) while another becomes increasingly small 
(regime S: right resonance; regime W: left resonance). In analogy to the parallel double 
dot, we define a scale Across(r,s) characterising the set-in of the crossover regime (for 
example by the height of the resonance that almost disappears). This scale depends on 
the dot parameters, but it is roughly given by Across ~ F/10 > Acir, ■^cir- If for some 
parameters it happens that either Acir ~ Across or Acir ~ Across, the vanishing of 
the correlation induced structures and the split-off of one of the peaks in the crossover 
regime might not be distinguishable. Finally, we obtain three Lorentzian resonances 
identical to the s = {++} case for A ^> F, only that now the phase that changes with 
the noninteracting lineshape over each resonance shows a jump by vr between those 
corresponding to level A and B while evolving continuously between B and C. In par- 
ticular, a only jumps in between the peaks related to those dots which couplings do 
not differ in sign (sgn(f^t^) = sgn(t^t|j)). The transition from small to large level 
detunings for s = {-\ — } is depicted in Fig. 14.141 (regime S) and Fig. 14.151 (regime W). 

Finally, the behaviour of the conductance for s = { } in between the limits of small 

and large A can be described as follows. In regime S (shown in Fig. I4.16|) . the central 
peak (the one 'corresponding' to the dot A) splits up into three at a scale Across while 
the outer ones become vanishingly small. In contrast, the central resonance gradually 
disappears in regime W (depicted in Fig. I4.17j) . and the outer one that one would 
relate to the more strongly coupled dot splits up into two. For extreme choices of the 

hybridisations, namely F^ <C F and F^ Fc, the resonance close to = in regime 
W always remains a maximum surrounded by transmission zeros and corresponding 
phase jumps though becoming unidentifiable on a scale G of the order of the unitary 
limit. This similarly occurs in the noninteracting case and we will delay a more detailed 
characterisation of the affected parameter region. In general, the central resonance 
completely vanishes in the sense that it becomes a minimum and the transmission 
zeros disappear. As mentioned above, an additional correlation induced sharp peak and 
especially the double phase jumps (lapses) on both of its sides appear if A exceeds a scale 
Adpl depending on the hybridisations. Normally, this happens before the crossover 
regime is reached, but if in particular for large left-right asymmetries it happens that 
Adpl ~ Across, these additional features will never be observed. Finally, in both regime 
S and W we end up with the well-known Lorentzian three-peak structure, only that now 
the phase does not exhibit any jumps between any of them. This is consistent with the 
cases considered before, since here the couplings of the levels corresponding to both 
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The same as Fig. I4l4t but for s = { }, and F = {0.06 0.14 0.3 0.4 0.07 0.03}, 

: the generic crossover from A <C F to A ^ F for regime S. 



u/r=2, A/r=o.02, s={-} 



u/r=2, A/r=o.i, s={-} 





1 




0.8 




0.6 


CD 




CD 


0.4 




0.2 









1 




0.8 


CM 


0.6 


(D 




o 


0.4 




0.2 




n 




-3 3 

u/r=2, A/r=o.7, s={-} 




2 


1 


1.6 


0.8 


1.2 


0.6 


0.8 


0.4 


0.4 


0.2 










1 


2 


1.6 


0.8 


1.2 


0.6 


0.8 


0.4 


0.4 


0.2 










-3 3 

u/r=2, A/r=5, s={-} 




FIG. 4.17 
{0.06 0.14 



The same as Fig. 14.161 but with generic hybridisations for regime W, F 
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Region of A 


Typical values 


Related phenomenon 


A ^ 
A ^ i 


A/i = U.Ui 


Three 'Coulomb' resonances of equal height, width of 
order T and separation U for arbitrary {s,r}; addi- 
tional correlation induced structures if the interaction 
is large enough. 


AciR 


A/r = 0.05 


Vanishing of the CIRs located at « ±J7/2. 


AciR 


A/r = 0.05 


Vanishing of the sharp features located close to half 
filling. 


Adpl 


Depends on 

left-right 
asymmetry. 


Appearance of the additional transmission zeros and 
corresponding double phase jumps (lapses) for s = 

{ — } and Vb <rA + rc. 


Across 


A/r = 0.1 


Set-in of the crossover regime where some peaks be- 
come increasingly small while others split up. 


A>r 


A/r = 5.0 


Three Lorentzian resonances separated by J7 + A 
which height and width are determined by the cor- 
responding hybridisations Y\ and rp. The transmis- 
sion phase jumps by tt between those peaks related to 
levels which couplings do not differ in sign. 



TABLE 4.1: The different A - regimes of the triple dot geometry. 



pairs of subsequent resonances differ in sign, sgn(i^t^) / sgn(t^t^), and sgn(i^t^) 7^ 
sgn(i^tg). 

Effective Level Interpretation 

As in the double dot case, further insight in the underlying physics might be gained by 
considering the effective energies (the eigenvalues of the isolated dot Hamiltonian with 
renormalized parameters) and the new hybridisations 7^ (obtained from the correspond- 
ing eigenvectors) of our system in the limits A <C F and A » P. For small level spacings, 
large C/, and arbitrary s we find (again in analogy to [Silvestrov&Imry 2000]) that at 
each of the three Coulomb peaks the most strongly coupled level crosses the chemical 
potential of the leads // = upwards leading to a resonance (Fig. 14.151 left panel). This 
is possible since the effective levels do not depend monotonically on the gate voltage 
(which is again related to population inversion). In particular, near the transmission 
zeros located in between the peaks, the level with the largest 7^ crosses the chemical 
potential downwards, but roughly at the same point it crosses another level so that a 
zero instead of a peak is observed. This can be interpreted as a Fano-antiresonance, 
especially since it shows that the interaction enhances the bare asymmetry of the P; 
such that 7i of the most strongly coupled effective level is significantly larger than the 
largest P;. 

For large level spacings every depends linearly on the gate voltage and is only shifted 
by \J when one of the other levels is filled (Fig. 14.181 right panel). The hybridisations 
7i almost coincide with the original ones since virtually no hopping is generated by the 
flow. However, well-separated levels linearly crossing the chemical potential just lead 
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FIG. 4.18: Conductance G (black), effective eigenenergies uoi (solid lines) and hybridisations 
7,; = 7f + 7f (dashed lines) for parallel triple dots with T = {0.06 0.14 0.3 0.4 0.07 0.03} for 
two different level spacing. At A ^ F, the most strongly coupled effective level (red) crosses 
the chemical potential upwards at each resonance, while close to the transmission zeros this 
level crosses /i downwards but at the same time it crosses another level, giving rise to a Fano- 
antiresonance. The interaction enhances the bare asymmetry of the hybridisations. Note that 
due to the degeneracy of the levels for \Vg\ ^ oo, there is no need to recover the original F; in 
this Hmit. For A » F, all tj^ cross ^ linearly, hence s governs the behaviour of a as for U 



to the noninteracting Lorentzian resonances which width and height are determined by 
the corresponding hybridisations Vf and T^, only that now their separation is ?7 + A 
rather than A. Furthermore it follows that the relation between the relative signs of 
the level-lead couplings s and the behaviour of the transmission phase in between the 
peaks has to be the same as for [7 = 0. 



The Noninteracting Case 

Finally, we will consider the parallel triple dot geometry at [/ = 0. Similar to the double 
dot, the behaviour of the conductance in the region of large level spacings and in the 
crossover regime does not change if the interaction is gradually switched off. However, 
one observes the width of the V^-region where the conductance is significantly larger 
than zero to decrease when the level spacing is lowered. Hence in the limit A <C F we 
would expect fairly sharp structures on a scale of F arising from simultaneous transport 
through several levels in analogy to the sharp dip close to = in GiVg) of the double 
dot with ?7 = 0. Therefore it is reasonable to seek for analytical statements that follow 
from the exact solution of the noninteracting problem. Unfortunately, evaluating the ex- 
pression for the conductance (I4.14p with Q replaced by the free propagator yields a very 
lengthy expression that is hardly suited to reveal interesting information. Nevertheless, 
in the special case of left-right symmetry some insight in the number and position of the 
transmission zeros can be gained. For the double dot with finite level spacing, sticking 
to F;^ = F^ does not lead to non-generic behaviour and thus it seems reasonable that 
the same holds for the triple dot geometry. Regardless of that, we have numerically 
confirmed that everything stated here indeed represents the generic behaviour. Solving 
for the zeros of (I4.14p then yields 



^.«r.,,.^, (4.15, 
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with 

( 



/i(r,s) 



-4r^-2rij + o.5 



4rB-i 



i6rA+i6rB-2 

8rA+4rB-i 
i6rB-2 



/2(r,s) 



2 for s = ++ 

ler^ + ler^ - 2 for s = +- 

leTij - 2 for s = — 



and, most important, 



/3(r,5) 



64r^ - ler^ + 64rArB + sr^ + i6r| + 1 for s = ++ 

144r|-48rB + 132rArB + l fors = +- 

i44r| - 48rs + 64rArij + uy\ - ler^ + 1 for s = — 



We have impHcitly assumed that V is chosen as the unit of energy, implying that we 
can substitute Vc = 1 — — T^. The first imphcation of (|4.15|1 is that the number of 
transmission zeros is independent of the level spacing, and it turns out (by multidimen- 
sional analysis or simply by performing a 3d plot) that for s = {++, H — } the function 
/3({s,r} is larger than zero for [ryi,rs] = [0 . . . 1, . . . T^] so that their number is 
indeed two. This is consistent with the general (numerical) observation that the region 
of large A as well as the crossover regime are identical for both U = and ?7 > for 
arbitrary choice of {T} if one keeps in mind that for s = {H — } in both regime S and 
W the transmission zeros present for small level spacings lie outside the Lorentzian res- 
onances but do not vanish in the limit A ^> F. However, the second important insight 
that can be gained from the analytical solution is that the gate voltages with GiVg) = 
depend linearly on the the level spacing with a prefactor that is roughly of the order of 
r. Hence for A ^ F though still having the correct number of transmission zeros we 
expect to find a very sharp three-peak structure on the usual scale of the hybridisations 
(Fig. 14.191 upper left panel). The structures become broader when increasing A, but 
for s = {-I — } the crossover regime sets in at scale A < F, implying that the heights of 
the resonances begin to differ (Fig. 14.191 upper right panel). 

However, for s = {-\ — } (and of course for s = {++}) the behaviour of the conductance 
is identical within both the interacting and the U = case if one is only concerned about 
the number of transmission zeros and the position of the resonances. This turns out 

to be different if we choose s = { }, because here the function /3(s,F) is no longer 

positive for all hybridisations of interest. By solving /3(s,F) = we can compute the 
region of F where the conductance does not exhibit any zeros (Fig. 14.191 lower left panel). 
We find that there are two solutions of GiVg) = if dot B is strongly coupled, implying 
that the behaviour of the conductance in regime S is identical to the interacting case for 
all A. On the other hand, in most parts of regime W we observe no transmission zeros 
which is still consistent with the claim that the crossover regime and region of large 
level spacings are identical for both U ^ and U = 0. However, since their number 
is independent of A, additional zeros do not appear in the limit A <C F in contrast to 
the interacting case. Hence for small level spacings, we observe a very sharp three-peak 
structure but the conductance stays finite in between them, and, more important, the 
phase evolves continuously (Fig. 14.191 lower right panel). One exception to this is found 
for Ta,Tb Fc or Fc,F^ ^ F^ (the upper and lower left corner in the lower left 
panel of Fig. 14.19]) . Here, we find indeed two transmission zeros within the three-peak 
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FIG. 4.19: Upper row: Gate voltage dependence of the conductance G and transmission phase 
a (red) of noninteracting parallel triple dots with s = {+-}, T = {0.06 0.14 0.07 0.03 0.3 0.4}, 
A/r — 0.05 (left panel), and A/F — 0.2 (right panel). Lower row, left panel: Appearance of 
transmission zeros as a function of the hybridisations for the left-right symmetric, noninteracting 

case with s ~ { } for arbitrary level spacing. Lower row, right panel: The same as in the 

upper row, but for s = { }, F = {0.15 0.25 0.1 0.2 0.12 0.18}, and A/F = 0.05. 



structure for A ^ F. As the hybridisations are from regime W, one of these peaks splits 
up while the central one becomes small when the level spacing is increased. In contrast 
to the other parameters in regime W, the resonance close to = always remains a 
maximum surrounded by zeros and phase jumps, but it becomes too small to identify 
on scale of the order of the unitary limit. As mentioned above, identical behaviour is 
observed in presence of the interactions. 



4.1.4 Parallel Geometries with More Than Three Dots 
Four Levels 

In this section, we will generalise the previously considered models by adding further 
levels. We will focus on the A regions of most interest, which are in particular the 
limits A ^ F and A ^ F, and we refrain from providing a complete description of the 
crossover regime. 

First, we will study a four-level dot where the projected free propagator reads 

/1.5A \ 



0.5A 



-0.5A 



H,s + Hpp, (4.16) 



1.5A/ 
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and the projected part i^eff — Hpp has the well-known form (|3.11|) . 



with the definition 



s=L,R LI' 





Si 


S2 


S3\ 


Si 


1 


S1S2 


S1S3 


S3 


SlS2 


1 


S2S3 




S1S3 


S2S3 


1 / 



The conductance is computed via (I3.33p . which here reads 

.2 , 2 



G 



h 



(4.17) 



where we have defined 



SLl' ■= \ 



1,1' 



51 for I = B and arbitrary I' 

52 for Z = C and arbitrary I' 

53 ioi I = D and arbitrary I' 
1 otherwise. 



It is straight-forward to convince oneself that the six independent choices for the 
relative signs of the level-lead couplings s = {S1S2S3} are given by s = {-|- -|- -|-, -|- -|- 

— ,-| , -I h, — I — , }. As before, we introduce U as the strength of the 

two-particle interaction assumed to be equal between the electrons of all dots. The fRG 
flow equations comprise four for the on-site energies, six for the inter-dot hoppings and 
21 for the independent components of the effective interaction. 

For nearly degenerate levels, we observe four conductance resonances separated by 
U of almost identical height and width (Fig. I4.20p . The transmission phase changes 
approximately by vr when crossing each of them and jumps by tt at the three trans- 
mission zeros located in between. For increasing interaction strength, the peaks exhibit 
additional correlation induced structures at which the phase changes steeply. In par- 
ticular, we find features similar to the CIRs of the parallel two-level dot (Fig. 14.201 
upper right panel) as well as double phase lapses (jumps), the latter depending on the 
left-right symmetry of the hybridisations (compare the lower left panel of Fig. [^?20] with 
the lower left panel of Fig. I4.22p . 

In the limit of large level spacings, we recover the usual structure of six Lorentzian 
resonances of separation U + A corresponding to transport through the individual levels 
(Fig. I4.21|) . The conductance vanishes for one gate voltage in between those peaks 
associated with levels which relative couplings do not differ in sign while remaining finite 
otherwise. Correspondingly, the transmission phase exhibits jumps by vr in the former 
case and evolves continuously in the latter. If the level spacing is chosen somewhere 
in between the limits A ^ F and A » F, a crossover regime similar to the two- and 
three-level dot is observed, but we will refrain from giving a detailed description here. 

Just like before, the noninteracting lineshape G{Vg) differs strongly from that with 
[/ > in the limit of nearly degenerate levels. It shows very sharp structures on a scale 
F, and, more important, it does not exhibit four transmission resonances with a zero in 
between (Fig. 14.221 lower right panel) for arbitrary choice of the dot's parameters (the 
hybridisations and the relative signs of the couplings). 
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FIG. 4.20: Gate voltage dependence of the conductance G (black) and transmission phase a 
(red) for a four-level dot with nearly degenerate levels and generic level-lead hybridisations F = 
{0.05 0.15 0.075 0.225 0.025 0.075 0.1 0.3} for various parameters. For s = {+ - +}, these hy- 
bridisations are not generic, and hence we have chosen F = {0.12 0.28 0.1 0.3 0.04 0.06 0.03 0.07}. 
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FIG. 4.21: The same as Fig. 14.201 but for large level spacings. 
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FIG. 4.22: The same as the the lower left panel of Fig. 14.201 but for various combinations of 
parameters, illustrating the appearance of DPLs (left panel) and the lack of three transmission 
zeros in the noninteracting case (right panel). 
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FIG. 4.23: Conductance G (black) and transmission phase a (red) for a six-level dot with 
r = {0.02 0.08 0.05 0.15 0.03 0.12 0.07 0.23 0.02 0.08 0.04 0.11} and s = {+ - - + -}. For 
A <g; r, one would expect the height and width of the resonances to become more comparable 
and the change of a in between the transmission zeros to become precisely tt when A is slightly 
decreased down to A/F = 0.01. Due to the Hmitations of the fRG, we cannot confirm this 
scenario explicitly (as we encounter serious problems if both [/ » F and A <C F and if many 
levels are considered; this will be explained in detail in the last chapter). 



Six Levels 

We end this section with considering six parallel dots. The projected free propagator 
and the expression for the conductance can be written down in complete analogy with 
the four-level case, and also our results in the regimes A <C T and A » T turn out to be 
the same. For nearly degenerate levels, we observe six resonances and five transmission 
zeros in between them in the GiVg) curve (Fig. 14.231 left panel). The transmission 
phase jumps by vr at each zero and approximately changes by tt over each resonance. 
Further correlation induced structures are observed if the strength of the interaction is 
increased. If on the other hand the level spacing is the dominant energy scale, we find six 
Lorentzian resonances of separation U + /S. and transmission zeros (and corresponding 
phase jumps) between those associated with transport through levels which couplings 
have the same relative sign (Fig. 14.231 right panel). 
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FIG. 4.24: Gate voltage dependence of the conductance G (black) and transmission phase a 
(red) for parallel double dots with U/T = 6.0, degenerate levels, s = — , and generic level- 
lead hybridisations T — {0.27 0.33 0.16 0.24} for finite temperatures. Increasing T, the CIRs 
gradually disappear. 



4.1.5 Finite Temperatures 

Up to now we have only considered the T = limit. This might be a meaningful 
simplification since transport through quantum dots is generally measured at very low 
temperatures. Nevertheless, it is important to investigate whether the qualitative results 
derived in this section are stable against very small but finite temperatures, especially 
when one is concerned with the very sharp structures (for example the CIRs or the 
DPLs) one would expect to be smeared out at finite T. 

For the geometry which we are going to discuss in detail in the next section, a single 
level containing interacting spin up and down electrons, we have already argued that 
our fRG truncation scheme which neglects all frequency dependencies yields a spectral 
function of Lorentzian lineshape very different from the exact one. Since for arbitrary 
T the conductance is computed as an integral of the latter with the derivative of the 
Fermi function, we cannot expect to obtain reliable results beyond the limit of very 
small T (and we have indeed checked that this is the case), and the same should hold 
for other spinful geometries. Nevertheless, we can hope that things are different for 
spin-polarised models. Unfortunately, no other calculations for such models at finite 
temperatures have been performed up to now, rendering it impossible to check the 
trustworthiness of our results. Therefore we will only employ the T > calculations 
to check the aforementioned stability of the sharp correlation induced features against 
small temperatures T <T but refrain from a detailed study of arbitrary T. 
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FIG. 4.25: Critical temperature where the correlation induced resonances observed for par- 
allel double dots with degenerate levels vanish as a function of the two-particle interac- 
tion U. Left panel: Hybridisations F = {0.27 0.33 0.16 0.24}, and s = — . Right, panel: 
F = {0.5 0.25 0.07 0.18}, s = +. 



The conductance for finite temperatures is computed using the general expression 
(|3.28p , and we define the T / transmission phase as 



a^T) := arg [ / de g'{e) ^ tf-{tf,rGt^,{e + iO) 



(4.18) 



with g{e) being the Fermi function, and the expression for spin-polarised models follows 
by dropping all spin dependencies. First, we study the T - dependence of the CIRs 
observed for a parallel double dot with degenerate levels for sufficiently large interaction. 
It turns out that they are stable against temperatures T < F. To be more precise, the 
evolution of GiVg) with increasing temperature is similar to the one arising when U 
is lowered in that sense that the four-peak structure observed for T = merges into 
two single resonances if the temperature exceeds a critical Tent depending on all other 
dot parameters (Fig. I4.24|) . This temperature scale is roughly given by Tent ~ r/4, 
and it depends linearly on the interaction strength (Fig. 14.25]) . On the other hand, the 
transmission phase which shows a jump by vr at VJj = for zero temperature evolves 
continuously for T 7^ with a slope that decreases if T is increased. An interesting 
question is whether this 'smoothing' of both G{Vg) and aiVg) is caused only by smearing 
with the Fermi function. In fact, we have tested that the T - dependence of G'(V^) is 
qualitatively the same if one considers the finite temperature propagator at z = iO, or 
if one computes the latter at T = and calculates the conductance and transmission 
phase as an integral of this Q{e + id) with the Fermi function. 

Surprisingly, the sharp correlation induces features that appear for geometries with 
more than two dots (those similar to the ones observed for the triple dot at the central 
resonance) become more pronounced for T / 0. In particular, the double phase lapses do 
not vanish unless T is so high that the whole structure of GiVg) is smeared. The generic 
evolution of the conductance and transmission phase curves with increased temperature 
is shown in Fig. 14.261 The DPLs which show up because of the left-right symmetric 
hybridisations are more emphasised in the sense that they become less sharp. 

To understand the temperature dependence of G{Vg) for large level spacings, it is 
meaningful to consider a single spin-polarised impurity first. The conductance can then 
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FIG. 4.26: Conductance G (black) and transmission phase a (red) as a function of the gate 
voltage for a four-level dot with f//r = 1.0, A/T = 0.05, T = {0.1 0.1 0.15 0.15 0.05 0.05 0.2 0.2}, 
and s = { } for four finite temperatures. 



be computed exactly as an integral of the derivative of the Fermi with the (exact) 
spectral function. As one would intuitively expect, the transmission resonance centred 
at Vg = fj, becomes wider and its height decreases at temperatures T ^ T (Fig. 14.271 
left panel). The transmission phase evolves less steeply. Having this in mind, the 
temperature dependence of the peaks observed for interacting multi-level dots with large 
A is straight-forward to understand. The width (height) of each resonance increases 
(decreases) at a temperature scale given by the strength of the individual hybridisations, 
T ^ Ti. At the same time, the slope of the phase a becomes smaller and the jumps 
by TT at the transmission zeros are gradually smeared out, in particular if they are 
located close to a resonance. It is important to note that here it is no longer T/T but 
rather T/Ti which determines a typical temperature scale, implying that temperatures 
which are small for nearly degenerate levels have a considerable effect at large A if 
some levels are weakly coupled or if the number of dots becomes large such that 
T/T < T/Ftyp = T/{T/N). This is shown in the right panel of Fig. OTl 



4.1.6 Summary & Discussion 

In this section, we have discussed the generic transport properties of A^ quantum dots 
with level spacing A connected in parallel to source and drain leads. The latter were 
modelled by a noninteracting tight-binding approach while we allowed nearest-neighbour 
interactions U between the electrons on the different dots. Using the fRG approach, 
we calculated the full propagator of the system which determines the linear-response 
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FIG. 4.27: Left panel: Conductance G (black) and transmission phase a (red) for a single 
impurity with on-site energy Vg for T/T ~ 0.0 (solid lines), T/T = 0.5 (dashed lines), and 
T/r = 2.0 (dotted lines). Note that this is an exact result. Right panel: The same from fRG 
for a four level dot with U/T = 5.0, A/T = 5.0, T = {0.05 0.15 0.075 0.225 0.025 0.075 0.1 0.3}, 

and s = {H h} at zero temperature (black and red curves) and at T/T — 0.2 (orange and 

blue curves). 



conductance and the transmission phase. To fully describe transport, we computed G 
and a as a function of a variable gate voltage that simultaneously shifts the energies of 
all dots. 

For nearly degenerate levels A ^ F and nonzero interactions we found N transmission 
resonances separated by U of almost equal width and height, the former being of order 
F, with — 1 transmission zeros in between. The phase changes approximately by vr 
when crossing each peak and jumps by vr at the zeros. If the strength of the two-particle 
interactions exceeds a critical value depending on the dot parameters, additional cor- 
relation induced features in the A^-peak structure are observed. They vanish for level 
detunings larger than A {{s, F}) <C F, only the double phase lapses observed for geome- 
tries with three and more dots become more pronounced if A is increased. However, the 
scale Adpl that determines their appearance strongly depends on the left-right asym- 
metry of the hybridisations F^, and in particular it is no longer much smaller than F for 
generic asymmetric F^. Furthermore, the additional structures become exponentially 
sharp if the interaction strength is increased, and only the N usual 'Coulomb blockade 
peaks' remain visible. 

If the level spacing is increased above Across ({s, F}) < F, the lineshape of the A" 
peaks changes as the quantum dot enters the crossover regime. Some of the resonances 
split up while others become vanishingly small. In the limit of A F, we end up with 
A^ Lorentzian resonances separated by A + [/ corresponding to transport through the 
individual levels. The width and height of each peak is determined by the associated 
hybridisations F^ and F^^. The transmission phase a drops by vr over each resonance. 
Transmission zeros and corresponding phase jumps are located between those peaks 
related to levels whose relative level-lead couplings do not differ in sign, while the 
conductance stays finite and a evolves continuously otherwise. 

The regime of large level spacings is identical in both the U = and U > case, only 
that in the latter the separation of the resonances is increased due to Coulomb repulsion. 
For nearly degenerate levels, things are different. For some parameters, the conductance 
and transmission phase curves are similar in both cases if one is only concerned with 
the number of transmission zeros, but the observed structures are extremely sharp on 
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FIG. 4.28: Measurement of the linear-response conductance (blue lines; the y axis is scaled in 
arbitrary units) and transmission phase (green Hnes) of a quantum dot experimentally realised 
as in Fig. 13.11 as a function of a plunger gate voltage (x axis, scaled in mV) that lowers the 
energy of the dot region. The occupation TV of the dot was also measured and is show within the 
plots. (Reprinted by permission from Macmillan Publishers Ltd: Nature 436, 529, copyright 
2005.) 

a scale of T and the width of the individual peaks differ by orders of magnitudes of T. 
For other parameters, there are even less than — 1 transmission zeros in the {7 = 
case. Thus the two-particle interactions are a vital ingredient in order to find almost 
equal resonances with A^ — 1 phase jumps in between. 

As stated in the introduction, these results might serve as an explanation of the puz- 
zling behaviour of the transmission phase observed in a series of transport experiments 
through single quantum dots [AHZMU 2005, SBHMUS 1997, YHMS 1995], especially 
since in the very recent one ( [AHZMU 2005]) in addition to the transmission probabil- 
ity and phase also the average occupation of the dot was measured. In the conductance, 
resonances were found each time another electron was added to the system. It showed 
that for dots that were only occupied by a few electrons (~5; mesoscopic regime) the 
phase either evolved continuously or jumped by vr in between,^ depending on the peak 
considered and on the dot measured (Fig. 14.281 left panel), which is a typical 'meso- 
scopic' behaviour. Surprisingly, a universally jumped by vr between all resonances for 
every experimental realised dot under consideration if the occupation of the latter was 
already 'fairly large' (~15; universal regime. Fig. 14.281 right panel). This is the be- 
haviour which was observed in all previous experiments where the average occupation 
was not measured. As was already pointed out by [AHZMU 2005], the fundamental 
difference between dots populated by a different number of electrons is the level spacing 
between the states that are occupied if further electrons are added to the system. For 
dots from the mesoscopic regime (those with a few electrons) one would expect the lat- 
ter to be larger than for dots from the universal regime (which are occupied by a large 
number of electrons). Hence the size of the detuning between the levels should be taken 
as an input from the experiment if one seeks to set up a model to theoretically describe 
transport through such a single quantum dot. In this sense, our findings are consistent 
with the experimental observations. For A <C F (universal regime), our calculations 
yield N transmission resonances with A^ — 1 zeros and corresponding phase jumps in 
between, with A^ being the number of levels in the dot under consideration. This is pre- 
cisely the behaviour found in the experiment for dots from the universal regime. On the 
other hand, for large A ^> F the evolution of the phase in between the A^ conductance 



■^Note that in the experimental papers the gate voltage is defined with opposite sign. 
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peaks depends on the actual parameters of the dot, in particular on the relative signs 
of the couplings of the levels associated with consecutive resonances. These signs are 
connected to the parity of the underlying wavefunctions and cannot be controlled in an 
experimentally realised system. Hence the theoretical prediction would just be that for 
dots occupied by a few electrons the phase should jump in between some peaks while 
evolving continuously between other and that the observed behaviour should change 
with the system under consideration. Again, this is precisely what was observed in the 
experiment. 

A few things concerning our model calculations need to be pointed out. First, it is 
important to note that the behaviour in the universal regime follows from the presence 
of the two-particle interaction; it is not observed if one solves the noninteracting model. 
Second, the generic behaviour in both regimes is completely independent of all other dot 
parameters in the limits A <C T and A » F. The crossover regime and the scale Across 
where the latter sets in depend on the level-lead hybridisations and the relative signs 
of the couplings, but in the limits of small and large level spacings one always ends up 
with the same behaviour characterising the universal and the mesoscopic regime. The 
typical ratio used in the calculations isA^r/A^Fw 200, which corresponds to 
15 levels if one assumes a quadratic dispersion of the dot level energies, a value which 
is in agreement with the experiment. Of course, one might argue that in the latter 
one probed the phase evolution for arbitrary many electrons in between both limits 
occupying the quantum dot, and hence the level spacing should decrease continuously. 
Our calculations, however, yield resonances accompanied with non-universal phase be- 
haviour even for A from the crossover regime, consistent with what is measured for 
mesoscopic dots. Third, one has to contemplate about why the additional correlation 
induced structures were not observed in the experiment. This might be because the 
interaction is too small, the temperature too high, or because the spacing of the levels 
that are filled in the universal regime is still too large. However, this argumentation 
is flawed by the fact that it would lead to strongly non-universal behaviour depend- 
ing on the dot under consideration. In particular, there are no universal values of U, 
T or A such that the additional features do not appear (or are smeared) but at the 
same time the desired A^-peak structure with N — 1 fairly sharp phase jumps is still 
observed no matter how the couplings are chosen. For each individual dot, there are 
certainly choices of {U, T, A} such that no correlation induced phenomena but otherwise 
well-pronounced resonances and phase jumps are visible; however, another dot (with 
different couplings) might for the same choice {U, T, A} be in the crossover regime, or 
the temperature might be so high that the resonances or the jumps are already smeared 
in either the universal or the mesoscopic regime (remember that for the latter it is T/Ti 
which determines a typical temperature scale). A much more likely explanation why the 
additional features are not observed that would not spoil the universality for A <C F and 
A » F could be that the experimental resolution is much too small to observe these in 
general sharp structures which become even exponentially sharper with the interaction 
strength increased beyond Uc- Taking a look at the curves measured in the experiment 
(Fig. I4.28p . this does not seem to be too far-fetched. Fourth, one should note that it 
would of course be interesting to extend our model to include the spin degree of freedom 
(though the latter seems to play no role in the experiment (following [Silvestrov&Imry 
2000, Golosov&Gefen 2006]); note that, however, our findings for spinful models in the 
next section do not contradict the interpretation presented here), and to perform calcu- 



4.1. SPIN-POLARISED DOTS 



81 



lations at finite temperatures (though T is the smallest energy scale in the experiment; 
see [AHZMU 2005]). Finally, we want to point out very clearly that we do not seek 
to fit the experimental curves by our calculations. It is of course important that the 
transmission resonances we computed are of comparable width and height and that the 
phase almost changes by vr over each of them. Nevertheless, the actual lineshape of the 
curves measured in the experiment might deviate from our findings, maybe because we 
use an oversimplified model, or maybe because the resolution in the experiment does 
not allow for a unique fit of the measured points. The most important observation, 
however, is that in the two-level case at large interactions C/ » P, our calculations yield 
curves which quantitatively resemble those observed in the experiment, in particular if 
one is concerned with the separation of the resonances and the S-like structure of the 
phase. Employing NRG calculations [KHODM 2006] one can show that the same holds 
for multi-level dots where we are unable to tackle large interactions U ^ T. However, all 
this is of minor importance since it does not influence or main argument: for small level 
spacings (universal regime), independent of the dot's parameters we find transmission 
resonances with A^ — 1 phase jumps in between, while for large level spacings (meso- 
scopic regime) a evolves continuously or jumps by vr depending on the system under 
consideration. Thus we might hope that our calculations provide a natural explanation 
of the experiments. 
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4.2 Spinful Dots 

In this section we will study transport through dots including the spin degree of freedom, 
where the local interaction of spin up and down electrons leads to Kondo physics [Hew- 
son 1993] at sufficiently low temperatures. Signatures of the Kondo effect were first 
observed in 1934 as an increase of the resistivity of gold when the temperature was suf- 
ficiently lowered. The set-in of the regime where the resistivity does no longer decrease 
with decreasing temperature was found to depend on the quality of the sample under 
consideration, and it was only observed at all if the latter contained local magnetic 
moments. In 1964 it was finally Y. Kondo who succeeded in fitting the experiments by 
quantum mechanical calculations based on a Hamiltonian that describes free electrons 
antiferromagnetically coupled to a spin one-half magnetic impurity. Below a charac- 
teristically energy scale, the Kondo temperature, this magnetic moment gets screened 
by the delocalised free electrons which leads to an increased scattering and hence the 
resistivity increases as well. If, on the other hand, a single quantum dot has an average 
occupation close to an odd integer number, it shows a resonance in the conductance due 
to quantum fiuctuations between two spin-degenerate states. In fact, it can be shown 
that a quantum dot with a single electron in the state of highest energy modeled by an 
Anderson impurity type Hamiltonian is similar to a magnetic spin one-half impurity. In 
contrast to metals, an increased scattering leads to an increased conductance through 
the quantum dot (because additional transport channels are available). 

First, we will study the Kondo resonance in transport through a one-level quantum dot 
and especially show how the Kondo temperature, the typical energy scale that destroys 
this resonance, can be extracted even from the linear-response conductance at T = by 
the application of a magnetic field. We will then report on transport properties of up 
to four of these dots coupled in series, which can be viewed as a short Hubbard chain. 
Such systems have been studied recently using the numerical renormalization group 
(NRG) [ONE 2005,Oguri&;Hewson 2005, ZBRR 2006, Nisikawa&Oguri 2006]. Here, we 
will demonstrate the power of the fRG to scan the parameter space for regions usually 
not covered by NRG. In particular, we will show that choosing non-symmetric couplings 
to the leads (which in general would correspond to the generic case) results in a vanishing 
of transmission resonances. Next, we will describe the so called two-stage Kondo effect 
that is characteristic for transport through a quantum dot to which another dot is 
attached by a tunnelling barrier, the so-called side-coupled geometry. This situation 
has been of recent interest as well [Cornaglia&Grempel 2005, Zitko&Bonca 2006], and 







FIG. 4.29: Dot geometries considered in this section. 
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we will again demonstrate the power of the fRG to scan the whole parameter space, 
this time verifying that the results obtained by the aforementioned authors using the 
NRG indeed represent the generic situation. Finally we will again turn to the system of 
parallel double and triple dots studied in the previous section, but here we will describe 
the physics that arises if we include the spin degree of freedom. These geometries have 
not been investigated in the literature so far. 



4.2.1 A Single Impurity 

Application of the fRG - the Kondo Resonance 

The free propagator with the leads projected out (13. lip reads 

iLO — Vg + aB/2 + I sgn(u;) F 

with B being a magnetic field applied to the system, Vg the gate voltage that shifts 
the on-site energy of the dot, F the sum of left and right hybridisations, F = F^ + F/j, 
and we have (as usual) performed the wide band limit. The fRG flow equations for the 
self-energy and the effective interaction (|2.40|) and (I2.4ip take the form 



C/2(A) 

_2 U'^{A)V^{A)Vi{A) 

[(A + F)2 + F^2(A)][(A + F)2 + y/(A)] ' 

where G is obtained from (I4.19P by the replacement Vg — aB/2 V^, and the initial 
conditions read Va{A oo) = Vg - aB/2 and U{A oo) = U, with U := W'"' being 
the local interaction of spin up and down electrons. 

If we solve these differential equations numerically and compute the linear-response 
conductance, 

we recover precisely what was claimed in the introduction. The Lorentzian lineshape 
that corresponds to the noninteracting case is gradually transformed into a box of 
width U where the conductance reaches the unitary limit if the interaction is increased 
(this is shown in the left panel of Fig. 14.301 for the left-right symmetric case; choosing 
7^ changes only the height of the resonance). For gate voltages within this 
Kondo resonance the mean occupation of the dot is odd and the spin degeneracy of the 
states with one electron give rise to an increased scattering amplitude and can therefore 
be understood as the reason for the perfect conductance. Despite the fact that the 
transmission phase a and level occupation are directly related to G by a generalised 
Friedel sum rule [Hewson 1993] (which holds exactly in our approximation since we are 
deriving an effective noninteracting model), the latter is shown as well in Fig. 14.301 
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FIG. 4.30: Left panel: Gate voltage dependence of the conductance G (solid lines) and average 
level occupation (dashed lines) for a single-level dot with Tl/T = Tn/T = 0.5 and zero magnetic 
field for different local interactions. If the interaction is large enough (U/T ^ 1), the Lorentzian 
lineshape of G(V^) observed in the noninteraction case is altered and the Kondo effect leads to 
a wide plateau of almost perfect conductance in a region of width U centred around = 0. 
Within this plateau, the average occupation of the dot is constant (and in particular odd). Right 
panel: The same as the left panel, but for fixed U/T ~ 15.0 and various magnetic fields. The 
spin degeneracy is lifted and the Kondo resonance gets destroyed at B w Tk- This implies that 
the Kondo temperature can be extracted from the linear- response conductance at T = 0. For 
B/T ~ 0.1, the partial conductance of the spin up (dashed) and spin down (dotted) electrons 
is shown as well. 



Since we have dropped all frequency dependencies in the fRG flow equations and 
mapped the interacting system to an effective noninteracting one, for arbitrary Vg the 
spectral function will always be a Lorentzian centred around ^^(0). This lineshape 
differs strongly from the exact one (obtained from the analytic Bethe ansatz solution 
[Hewson 1993], or very precise NRG data [GDCO 2000]). In particular, neither the 
resonance at u; = is very sharp, nor the spectral function exhibits the signatures of 
atomic physics for large U, the Hubbard bands. Hence, we cannot expect to obtain 
reliable results for the conductance at high temperatures. For T = 0, however, all that 
enters in the computation of the latter is the spectral function at zero frequency, and it 
turns out that the fRG successfully captures its pinning at u; = for gate voltages close 
to = 0. In particular, the renormalized level energy Va{0) is approximately zero in a 
region of size U ^ T around Vg = 0. For the fRG truncation scheme that keeps only the 
flow of the self-energy an exponential pinning of the spectral weight at a; = can even 
be shown analytically [KEM 2006]. If on the other hand one keeps the whole frequency 
dependence of the two-particle vertex, it is possible to compute a spectral function that 
shows a the sharp Kondo resonance as well as the Hubbard bands, although with much 
higher numerical effort [HMPS 2004]. 

Adding a magnetic field B to the system lifts the spin degeneracy at the particle- 
hole symmetric point and thus destroys the Kondo resonance. The wide plateau in 
G{Vg) present for i? = splits up and for large fields (in particular when B is the 
dominant energy scale) we recover two Lorentz-like resonances separated approximately 
hy U + B which correspond to individual transport of the spin up and spin down 
electrons (Fig. 14.301 right panel). This is similar to the spinless two-level dot with large 
level spacing. 
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FIG. 4.31: Left panel: The Kondo energy scale Bk of a single-level dot as a function of U for 
different gate voltages within the Kondo resonance, Vg = 0.0 (black), Vg ~ 2.0 (red), and Vg = 
4.0 (blue). The data points obtained from fRG where fitted as Bk = acxp (— \bU/T — cT/U\). 
The values computed for the fitting parameters are shown in the inset. Their error is less than 
one percent, and the fitting curves cannot be distinguished from the original data on the scale 
of this plot. Right panel: Gate voltage dependence of the fitting paramater c. 



Determination of the Kondo Energy Scale 

A meaningful definition of the energy scale on which the Kondo resonance survives 
could be the magnetic field strength Bk required to reduce the conductance down to 
half the unitary limit (with Vg implicitly assumed to be chosen such that G{Vg) is of 
almost unitary height, i.e. within the conductance plateau). If on the other hand one 
tackles the system by an NRG approach that gives the correct spectral function, one 
usually defines this energy scale, which is then called Kondo temperature Tk, as the 
width of the Kondo resonance at zero frequency. In the exact solution, its dependence 
on the dot paramaters is given by 



TK = C{U,Vg)exp 



ttU 

8 r 



2 r2 [/ 



(4.21) 



where the prefactor C{U,Vg) is constant to leading order in Vg and U [Hewson 1993]. 
Hence it is meaningful to probe whether the fRG results for the Kondo energy scale 
defined by the magnetic field dependence of the conductance give reliable results by 
computing Bk for several values of U and Vg and fit the resulting curves Bk{U) by a 
function of the form 



f{U/V)=a{Vg)ei^V 



KVg)^-c{Vg)^ 



(4.22) 



This is done in Fig. 14.311 It turns out that the fitting parameter h 0.32 is almost 
independent of Vg and indeed close to the exact value 7r/8 « 0.39. Rirthermore, the 
factor c depends approximately quadratically on the gate voltage (Fig. 14.311 right panel), 
while the Vg dependence of prefactor a is only weak. Altogether, this shows that the 
exponentially small Kondo energy scale can indeed by extracted from the linear-response 
conductance computed by our fRG approximation scheme at zero temperature. 
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4.2.2 Linear Chains of Dots 
General Considerations 

In this section we study transport througii a system that contains N single-level dots 
connected in series, a short Hubbard chain. The projected free propagator for this 
geometry reads 



lUJ 



B 



/ i sgn(u;) Tl ti2 
t2i -Ai 



tlN 



\ 



-A 



N- 



_i + i sgn(cj) TrJ 

(4.23) 

a magnetic field B, the gate 
voltage Vg that simultaneously shifts the on-site energies of the different levels, and 



with real hoppings tij = tji, arbitrary level spacings Aj 



hybridisations with the left and right lead := Trltj^l piead and Tji := vr|t}^| piead, 
respectively. In the following we will mainly focus on the special case where only 
equal nearest-neighbour hoppings t := fj^j+i are present and all on-site energies are 
the same, implying Aj = 0. Furthermore, we will assume that there is only a local 
interaction U := U^f between spin up and down electrons, which we will also choose 
to be independent of the level considered. Later on, we will relax all these assumptions 
and show that all our main results remain unchanged. It is important to note again 
that within our fRG approach this neither causes conceptional problems nor extreme 
numerical effort since all parameters enter the flow equations only via the free propagator 
(or the initial conditions), leaving at maximum ((2A^)^ -|- (2A^)^ — symmetries) ordinary 
differential equations to solve. We will not consider the effects of a magnetic field here. 

If we solve the fRG flow equations (12.401 I2.4ip . we can as usual calculate the full 
propagator of the system ^ = ^(A = 0) and compute the linear-response conductance 
at zero temperature 



G 



(4.24) 



Even in the case that we start out only with nearest-neighbour hoppings and local 
interactions, it turns out that in general all elements of the self-energy and two-particle 
vertex (except those that violate spin-conservation and anti-symmetry) such as hoppings 
between arbitrary sites and long-range interactions are generated by the flow. 



A Chain of Three Dots 

We will start with three dots with nearest-neighbour hoppings t and equal hybridi- 
sations with the left and right lead, respectively. It is again instructive to consider 
the noninteracting case first, which is very simple if we choose a basis where the dot 
Hamiltonian is diagonal. Since the eigenvalues of the latter are given by u;2 = Vg and 
u;i^3 = Vg ^ \f2t and the corresponding normalized eigenvectors read (~"^)0) "^) and 
(^,±-7j,i), transport through the system can be understood by three parallel levels 

with level spacing ^/2t, hybridisations T^/T = T^/T = 0.25 and rf g/T = T^^g/T = 
0.125, and relative signs of the level-lead couplings s\ = — , and S2 = —■ Hence for 
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FIG. 4.32: Gate voltage dependence of the conductance G (black) and total average occupation 
(red) of a linear chain of three dots with Tl/T — Tn/T ~ 0.5 for two different nearest-neighbour 
hoppings t for the noninteracting case and for the largest U/T that can be tackled by our fRG 
approximation scheme. Well-pronounced plateaus of unitary conductance and especially the 
wide plateau at = caused by the suppression of charge fluctuations are only observed for 
larger interactions (see for example t/T « 8, U/T ^ IGtt in [ONH 2005]). 



large nearest-neighbour hoppings t compared to the strength of the level-lead hybridi- 
sation r = F/, -|- F/j (the molecular orbital regime) we expect to find three Lorentzian 
conductance resonances of width 2{Tf + F^) separated by V2t, and the total average 
occupancy n of the dot system to decrease by two (due to the two spin directions) over 
each peak. The transmission phase is no independent quantity since it is connected 
to n by tto- = 7rre/2 [Nisikawa&Oguri 2006], implying that it changes by vr over each 
peak and evolves continuously in between, as required by the two relative signs si^2 = — • 
If the nearest-neighbour hopping is decreased and becomes comparable to F, the three 
levels overlap and contribute simultaneously to the transport. The resonances merge 
and for t <C F only one peak at = is observed. 

If we turn on an local interaction U between the spin up and down electrons, the 
Kondo effect alters the GiVg) lineshape. For a chain comprising N dots, the sharp 
Lorentzian resonances away from the particle-hole symmetric point in the molecular 
orbital regime are transformed into box-like peaks of width U /N, while the resonance 
at = becomes a box even wider than the others. [Nisikawa&Oguri 2006] suggested 
to interpret this as a strong suppression of charge fluctuations at half-filling of the chain. 
The average number of electrons in the dots is consistent with the single-level Kondo 
behaviour. For gate voltages on-resonance, the occupation is odd and almost constant 
which leads to a half-integer total spin and hence a conductance plateau at the unitary 
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FIG. 4.33: Gate voltage dependence of the conductance G (solid lines) and total average occu- 
pation (dashed lines) of a linear chain of three dots with Tl/T = Tr/T = 0.5 for six different 
choices of the nearest-neighbour hoppings ty, level spacings A^, local interaction U\oc and 
inter-level interaction Un '■= U^'^ = U^'j' present between all electrons of different levels. 
Left panel: iia/F = ias/r = 3.6, AJT '= 0.0, i7ioc/r = I.Stt, Un/T = 0.0 (black curves), 
ti2/r = 2.7, t23/r = 3.5, Ai/r = -0.7, As/r = -0.3, C/ioc/r = 1.5^, Un/T = 0.0 (red 
curves), t^/T = tas/r = 3.0, A,/r = 0.0, C/ioc/r = Un/T = 1.5tt (green curves). Right 
panel: iia/T = taa/r = 8.0, A,/r = 0.0, C/ioc/r = 2.07r, Un/T = 0.0 (black curves), 
ti2/r = 7.4, i23/r = 9.5, Ai/r = 0.9, As/r ^ -0.9, C/ioc/r = 2.07r, Un/T = 0.0 (red 
curves), t^/T = tas/r = 8.0, A^/T = 0.0, C/ioc/r = Un/T = 2.0n (green curves). 



limit. Between the resonances the occupation is even and transport is suppressed. Since 
for a chain with an odd number of sites half-filling corresponds to an odd number of 
electrons, the suppression of charge fluctuations around Vg = leads to a wide region 
of almost perfect conductance. 

If the interaction becomes much larger than the separation of the peaks (which is 
still of order of the nearest-neighbour hopping t), or t is smaller than the hybridisations 
with the leads such that already the noninteracting chain is no longer in the molecular 
orbital regime, the resonances begin to merge, and only one wide plateau of perfect 
conductance around Vg = is observed (the local regime). Within our approximation 
scheme, however, this parameter region is out of reach. A much more detailed discussion 
of the limitations of the fRG truncation scheme that neglects the frequency dependence 
of the two-particle vertex as well as the flow of all higher-order vertex functions will 
be given in the next chapter. In particular, by explicit comparison to very accurate 
NRG data we will establish upper boundaries for the interaction strength that can still 
be tackled within this approach for every system under consideration. A comparison 
between the noninteracting case and the largest U that can be treated for a chain with 
three sites is shown in Fig. 14.321 

Relaxing the assumption of direct hoppings being present only between nearest- 
neighbours and that all levels have the same on-site energy does not lead to a qual- 
itative change of the results. Adding an inter-level interaction between all electrons 
of different levels shifts the resonances further outwards, as should be expected due to 
the additional Coulomb repulsion. Surprisingly, now the outer peaks become wider and 
their box-like structure evolves earlier than for the resonance at half-filling, at least for 
the interactions we can treat here. This is shown in Fig. 14.331 

Applying a left-right asymmetry A := Tji/Tl between the hybridisations with the 
left and right lead has a more dramatic effect if the hopping t is comparable to the 
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FIG. 4.34: Left panels: Gate voltage dependence of the conductance through a Unear chain of 
three dots for the noninteracting case and for U /T = 1.0 for four different left-right asymmetries 
A := Tr/Tl. For J7 > 0, the central resonance vanishes if A is increased. Right panels: 
Eigenenergies Ui (solid lines) and corresponding hybridisations with the leads 7^ and 7/^ of the 
isolated dot system with renormalized parameters for U/ttT = 1.0 as a function of Vg for two 
different asymmetries A ^ 1 (dashed lines: Tf + Tf-) and A = 9 (dotted lines: Tf, dashed 
lines: Ff ). For reference, the conductance is shown as well (black long-dashed Hne; the y-axes 
is scaled to one). 



hybridisation strength, t ^ T. While for the noninteracting case increasing A basically 
leads to an overall decrease of the conductance (which is more or less obvious for not too 
small t since the height of each peak arising from one level crossing the chemical potential 
decreases with increased A), the peak at = completely vanishes in presence of a local 
interaction if one moves away from the left-right symmetric case that was treated before 
(Fig. 14.341 left panel). This might be understood as follows. Due to spin degeneracy, we 
can describe our system by three effective energies u!i{Vg) (the eigenenergies of the dot 
system decoupled from the leads at the end of the flow) , which for the noninteracting 
case cross the chemical potential separately, leading to a transmission resonance no 
matter how A is chosen. One would not expect this picture to break down if a small 
interaction is turned on, and indeed it turns out that also in this case for j4 = 1 one 
effective level crosses the chemical potential approximately at the position of each peak 
(Fig. 14.3^ upper right panel). Surprisingly, the same holds for very large asymmetries 
(Fig. 14.341 lower right panel). In fact, the effective energies for this case nearly coincide 
with the ones for A = 1 (which implies that the total average occupation of the chain 
barely changes as well). What causes the resonance at = to vanish is that the new 
hybridisations 7/" (V^) and 7/^(V^), which almost coincide with the noninteracting ones 
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FIG. 4.35: Gate voltage dependence of the conductance G (black) and total average occupation 
(red) of a linear chain of four dots with Tl/T = Tn/T ^ 0.5 for two different nearest-neighbour 
hoppings t for the noninteracting case and for the largest U/T that can be tackled by our 
fRG approximation scheme. To observe well-pronounced plateaus of unitary conductance, the 
interaction needs to be even larger (see for example t/T « 8, U/T ^ IGtt in [Nisikawa&Oguri 
2006]), while the 'Mott-insulating' type of behaviour at half-filling is already visible. 



in the left-right symmetric case, are strongly altered by the interaction for A ^ 1. In 
particular, the asymmetry and coupling strength of the level that crosses the chemical 
potential at half-filling are both increased in those regions of Vg where it significantly 
contributes to the transport. In contrast, the left-right asymmetry of the levels that 
cross at the outer peaks is slightly reduced while their hybridisation strength remains 
almost unchanged compared to the A = 1 case, at least close to the point where they 
cross fj,. The combined effect of the increased asymmetry and level broadening explains 
the vanishing of the resonance at half-filling. For a nearest-neighbour hopping from the 
molecular orbital regime, this effect is no longer observed. 



A Chain of Four Dots 

The situation of a chain containing four dots is very similar. For [7 = and large 
hoppings t ^ T we observe four resonances that merge if the nearest-neighbour hopping 
t becomes comparable to the hybridisation with the leads T = Tl+Tr. In the molecular 
orbital regime (large t) a local interaction U transforms the Lorentzian lineshape of 
the peaks into boxes of width U/N due to the set-in of the Kondo effect. The total 
occupation, which in the noninteracting case changes by two at each resonance and 
remains constant in between, gradually becomes step-like, in particular odd for gate 
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voltages within the resonances and even in between. Again, charge fluctuations are 
strongly suppressed close to = 0, but this time half-filling corresponds to an even 
number of electrons in the chain so that a conductance valley is observed instead of a 
plateau of unitary height. Even if t is lowered such that in the noninteracting case the 
conductance is not small at = due to the overlap and decreased separation of the 
two nearest levels, the interaction leads to an exponential reduction of the conductance. 
Hence due to the suppression of charge fluctuations at half-filling the four-site chain 
tends towards a Mott insulator. 

One should note that the behaviour of GiVg) at = observed here holds for 
arbitrary long chains. If the latter comprises an odd number of sites, there is a wide 
conduction plateau of unitary height at the particle-hole symmetric point, while for 
an even number of dots half-filling corresponds to an even number of electrons and 
transport is strongly suppressed. 

The effects of the interaction in a four-site chain were computed by [Nisikawa&Oguri 
2006] employing an NRG approach. Our fRG approximation scheme allows to treat 
values of U as large such that the strong suppression of charge fiuctuations at Vg = 
shows up, but we are unable to treat U large enough to observe strongly-developed 
box-like resonances (see next chapter). The conductance G{Vg) for the noninteracting 
case and the largest interaction that can still be tackled by our approach is depicted 
in Fig. 14.351 Introducing arbitrary hoppings, level spacings, left-right asymmetric hy- 
bridisations, or an additional inter-level interaction between the dots does not alter the 
results qualitatively. 



4.2.3 The Side-Coupled Geometry 
General Considerations 

The side-coupled geometry consists of a single-level dot connected to the leads (the 
embedded dot) to which another dot with one level (the side-coupled dot) is attached 
by a tunnelling barrier t (see Fig. 14.36]) . The projected free propagator reads 

[g^A^^)] - = - + 4 + ^S^(-) f - + I) , (4.25) 

where i? is a magnetic field applied to the system, Vg denotes the gate voltage that 
controls the on-site energies of both dots, and T^^ji are the hybridisations with the 
left and right lead in the wide-band limit. As usual, we choose the total hybridisation 
strength F = F^, -|- as the unit of energy. The flow equations (|2.40p and (|2.41l) yield 
the full propagator ^ = ^(A = 0) of the system in presence of the interactions, which in 
general will be a local interaction U := U^j^^ between spin up and spin down electrons 
on both dots and a nearest-neighbour interaction Un '■= U^f = U^f between electrons 
on the two different dots. The zero temperature conductance can then be computed as 
(13:341) . 

Here Pa{i'^) denotes the spectral function of the embedded dot. 

When considering the parallel double dot in the next section, it will prove helpful 
to relate this problem back to the side-coupled one. In fact, it is possible to map the 



92 CHAPTER 4. NUMERICAL RESULTS 



FIG. 4.36: The side-coupled geometry. 

side-coupled dot onto a parallel double dot by a basis transformation to bonding and 
antibonding states, 



{|e,cr), \s,a)} 



\A,a) 



V2 



(|e,cr) + \s,a)), \B,a) 



V2 



(|e,cj) - \s,a)) 



with \e{s),a) being the single-particle wavefunction of the embedded (side-coupled) dot. 
The noninteracting part of the Hamiltonian Hip describing the side-coupled geometry, 



H 



IP 



E 



l=e,s 

crCO,a,L + tRdl^„Co^a,R + H.C 



then maps onto 



H 



IP 



E 



{Vg + aB/2 + t)nA,a + (Vg + aB/2 - t)nB,a - 



s=L,R 



^4<xC0,a,. + H.C. 



This is just the noninteracting part of a Hamiltonian describing a parallel system of 
single-level dots with level spacing 2t and level-lead couplings ij^R = if, r = The 
mapping of the interaction part of the side-coupled Hamiltonian, 



is very simple for equal local and nearest-neighbour interactions U and Un- Namely, 

= ne,^{n'A,l + nB,i + n^.t) + ?^s,T('^A,i + nB,i) + ne,ins,i 
= (nA,T + nB,^){nA,i + nB,i) + ng^-fn^,! + ne,jns,j 
= (nA,t + nB,^)inA,i + nB,i) + nA,inB,] + nA,inB,i 
= H;^t/U, 

where we have used that for both spin directions the following identity holds, 

4cg Cg Cg 

(ca + Cij)(cA + cb)(4 - 4)(CA - cb) 
-ic\ + 4)(c|4 - 4)(CA + cb)(ca - cb) 

('^a'^b ~ c^c|4)(cBCA - cacb) 
4nAnB. 
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Alltogether this shows that an equal integration between all electrons of the side-coupled 
dot maps to the same kind of interaction of the parallel two-level dot. Unfortunately, 
this becomes more complicated if initially only a local interaction U in the side-coupled 
geometry is present. In this case, also correlated hoppings are generated by the trans- 
formation, 



U 
2" 



Small Inter-Dot Hoppings: Two- Stage Kondo Regime 

The mapping between both geometries will be helpful to gain further insight into the 
interacting double dot studied in the next section. In the noninteracting case, however, 
the transport through the parallel two-level dot with spin is by definition equivalent to 
transport through the spin-polarised one which has been studied in Sec. 14.1.21 This im- 
mediately gives insight into the behaviour of the side-coupled geometry if no interaction 
is present. In particular, if the inter-dot hopping t is small compared to the hybridisa- 
tion strength F, the conductance as a function of the gate voltage Vg has a Lorentzian 
lineshape with a dip close to the particle-hole symmetric point. Since for t = the 
side-coupled geometry maps onto a non-generic situation of the double dot which in 
particular does not exhibit a dip in G{Vg) at = 0, the separation of the 'resonances' 
for small but finite t is very poor (for example compared with the situation shown in 
Fig. 14.21 for degenerate levels). This dramatically changes if we add a local interaction 
U of the spin up and spin down electrons on the embedded and side-coupled dot to the 
system. At zero temperature, in absence of the side-coupled dot the spin on the em- 
bedded dot gets screened by the lead electrons which gives rise to the Kondo resonance 
in the conductance (the first stage Kondo effect). Hence, coming from large negative 
gate voltages, G increases with a slope determined by U/T similar to the single-level 
case (Fig. 14.381 left panel). If the side-coupled dot is present its spin gets screened as 
well (the second stage Kondo effect), so that at some point the conductance falls off 
and merges into a valley where transport is strongly suppressed (Fig. 14.371 left panel). 
Again in analogy to the single-level case, it is roughly U/t'^ that determines how steeply 
G{Vg) falls off (as slightly indicated in Fig. I4.38|) . Thus, if the inter-dot hopping is 
small compared to the hybridisation strength, already for [/ ^ T we would expect a 
well-pronounced valley around Vg = 0, while the Kondo resonance on the embedded dot 
has not evolved yet and hence the conductance should vanish slowly when increasing 
\Vg\. If the interaction is large compared to the hopping t, the width of the region where 
transport is suppressed is given by U (Fig. 14.381 left panel) and is independent of the 
actual choice of t (Fig. 14.381 right panel), which is consistent with the single-level case. 

The behaviour of the average level occupancies is in agreement with this picture. For 
U ^ t^, the side-coupled dot is either occupied by two electrons or unoccupied outside a 
region of width U around Vg = 0, while at Vg = ±U/2 the occupation abruptly changes. 
The occupation of the embedded dot falls off much more continuously from Vg = —00 
to = -|-oo because U /F is small. All this is similar to the single-level case, but the 
behaviour within the Kondo valley is different. The occupancy of the side-coupled dot 
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FIG. 4.37: Gate voltage dependence of the conductance G (black), transmission phase a (red) 
and average level occupancies (green: embedded dot, blue: side-coupled dot) of the side-coupled 
geometry with Tl/T = Tn/T = 0.5 and zero magnetic field for parameters from the two-stage 
Kondo (left panel) and the local spin-singlet (right panel) regime. For the former, the two-stage 
Kondo effect causes a valley of width U where transport is strongly suppressed. In contrast, 
the local spin-singlet phase can be understood as individual transport through the bonding and 
anti-bonding states of the isolated molecule. 



does not show a plateau but depends non-monotonically on the gate voltage. This also 
holds for the embedded dot. In contrast, the total occupation depends monotonically 
on Vg which is necessary because of the behaviour of the conductance to which it can be 
directly related by a generalised Priedel sum rule [Cornaglia&Grempel 2005]. For the 
side-coupled geometry, this relation reads G/e^h~^ = 2 sin^(n7r/2), which implies that 
G almost reaches the unitary limit for an occupation close to an odd integer number 
while it vanishes if there is an even number of electrons in the dot region. This again 
justifies why transport is strongly suppressed if the Kondo effect is active on both the 
embedded and the side-coupled dot. The transmission phase a, which is also related to 
the conductance by a Friedel sum rule, changes by vr over each of the two 'resonances' 
and jumps by vr in between, which is similar to the noninteracting case. Both the 
average level occupancies and the phase for a small hopping from the two-stage Kondo 
regime are also shown in Fig. 14.371 

Choosing asymmetric couplings to the left and right lead only lowers the overall height 
of the conductance curve G(V^) but does not lead to new physics (Fig. 14.391 upper left 
panel), which is similar to the single-level case. Relaxing the assumption of equal local 
interactions on both dots only leads to quantitative changes. An important observation, 
however, is that the strength of the interaction on the side-coupled dot determines the 
width of the Kondo valley for small CZ/F, whereas its dependence on the the local 
interaction strength on the embedded dot is only weak (Fig. 14.391 upper right panel). 
This manifests the interpretation that at T = the Kondo effect on the side-coupled dot 
is active and causes the wide valley of small conductance. Adding a nearest-neighbour 
interaction to the system spreads up the whole curve but does not lead to new physics. 

The energy scale that determines the breakdown of Kondo physics for this geometry 
can be determined from the (numerically) exact T = spectral function obtained by 
NRG, which shows a sharp resonance at zero frequency characteristic for the (first stage) 
Kondo effect of the embedded dot and an even sharper dip at w = within this resonance 
due to the (second stage) Kondo effect on the side-coupled dot. The width of these two 
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FIG. 4.38: Gate voltage dependence of the conductance G of the side-coupled geometry with 
left-right symmetric hybridisations and zero magnetic field for various parameters from the 
two-stage Kondo regime. 



resonances then define the energy scales and at which the corresponding Kondo 
resonance gets destroyed. As indicated, it turns out that » T^. This is consistent 
with the single-level dot where the Kondo temperature is proportional to exp(— [//F), 
and hence for small inter-level hoppings we would expect the energy scale determined 
by U /T to much larger than T|. determined hyU/t^. Thus, if we move away from T = 0, 
the conductance at = should increase because the spin on the side-coupled dot is 
no longer screened above a temperature T|-. If the T exceeds the second characteristic 
energy scale, Tj^^ also the electron on the embedded dot is no longer screened and thus 
one would expect the conductance to decrease again. All this was shown in an NRG 
calculation by [Cornaglia&Grempel 2005]. As explained above, since the fRG truncation 
scheme we use neglects all frequency dependencies we cannot expect to obtain reliable 
results for finite temperatures. In contrast, it is no problem to tackle another quantity 
that destroys the Kondo resonances, a magnetic field, to extract the two Kondo energy 
scales. In analogy to the single-level dot, we define the Kondo 'temperature' T^- of the 
second stage of the Kondo effect as the field required to raise the conductance at = 
to half the unitary limit, while the Kondo scale is just that of the single-level dot 
obtained from the side-coupled geometry if we set i = 0. The fRG results obtained from 
the magnetic field definition of and TJ- are consistent with those obtained from the 
spectral function within the NRG approach. A more detailed comparison will be given 
in the next chapter. 

The magnetic field dependence of the conductance within the two-stage Kondo regime 
is shown in the left panel of Fig. 14.401 As explained above, for S > the conductance 
increases and the valley around Vg = Q evolves into a maximum since the second stage 
of the Kondo effect is no longer active above an energy scale T"^. For small magnetic 
fields B 0...10T|^, the conductance curves are qualitatively identical with those 
obtained by [Cornaglia&Grempel 2005] for = but finite T. This shows again that a 
magnetic field destroys the Kondo resonance similar to the temperature, justifying our 
definition of the Kondo energy scales. Further increase of the magnetic field strength 
lets the conductance around = decrease again since also the Kondo resonance on 
the embedded dot gets destroyed for B w T^. In the limit of large fields we finally 
recover twice the spin-polarised structure in the G{Vg) curve, which corresponds to the 
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FIG. 4.39: Left panel: Conductance of the side-coupled geometry as a function of Vg for B = 
and different left-right asymmetric hybridisations with the leads A := Tr/Tl for one parame- 
ter set from the two-stage Kondo regime (upper row) and one from the local spin-singlet phase 
(lower row). Right panel: The same for symmetric hybridisations, but different local interac- 
tions on the embedded and side-coupled dot Ue and Us and an additional nearest-neighbour 
interaction Un (denoted as U = {Ue, Us, Un})- 



noninteracting one if only a local interaction is present. 



Large Inter-Dot Hoppings: Local Spin-Singlet Phase 

The physics in the limit of large hopping between the embedded and the side-coupled 
dot is very different from the two-stage Kondo regime. In the noninteracting case, which 
can again easily be understood by the mapping on the double dot geometry, and at zero 
magnetic field the conductance shows two Lorentzian peaks of width T separated by 2t. 
These resonances correspond to transport through the bonding and antibonding level of 
the isolated two-dot 'molecule'. The average occupancy drops by two at the two peaks. 
The transmission phase changes by vr over each of them and jumps by vr in between 
since side-coupled dot maps onto a double dot with relative sign of the hybridisations 
s = +. Turning on a local interaction on both dots gradually transforms the Lorentzian 
resonances into box-like structures and enlarges their separation (Fig. 14.371 right panel) 
due to the set-in of the Kondo effect on the bonding and antibonding level. For large U, 
the average occupation within each resonance is constant and close to an odd integer. 

For large interactions, it is again meaningful to treat the isolated two-dot 'molecule' 
decoupled from the leads to understand the width and the separation of the resonances. 
The single-particle part of the Hamiltonian of the single-level double dot with param- 
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FIG. 4.40: Conductance G as a function of Vg of the side-coupled geometry with Tl/T ^ 
Ffl/r = 0.5 for various zero magnetic fields. Left panel: Two-stage Kondo regime with T^/T = 
0.0096 and T^/F = L16 for B/T|. = 0.0 (black), B/Tj. = 1.0 (red), B/T^ = 10.0 (green), 
B/T^ = 100 w 1.2T^/r|. (blue), and B/T^. = 500 « &T]^/Tj. (orange). Right panel: Local 
spin-singlet phase with T^/F = 0.11 for B/T]^ = 0.0 (black), B/T}^ = 0.5 (red), B/T^ 1.0 
(green), B/Tj^ = 5.0 (blue), and B/T}^ = 50.0 (orange). 



eters arising from the mapping of the side-coupled geometry is diagonal with entries 
Vg — U /2 + t and Vg — U/2 — t. Hence the vacuum and the single-particle state of lowest 
energy are degenerate \iVg = t + U/2. For higher gate voltages, the isolated molecule 
is unoccupied. In the occupation number basis, 

{\nA,^,nA,i,nB,^,nB,i)} = 

{|1,1,0,0),|0,0,1,1),|0,1,0,1),|1,0,0,1),|0,1,1,0),|1,0,1,0)}, 

the two-particle part of the Hamiltonian reads 



2Va-U + 



(2t + U/2 

+U/2 



V 



+U/2 
-2t + U/2 



+U/2 -U/2 
-U/2 +U/2 



0/ 



The smallest eigenvalue, which is obtained from diagonalising the first block of this 
matrix, is given by 2Vg — 0.5U — O.bVU"^ + 16t^. This eigenvalues and the lowest of 
the single-particle part of the Hamiltonian are degenerate if Vg = - t + OWU^ + IQt^, 
hence for gate voltages in the interval {t + U /2, —t + 0.5\/C/^ -|- 16t^) the single-particle 
state is that of lowest energy and thus at T = the molecule is occupied by one 
electron. Alltogether this means that we would expect a transmission resonance of 
width 2t + U/2 - 0.5Vt/2 + 16*2 centred at Q.2b{^^TPTlGfi + U) (and likewise for 
negative gate voltages due to particle-hole symmetry), which is indeed observed for 
large interactions [Zitko&Bonca 2006]. As for the chain of dots studied in the previous 
section, we are unable to reach regions of U where the conductance shows completely 
box-like structures which width and position is given by these values within our fRG 
approach. 
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The eigenvector that corresponds to the two-particle state of lowest energy at half- 
filing is of the form 



which is a spin-singlet state. Therefore the 'molecular' regime for large inter-dot hop- 
pings is also called local spin-singlet phase. 

An asymmetry in the hybridisations with the left and right lead only results in an over- 
all decrease of the conductance (Fig. 14.391 lower left panel). Relaxing the assumption 
of equal local interactions or introducing an additional nearest-neighbour interaction to 
the system leads only to qualitative changes but no new physics (Fig. 14.391 lower right 
panel). 

Applying a magnetic field B to the system lifts the spin degeneracy and hence the 
Kondo resonances of the bonding and antibonding levels get destroyed. Similar to the 
single-level case the box-like resonances split up, whereas the conductance at half-filling 
always remains a minimum. As before, the Kondo energy scale Tk can be defined 
by the magnetic field required to lower G to half the unitary limit within the peaks, 
and again the results obtained from this definition yield values consistent with those 
calculated in an NRG approach [Cornaglia&Grempel 2005] (see the next chapter for 
more details). Furthermore, the GiVg) curves for small magnetic fields B G [0,Tk] are 
similar to those computed by the aforementioned authors for finite temperatures, which 
again illustrates that a small magnetic field and small nonzero temperatures indeed 
similarly destroy the Kondo resonance, justifying our definition of the energy scale Tk- 
For large magnetic fields we recover twice the spin-polarised structure in G{Vg). Hence, 
if only local correlations are present we recover four Lorentzian peaks. The influence 
of a magnetic field in the molecular-orbital regime is again shown in the right panel of 



As the evolution of the conductance curves for small magnetic fields at T = and finite 
temperatures in both regimes is similar, the former provides a zero temperature criterion 
whether a certain choice of parameters corresponds to the two-stage Kondo regime or 
the local spin-singlet phase. The distinction between the two is not unambiguously 
possible if one only considers the gate voltage dependence of the conductance since the 
latter is strongly suppressed at half-filling in both cases. For all parameters shown in 
this section, we have tested the evolution with finite magnetic fields of G{Vg) to justify 
that they indeed correspond to the regime claimed. 

4.2.4 Parallel Double Dots 

In this section, we will generalise the results derived for the parallel double dot geometry 
in Sec. 14.1.21 to the case including the spin degree of freedom. The geometry is again 
depicted in Fig. 14.411 The free propagator, which is by definition diagonal in spin space, 
reads 



A|l, 1,0, 0) + /x|0, 0, 1, 1) = (A|^) + fx\B)) (in) - lit)) 

= (A|e)+/i|5))®(|Ti)-|iT)), 



Fig. 11201 



iOJ — Vg + 



A+aB 
2 



+ i sgn 



{u;)Ta i sgn(^) [yr^ + ^yrM 



isgn(^) yr^+.^rK iu;-Vg 



A-aB 
2 



V 



+ i sgn(cj) 



(4.27) 
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FIG. 4.41: The double dot geometry. 

where Vg is the gate voltage that shifts the on-site energies of the dots, A a level 
detuning, and is a magnetic field applied to the system. The hybridisations with the 
left and right leads are denoted by T^'^, the relative sign of the level-lead couplings by 
s, and as usual we define T/ :=rf + rf, and T := F^i + rB. If we introduce correlations 
to the systems, which will in general be a local interaction U[f between the spin up 
and spin down electrons on each dot as well as nearest-neighbour interactions U^j^ and 
U^f, we have to solve the flow equations (I2.40p and (|2.41l) to obtain (an approximation 

to) the full propagator ^ = ^(A = 0). The conductance for each spin direction can 
than be calculated in analogy to the spin-polarised case, (I4.10p . 



,e2 



= 4- 
n 

(4.2^ 

The parallel double dot geometry allows for a large amount of parameters to be 
varied. Fortunately, the noninteracting case is well understood since it is by definition 
equivalent to the spin-polarised case which was discussed in detail in Sec. 14.1.21 In 
presence of correlations, it seems reasonable to stick to two special choices out of all 
possible interactions, namely the one where only local correlations on each of the dots 
are present and that with equal interactions between all electrons. The first situation 
can be interpreted as two spatially well-separated dots with one level each, while the 
second would be a reasonable model for a single dot containing two levels. Despite the 
fact that it might be rational to introduce a direct hopping t between two separated dots, 
we will consider t = here since such a hopping could always be tuned away by a basis 
transformation of the dot states. As for spin-polarised dots, we will focus on studying 
the generic behaviour of all quantities of interest. In particular, we would expect that 
all the G{Vg) curves in dependence of the other parameters are qualitatively similar for 
all level-lead hybridisations up to a few pathological cases. As in the spin-polarised case 
we have checked that this is indeed true by computing roughly 5000 points within the 
parameter space {U, A, B, F^'^, s}. 

Purely Local Interactions, s = -|- 

We will start with s = small level detunings A <C F and the case where only local 
correlations U between the electrons on each of the two dots are present. Already 
for small interactions C/ w F, the 'Fano-antiresonance' observed in the noninteracting 
case is transformed into a wide valley of width U where the conductance is strongly 
suppressed. At Vg ^ ±U/2 two resonances show up whose wings become smaller for 
\Vg\ > U/2 if the interaction is increased (Fig. 14.421 lower left panel). An increased 
left-right asymmetry in the hybridisations leads to an overall decrease of G (Fig. 14.421 
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FIG. 4.42: Comparison of the double dot at s = + and purely local interactions U with the side- 
coupled geometry. Upper row: Gate voltage dependence of the conductance for the side-coupled 
dots with Tl = Tfl = 0.5, U = 1.5, t = 0.3 (black) and of the parallel double dot (red) with 
U = 1.5, A = 0.6 for T = {0.125 0.125 0.125 0.125} (left panel), and T = {0.05 0.1 0.25 0.1} 
(right panel) . For clarity, all energies are given in arbitrary units. Lower left panel: Conductance 
through the double dot geometry with degenerate levels, F = {0.1 0.2 0.5 0.2} for U/T = 1.0 
(black), U/T = 1.5 (red), and U/T = 2.0 (green). Lower right panel: The same but for 
U/T = 2.0 and F = {0.1 0.2 0.5 0.2} (black), F = {0.1 0.2 0.6 0.1} (red), F = {0.04 0.26 0.6 0.1} 
(green), and F {0.26 0.04 0.6 0.1} (blue). 



lower right panel). Not surprisingly, all this is in complete analogy to the side-coupled 
dots in the two-stage Kondo regime. As explained in the previous section, the latter can 
be mapped onto the parallel double dot geometry with ^-i?-symmetric hybridisations 
and additional interaction terms. However, by going over to arbitrary hybridisations 
and slowly changing the interaction to a purely local one we convinced ourselves that 
the physics does not change. The conductance as a function of the gate voltage has 
qualitatively the same structure (Fig. 14.421 upper row), and, more important, the valley 
at half filling disappears when turning on small magnetic fields (see below). Alltogether, 
this means that the physics of the double dot geometry with s = +, small level spacings 
and purely local interactions can be understood from the side-coupled dot in the two- 
stage Kondo regime. 

Since we know that for the side-coupled geometry with direct hopping t it is 
that determines whether the second stage of the Kondo effect is active or not, for 
small t/T the interaction already becomes important for fairly small U/T and hence we 
should expect the fRG only to give reliable results for [/ < F. At the beginning of this 
chapter, we have claimed that we were only going to show data for interaction strength 
where the fRG still produces assured results. In all situations under consideration 
up to now this was indeed fullfilled, as will be shown in the next chapter by explicit 
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FIG. 4.43: Gate voltage dependence of the conductance G (black), transmission phase a (red) 
and average level occupancies (dot A: green, dot B: blue) for the parallel double dot with generic 
level-lead hybridisations F = {0.1 0.2 0.5 0.2} and purely local interactions U for four different 
level spacings A. For small level detunings, the phase was hidden within the region where the 
conductance is strongly suppressed (see text). For A/F = 10.0, the conductance for U/T = 5.0 
is shown as well (black dotted line). 



comparison to very accurate NRG data. However, for the parallel double dot geometry 
an interesting parameter regime is certainly that with degenerate levels, where by the 
aforementioned argumentation the fRG results should be questionable already for fairly 
small interactions. This indeed turns out to be the case, but again we will delay an 
appropriate discussion to the next chapter, since it proves that the fRG still covers the 
essential physics. Here, we should only keep in mind, that for s = + and small level 
detunings close to (away from) half filling we will capture the behaviour of G{Vg) and 
all other quantities of interest qualitatively (quantitatively). 

Similar to the noninteracting case, the transmission phase a changes by vr over each 
of the resonances. Within the valley, it is almost constant up to several jumps by vr. 
They are caused by additional transmission zeros numerically appearing because the 
conductance is far underestimated. Hence these phase jumps are almost certainly an 
artifact caused by the application of the fRG at too large U. The same holds for the 
average level occupancies. For large negative gate voltages, both dots are occupied by 
two electrons. The more weakly a level is coupled, the more steeply its occupation falls 
off close to Vg = —U/2 and merges into a plateau within the valley, which also seems 
to be overestimated. All this is depicted in the upper left panel of Fig. 14.431 

Increasing the level spacing up to A w F does not change the lineshape of G{Vg) 
substantially (Fig. 14.43]) . which is plausible since even for t/T w 0.5 we found the side- 
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FIG. 4.44: Gate voltage dependence of the conductance G (black) for the double dot geometry 
with purely local interactions, left-right symmetric hybridisations F = {0.15 0.15 0.35 0.35}, 
and s = — . For these parameters, G{Vg) (and in particular its four-peak structure) can be 
understood from the behaviour of the conductance Gg (dashed lines) and transmission phase as 
(dotted lines) of two single-level dots with F^/F — Tr/T = 0.15 (red) and F^/F = Fj^/F = 0.35 
(green) separated by A (note that F still refers to the hybridisation strength of the double dot) . 



coupled geometry to be within the two-stage Kendo regime (Fig. 14.38]) . An interesting 
observation, however, is that for such level detunings (where the fRG is known to give 
reliable results) the average occupancies depend non-monotonically on the gate voltage 
within the conductance valley (Fig. 14.431 lower left panel), a scenario similar to the 
side-coupled geometry. For large A » F, we finally recover two separated resonances. 
Their Lorentzian lineshape exhibited at C/ = is transformed into box-like structures 
for large interactions U /Ta,b ^ 1 because of the Kondo effect that is active on one dot 
at each peak. This picture is supplemented by the fact that each level occupancy shows 
the typical Kondo behaviour at one of the resonances while being constant otherwise. 
The same holds for the transmission phase (up to a phase jump of vr in between the 
resonances). In addition, the height of the individual peaks is determined by the left- 
right asymmetry of the hybridisations F^/F^ and F^/F|j. For large U and A (and 
[/ < A) the position of the two resonances can be understood as usual by diagonalising 
the isolated dot system. The lowest energies of states with one, two, three and four 
electrons read Vg - A/2 - C//2, 2Vg - A U/2, 3Vg - A/2 + C//2, and AVg + 3[//2, 
and hence at zero temperature the dot is occupied by an odd number of electrons in 
the Vg-regions [A/2, A/2 + U] and [-A/2 - [/, -A/2]. Assuming that in analogy to 
the single-level case transport is possible in these situations, it follows that we would 
expect two peaks of width U separated by A. 

Applying a magnetic field to the system closes the analogy between the side-coupled 
geometry and the parallel double dot. For small A, the conductance close to half filling 
is first increased for finite fields B <^T due to the suppression of the second stage Kondo 
effect. At some i? < F, it starts to decrease again and for large i? » F we recover twice 
the spin-polarised structure in the lineshape of G(V^), which is here identical with the 
noninteracting one (Fig. 14.471 upper left panel). For large A, the situation for finite 
magnetic fields is again completely similar to the single-level case at each resonance. 
The box-like structures split up on an energy scale determined by the corresponding 
Kondo temperatures proportional to exp(— [//Fyi^s) (Fig. 14.471 lower left panel). For 
large A, we recover the noninteracting lineshape for the partial conductance and spin- 
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FIG. 4.45: Gate voltage dependence of the conductance G for the double dot geometry with 
degenerate levels, s = — , and generic hybridisations F = {0.1 0.2 0.5 0.2} for purely local 
interactions U (left panel) and in presence of additional nearest-neighbour correlations Un 
(right panel) . The evolution of four resonances in the left panel can be easily understood in the 
left-right symmetric case (Fig. I4.44| . 



dependent transmission phase of the spin up and down electrons (Fig. 14.471 lower 
right panel). 

Purely Local Interactions, s = — 

If the relative sign of the level-lead couplings is s = — , we can no longer map the 
parallel double dot geometry to the side-coupled one. Fortunately, for the special case 
of left-right symmetric hybridisations the flow equations for the self-energy and the 
effective interaction can be related to those of a single-level, allowing for an analogy 
between both situations within our approximation. In particular, if Vf = Vf, the free 
propagator (I4.27P is diagonal, and if only local correlations are present, no off-diagonal 
terms are generated by the flow. This means that the full propagator of the system is 
diagonal with entries that correspond to two individual single-dot problems each with 
interaction between spin up and down electrons U and level-lead hybridisations Ta 
and Tb- The conductance G of the system is then proportional to the absolute square 
of the difference of both single-level propagators evaluated at z = iO, implying that 
transport through the system can be interpreted as quantum interference between two 
independent channels each corresponding to a single spinful impurity. 

In Fig. 14.441 we show the conductance through the double dot system with purely local 
interactions and s = —. For degenerate levels and U exceeding a critical value depending 
on the hybridisations, we observe four transmission resonances instead of two. They 
gradually disappear when A is increased. As stated above, this can be understood by 
considering the G{Vg) curves of the individual levels. In general, the box-like structure 
of these curves is developed with different strength because of the different fZ/F/, and 
hence both curves will intersect somewhere around Vg « ±C/ /2, at least if the interaction 
is large enough for the width of both peaks to be determined by U rather than by F^^^. 
Since the same holds for the transmission phase, the conductance of the double dot 
system shows three transmission zeros at VJ^ = and Vg ±fZ/2, which implies a 
total of four resonances whose separation should increase with U (see also Fig. 14.451 left 
panel). A nonzero level spacing shifts two of the zeros closer together and they finally 
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FIG. 4.46: The same as Fig. 14.431 but for s = — and purely local interactions U jV = 4.0. In 
each plot, the conductance for U jV = 1.0 but otherwise identical parameters is shown as well 
(black dotted Hues). The crossover from small to large A is analogous to the spin-polarised 
case. 



disappear for A = 0{U). At first sight this might sound surprisingly since the G{Vg) 
curves of the individual levels shifted by A against each other still intersect twice for 
arbitrary large level level spacings. However, this does not hold for the transmission 
phase. In fact, the latter determines the conductance by a Friedel sum rule, G = 2sin^ a, 
but this does not imply that on the other hand the conductance fixes the phase uniquely. 
Hence, an intersection of the two GiVg) curves is not necessarily accompanied by equal 
transmission phases, so that it is no contradiction that for A = 0{U) only one of the 
zeros in the conductance remains. 

Relaxing the assumption of left-right symmetric hybridisations leads to a small but 
finite diagonal entry in the free propagator which stays small during the fRG flow. Hence 
one would not expect to find qualitatively different behaviour for such more general 
hybridisations, and this is indeed the case (Fig. 14.451 left panel). The resonances prove 
to be robust against fairly large additional nearest-neighbour interactions and they only 
disappear if all interactions in the system are approximately of equal strength (Fig. 14.451 
right panel). 

The transmission phase a jumps by vr at each of the three transmission zeros, but its 
change over each resonance is less than vr. In the left-right symmetric case, the average 
level occupancies are exactly those of two single-level dots with hybridisations Ta,b, and 
again for arbitrary level-lead couplings the results are qualitatively the same (Fig. 14.461 
upper left panel). 

As mentioned above, two of the three transmission zeros present if the interaction is 



4.2. SPINFUL DOTS 105 




FIG. 4.47: Conductance G as a function of the gate voltage for a parallel double dot with 
generic level-lead hybridisations T — {0.1 0.2 0.5 0.2} and purely local interactions U for various 
magnetic fields. Lower panels: Note that the Kondo temperatures for the two Kondo resonances 
determined from the corresponding single-level problems read T^/T = 0.01 and /T = 0.25. 
The right panel shows G{Vg) for B /T ~ 30.0 (black curve). The partial conductances (red: spin 
up electrons, green: spin down electrons) and the transmission phases (a^: blue, a^: violet) 
are presented as well. 



large enough vanish if the level spacing becomes of order U . The third zero is shifted to 
larger Vg and the corresponding resonance that is located at even larger gate voltages 
becomes vanishingly small. The peak that is moved outwards and gradually disappears 
is always the one one would associate with the more weakly coupled level (better: it is 
the right resonance that vanishes if Ta is more weakly coupled) . For further increase of 
the level spacing, the conductance in the region where originally the other three peaks 
were located evolves into one wide resonance which splits up and finally we recover two 
Kondo resonances for A » F. This is completely similar to the s = -|- case except for 
the fact that there is no transmission zero and corresponding jump by vr in between the 
peaks. If U is so small that for degenerate levels only two resonances are observed, one 
of them splits up while the other is shifted outwards and becomes small if A is increased. 
One should note that this whole scenario is completely identical to the crossover regime 
of the spin-polarised double dot studied in Sec. 14.1.21 (aside from the fact that it starts 
out with four instead of two resonances for U sufficiently large) . 

The magnetic field dependence of the GiVg) curves is similar to the s = + case. 
For degenerate levels, the conductance close to half-filling first increases and evolves 
into a maximum when a field is turned on (Fig. 14.471 upper left panel). For left-right 
symmetric hybridisations, this can again be understood from the interfering single-levels 
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FIG. 4.48: The appearance of the novel resonances for a parallel double dot with equal local 
and nearest-neighbour interactions J7, degenerate levels, and s = + in dependence of U for fixed 
level-lead hybridisations F = {0.1 0.2 0.5 0.2} (left panel) and in dependence of F for fixed U 
(right panel). 



dots. Since there the conductance at half filling is suppressed at an energy scale that 
exponentially depends on the strength of T^ and r^, one of the single-level Giyg) curves 
will fall off much earlier than the other one when a magnetic field is turned on. If i? 
exceeds a certain strength, the conductance at = falls off again and for large B we 
recover twice the spin-polarised, i.e. the noninteracting, structure. For equal strength 
of the interaction, the scales on which the conductance first increases and then begins 
to decrease again are of the same order of magnitude for both the s = -|- and s = — 
case (compare the upper panels in Fig. 14.47]) . This suggests a small energy scale, or, in 
other words, a large effective interaction within the system. This is surprising since one 
would expect a large effective interaction to be recognised by a breakdown of the fRG 
scheme for fairly small \J compared to the 'real' energy scales of the system. However, 
for s = — no such breakdown was observed (see the next chapter for the clarification of 
the word 'breakdown' in cases where no reference data known to be precise is available). 

Interactions Between All Electrons, s=+ 

Next, we study the case where we have equal local and nearest neighbour interactions \J 
within the system. This should be an appropriate model for a single dot containing two 
levels. The easiest case is that with relative sign of the level-lead couplings s = -|- since 
for A-5-symmetric hybridisations it can be mapped onto the side-coupled geometry with 
equal local and nearest neighbour interactions. For small level detunings A <C F close to 
half-filling we find the well-known valley of width \J of suppressed transport associated 
with the second stage Kondo effect. For degenerate levels, by increasing the interaction 
strength the wings of the conductance peaks become flat and at a certain critical \Jc 
depending on the hybridisations two additional peaks split off so that we find a total 
of four transmission resonances located at « ±?7/2 and Yg « ±?7. The additional 
peaks close to the valley of suppressed conductance become increasingly sharp when 
the interaction gets stronger. This whole scenario is quite similar to the appearance of 
additional resonances (the CIRs) within the spin-polarised model discussed in Sec. l4.1.2l 
The conductance valley at half-filling and the evolution of the four-peak structure when 
increasing the interaction strength is depicted in the left panel of Fig. 14.481 The critical 
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FIG. 4.49: Stability of the new resonances against small level detunings (left panel) and asym- 
metric interactions (right panel). The hybridisations read T = {0.1 0.2 0.5 0.2}. In the left 
panel, the features of G'(V^) in the valley are artifacts of the fRG approximation scheme (see next 



chapter). The parameters of the right panel are [/^''^/F = U'^^'^g/T ~ U'^j'/T ~ U^j^ /T ~ 3.0 



(back), C/;;^/F = 3.1, [/^-"s/F = 2.9, U;'f /T = 2.9, 
C/^'* /F = 2.6, C/;f /F = 2.7, U^f /T = 3.2 (green). 



(T,cr 
I.] 

U^f/T = 3.0 (red) and C/^^/F = 3.2, 



Uc only depends on the difference between the left-right asymmetries of the different 
dots and on the bare ^--B-asymmetry. The larger the former and the smaller the latter, 
the larger is Uc (|4.48l right panel). In general, it was always found to be of the order of 
magnitude of F for all systems considered. The additional resonances are robust against 
a very small level detuning A ^ F and against asymmetric interactions (Fig. 14.49]) . The 
former explains why they were not observed when considering the side-coupled geometry 
since there we only focused on t <T. 

For A ^ F, from large \Vg\ to the valley the transmission phase changes by vr over 
the resonances and jumps by vr close to = 0. If present, a evolves steeply when 
crossing one of the additional structures located at Vg ^ ±U/2 which is again similar 
to the CIRs observed for the spin-polarised case. The average level occupancies both 
fall off with a slope determined by the corresponding F/ when increasing Vg —U/2 
and merge into a plateau of width U with indications of a nonmonotonic dependence 
on the gate voltage. The details are shown in the upper left panel of Fig. 14.501 

Increasing the level spacing, only the outer resonances remain and their separation 
becomes larger, in complete analogy to the case where only local interactions are present. 
For A 3> F we recover the two Kondo resonances associated with the individual levels. 
Their position can be estimated by diagonalising the dot system decoupled from the 
leads. In this case, the lowest energies of states with one, two, three and four electrons 
read Vg - A/2 - 3U/2, 2Vg - A - 2U, 3Vg - A/2 - 3U/2, and 4Vg, and hence at zero 
temperature the V^-regions where the dot is occupied by one (three) electrons are given 
by [(A + U)/2, (A + 3f/)/2] ([-(A 3U)/2, -{A + U)/2]). Therefore we expect two 
conductance peaks of width U separated by A -|- 2U. The evolution from small to large 
A is shown in Fig. 14.501 

The magnetic field dependence of the curves G{Vg) is identical to the case with local 
correlations. The conductance at half filling first increases for fields much smaller than 
in the case considered before, which seems plausible because of the additional interac- 
tion terms. Further increase of B destroys the first stage of the Kondo effect and G 
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FIG. 4.50: The same as Fig. 14.43] but with equal local and nearest-neighbour interactions. 
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FIG. 4.51: The same as Fig. 14.461 but for equal local and nearest-neighbour interactions. The 
crossover from A ^ F to A F is similar to the spin-polarised case. 
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FIG. 4.52: The same as Fig. 14.47] but for equal local and nearest-neighbour interactions. Note 
that for the parameters of the upper left panel U is large enough to recover the CIRs predicted 
by the spin-polarised model in the limit of large B. 



decreases again. For large fields, we finally recover twice the spin-polarised structure 
in the lineshape of the conductance. In particular, we observe the correlation induced 
resonances in quantitative agreement with the spinless model if the interaction is large 
enough (Fig. 14.52]) . For large level spacings, the magnetic field dependence of the con- 
ductance is identical to the case where only local correlations are present. 

Interactions Between All Electrons, s = — 

As described above, the additional resonances observed for s = — , degenerate levels and 
large enough local correlations (see Fig. 14.461 upper left panel) gradually disappear and 
for equal local and nearest-neighbour interactions only two peaks located at Vg ^ ±U/2 
are found (Fig. 14.511 upper right panel). The transmission phase changes by vr over 
each of them and jumps by vr at the transmission zero in between. Surprisingly, the 
behaviour of the average level occupancies is very different from all previous situations. 
The occupation of the more strongly coupled level drops by one at each of the two 
resonances while being constant otherwise. The average occupancy of the more weakly 
coupled dot falls off continuously when the gate voltage is raised and the dots are 
heightened in energy. 

Increasing the level spacing the evolution of the conductance is similar to the case 
with only local correlations. In particular, it is again analogous to the crossover regime 
of the spin-polarised case. The peak that 'corresponds' to the more strongly coupled 
level splits up while the other is shifted outwards and gradually vanishes at a crossover 
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FIG. 4.53: Conductance G (black) and transmission phase a (red) as a function of the gate 
voltage for parallel triple dots with equal local- and nearest neighbour interactions and level- 
lead hybridisations from regime S, T = {0.06 0.14 0.3 0.4 0.07 0.03}. The crossover from A < T 
to A r is illustrated by four different level spacings. It is analogous with the spin-polarised 
case. The smallest detuning is chosen such that the fRG results are still reliable. 



scale depending on the dot parameters. For large A, we find the two Kondo resonances 
of the individual levels separated by A -|- 2U as for the s = -|- case. In contrast to the 
latter, the phase evolves continuously in between the peaks (Fig. 14. 51]) . 

The magnetic field dependence of GiVg) for degenerate levels is puzzling as well. 
Its evolution with increasing field strength is identical to all other cases, but the con- 
ductance at half-filling is significantly raised at a scale much larger then in the situ- 
ation where only local correlations are present (compare Fig. 14.471 upper right panel 
to Fig. 14.521 lower left panel). For large magnetic fields, however, we recover twice the 
spin-polarised structure and in particular the CIRs if the interaction exceeds the critical 
strength depending on the hybridisations. For A ^ F we observe the usual splitting of 
the two single-level resonances. 

A brief summary of the variety of phenomena arising for the different parameter 
regimes of the double dot (in particular for nearly degenerate levels) is given in Sec. 14.2.61 

4.2.5 Parallel Triple Dots 

Finally, we will study the triple dot geometry similarly to the one of Sec. 14.1.31 only 
that now we are accounting for the spin degree of freedom. The free propagator and 
the expression for the conductance follow similar to the spin-polarised case. We intro- 
duce equal local and nearest-neighbour interactions of strength U and assume the level 
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FIG. 4.54: The same as Fig. 14.531 but for finite magnetic fields. The left panel suggests a large 
effective interaction within the system. The right panel shows that for large B /T = 5.0 we 
recover the spin-polarised structures in GiVg) and a^iVg). 



spacing A to be the same for all levels. 

For large level spacings, the behaviour of the conductance is completely analogous to 
the double dot case. We observe three resonances of width U (Fig. 14.531 lower right 
panel) . The Kondo effect is active at each of them and the noninteracting Lorentzian 
lineshape is gradually transformed into a box-like structure depending on the size of 
the interaction compared to the corresponding hybridisations, C//F;. The transmission 
phase changes by tt over each resonance with a lineshape similar to the single-level case. 
Transmission zeros and related phase jumps are located between those peaks associated 
to transport through levels whose relative couplings do not differ in sign. 

Surprisingly, for s = {H — , } the dot enters a crossover regime completely similar 

to the spin-polarised case if the level spacing is lowered. In particular, some of the 
Kondo resonances merge while the height of others which are located outside the three- 
peak structure and are unobservably small for A » F gradually increases (Fig. 14.531 

shows an example for s = { }). If one decreases the level spacing step by step, in the 

limit A <C F one ends up with three resonances of comparable width and height with 
transmission zeros in between. The phase changes by vr over each peak and jumps by vr 
at the zeros (Fig. 14.531 upper left panel). 

In the spinful double dot case, we experienced problems with our £RG approximation 
for s = + and small level spacings due to a strong effective correlation. It turns out that 
the same holds for the triple dot with s = {++}, and also for the other choices of the 
relative signs the fRG yields results that seriously deviate from exact NRG data, though 
at much smaller values of A. However, despite the fact that for some choices of the 
hybridisations we are not able to perform reliable calculations for level spacings smaller 
than Across, it seems reasonable to assume that the behaviour of the conductance and 
the transmission phase in between the limits A ^ F to A ^ F is similar within the 
spin-polarised and the model containing spin in what concerns the number of peaks and 
transmission zeros. 

The fact that we have to expect a strong correlation effect to be present for small 
level spacings is also signalled by the magnetic field dependence of GiVg). Applying 
very small fields B <^T leads to a significant change of the lineshape of the conductance 
(Fig. 14.541 left panel). In analogy to the single-level case, the scale of magnetic fields 
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destroying the Kondo resonance should be given by exp(— C//r). This is of the order 
of one for J7 « F, and hence the influence of small fields S <C F on G{Vg) and the 
breakdown of the fRG for small level spacings could consistently be interpreted as some 
large effective interaction within the system. For the double dot with s = +, this 
large interaction was identified to be the second stage Kondo effect of the side-coupled 
geometry. Unfortunately, we are unable to present such an easy explanation here. 

If we apply large magnetic fields to the system, we recover twice the spin-polarised 
structure in the lineshape of G{Vg). In particular, we can quantitatively reproduce all 
the results of Sec. 14.1.51 especially the three-peak structure of the conductance with the 
transmission zeros located in between for nearly degenerate levels (compare the right 
panel of Fig. 14.541 with the upper left panel of Fig. I4.16p . 

4.2.6 Summary 

In the first section, we have presented how the fRG can be used to compute the most 
common signature of Kondo physics in zero temperature transport through quantum 
dots, the Kondo resonance exhibited by the conductance of a single impurity coupled 
to noninteracting leads. We have shown how the characteristic energy scale on which 
fluctuations destroy this resonance, the Kondo temperature, can be extracted from the 
linear-response conductance at T = by considering the magnetic field dependence of 
G{Vg). 

Next, we have turned to short Hubbard chains which are constructed by coupling N 
of these single-level dots in series. In the limit of large nearest-neighbour hoppings, 
on which we mainly focused here, the conductance of these chains shows well-separated 
Kondo resonances each time they are occupied by an odd number of electrons. The most 
important physical observation was, however, the suppression of charge fiuctuations 
exhibited if is an even number. Furthermore, a vanishing of transmission resonances 
was found for generic asymmetries in the hybridisations with the leads. 

Next, we discussed transport through the side-coupled geometry where the size of 
the hopping t between the embedded and the side-coupled dot determines two distinct 
regimes of different physics. For small t the behaviour of the system is governed by 
a two-stage Kondo effect. In particular, due to the second stage Kondo effect taking 
place on the side-coupled dot the Kondo resonance that one would expect if only the 
embedded dot was present is transformed into a valley of width U where the conductance 
is strongly suppressed due to the large U/t'^. On the other hand, two well-separated 
Kondo resonances attributed to individual Kondo effects of the bonding and antibonding 
states are observed for large t (the molecular-orbital regime). 

Finally, we put much effort in describing transport through the most complex geome- 
try considered here, the spinful double dots. For nearly degenerate levels, the behaviour 
of G{Vg) depends on the system parameters, in particular on the relative sign s of the 
level-lead couplings. For s = -|- the conductance is strongly suppressed in a region of 
width U around Vg = 0. This can be understood as the two-stage Kondo effect of 
the side-coupled geometry. In presence of equal interactions between all electrons, the 
'peaks' surrounding the conductance valley split up at a critical strength of U such that 
for even larger interactions a total of four transmission resonances located at Vg = ±U 
and Vg = ±U/2 are found. In both cases the conductance at half-filling is signifi- 
cantly infiuenced by the application of extremely small magnetic fields (compared to 
exp(— {7/F)). For s = — and purely local interactions exceeding a critical scale of order 
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r, G{Vg) exhibits a four-peak structure of completely different nature. Within the fRG 
approximation scheme, this can be interpreted as interference between two independent 
single-level dots. If additional nearest-neighbour correlations are added to the system, 
two of the resonances gradually disappear. For purely local U, the application of a mag- 
netic field of size comparable to the s = + case has a strong impact on the conductance. 
Surprisingly this does no longer hold for equal interactions between all electrons in the 
dots. If the level spacing is increased above Across < the system enters a crossover 
regime completely similar to the spin-polarised case. The additional correlation induced 
features which are observed for s = -|- and equal interactions between all electrons as 
well as for s = — and purely local correlations vanish. For s = —, one of the remaining 
peaks splits up and the transmission zero is gradually moved outwards to larger \Vg\, 
while for s = -|- only the lineshape of the two resonances changes. In all cases, one ends 
up with two well-separated Kondo peaks of strength U /Ti associated with the individual 
levels in the limit A » F. For s = +, & transmission zero (and a corresponding phase 
jump) is located in between, while G stays finite and a evolves continuously for s = — . 

The most important things one should learn from the discussion of the double dot 
geometry are the following. First, it excellently demonstrates how a two-particle in- 
teraction leads to complex structures (in particular to more than ordinary Coulomb 
blockade peaks) in the conductance. Even more, we demonstrated that the noninter- 
acting limit does not exhibit the rich variety of different physics (and the same holds 
for the other geometries), stressing the importance of correlations. Second, one should 
note that in presence of large magnetic fields we precisely recover all the results of the 
spin-polarised model on a quantitative level, justifying that the latter is more than just 
a 'toy model'. Finally it is important to realise that our findings for the case of al- 
most equal interactions between all electrons (which is likely to be the experimentally 
realised situation) basically allow for an explanation of the 'phase lapse' experiments 
in complete analogy to the discussion at the end of the last section. There is reason to 
believe that the same holds if more than two dots are considered (which was illustrated 
for = 3), but alternative methods are required to check this more systematically. 



Chapter 5 

Comparison With Other Methods 



In this chapter, we will prove that the functional renormalization group indeed provides 
a powerfull tool to study transport through correlated quantum dot systems in the 
zero temperature limit. In order to do so, we will compare our T = results for 
almost every physical situation which has been previously investigated to numerical 
renormalization group (NRG, see [Wilson 1975, KWW 1980, CHZ 1994]) data or, if the 
actual problem allows for it in some special cases, to the exact solution. As mentioned 
in the introduction, the NRG has been frequently used in recent papers as a reliable 
method non-perturbative in the two-particle interaction to study spinful dot systems. 
Only recently, also spinless systems where treated by mapping them back to ones with 
spin [Meden&Marquardt 2006,Hecht 2006]. All data which has not been published so 
far has been checked against exact U = results ensuring that all NRG calculations 
presented here are very precise and therefore suited as a reference to compare the 
fRG against. In particular, for every system under consideration we will show up to 
which interaction strength we can still obtain reliable results by our second order fRG 
truncation scheme that neglects the frequency-dependence of the two-particle vertex as 
well as the flow of all higher-order vertex functions. 

Next, we will establish a criterion to judge the quality of fRG results which is inde- 
pendent of the availability of reference data from other methods. In particular, we will 
argue that a breakdown of the usual second order approximation is always signalised 
by diverging components of the effective interaction at the end of the flow. 

Finally, we will demonstrate that the fRG is superior to conventional Hartree-Fock 
calculations in describing the effects of the interaction between the electrons. For the 
simplest case of a spinless two-level dot, we will show that such a mean- field approach 
captures the Coulomb blockade physics, but it fails in correctly describing the phe- 
nomena connected with the appearance of the correlation induced resonances where 
the interaction is of particular importance. Furthermore, we will prove that taking 
into account quantum fluctuations remedies for some artifacts by which Hartree-Fock 
is plagued due to its mean-fieldish nature. 

Unless stated otherwise, the fRG data presented in this chapter was always obtained 
by the usual second order truncation scheme that accounts for the flow of the self-energy 
and the two-particle vertex evaluated at zero external frequency, and we will generally 
speak of it as 'the' fRG result, implying the linguistic sloppiness that in principle it is 
of course not forbidden to carry out arbitrary-order calculations. Furthermore, we will 
employ the same simplifications as in Chapter [H In particular, we will always perform 
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the wide-band limit (which leads to energy-independent level-lead hybridisations) and 
set the chemical potential of the noninteracting leads to zero. All data shown was 
obtained for T = 0. 



5.1 NRG Calculations 
5.1.1 Spin-Polarised Dots 

We will start comparing NRG calculations for parallel geometries of spinless quantum 
dots which we have investigated in great detail in Sec. 14.11 The NRG data we use was 
provided by [Hecht 2006] (see also [KHODM 2006]). Its rehabihty has been established 
by comparison with the exact U = case. 



Parallel Double Dots 



In Fig. 15.11 we consider the parallel double dot with generic parameters {s,r}. For 
both small and large level spacings, the NRG and fRG curves for the conductance, the 
transmission phase and the average level occupancies as a function of the gate voltage 
are almost indistinguishable even for very large U/T = 15. In particular, we can confirm 
the observed non-monotonic dependence of the occupancies on Vg close to the particle- 
hole symmetric point. Even more important, the NRG calculations explicitly show 
that the additional sharp correlation induced resonances (CIRs) are no artifact of the 
fRG approximation (a much more detailed comparison between fRG and NRG in what 
concerns the CIRs can be found in [Meden&Marquardt 2006]). 

In fact, there is one special choice of parameters that allows for an exact solution of 
the interacting many-particle problem. Similar to Sec. 14.2.31 the general double dot 
Hamiltonian with ^-i?-symmetric hybridisations and s = -|-, 



H = {Vg- U/2 - A/2)?iA + {Vg - U/2 - A/2)nB 



l=A,B 
s=L,R 



+ H.C. 



+ UnAUB, 



can by introducing bonding and antibonding states, |±) := {\A) it \B)), be mapped 
onto 



- V2tLd\_co,L - tRd\_co,R + H.c. 



-|- Un+ri-. 



For degenerate levels (A = 0), we have [iJ, n_] = 0, and hence the particle number 
operator n_ is a conserved quantity. Therefore the conductance G of the system is 
exactly given by 



G 



4e^2tl2tl 
h 



Vg-U/2-2i[tl+tl]+f{Vg)U 



and particle-hole symmetry stipulates that the function f{Vg) G {0, 1} which determines 
whether the antibonding state is occupied or not has to read f{Vg) 

fiVg) 



for Vg> Q and 

1 otherwise. Therefore the conductance exhibits two Lorentzian peaks of 



width 2r separated by U but no correlation induced resonances no matter how big the 
interaction is chosen. Therefore A-S-symmetric couplings together with s = -|- do not 
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FIG. 5.1: Comparison between fRG (solid lines) and NRG (dashed lines) data of the gate voltage 
dependence of the conductance G, transmission phase a, and average level occupations (dot 
A: green, dot B: blue) of a spinless two-level dot with hybridisations F = {0.15 0.15 0.35 0.35} 
for various parameters. The features in the lower left panel are remnants of the CIRs not 
observable at this scale. The NRG data was taken from [Hecht 2006]. 
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FIG. 5.2: Comparison of the conductance G calculated by fRG (black) with the exact solution 
(red) for a double dot with F = {0.25 0.25 0.25 0.25} and degenerate levels. Note that due to 
the non-generic hybridisations no CIRs and transmission zero at Vg = are observed. 



118 CHAPTER 5. COMPARISON WITH OTHER METHODS 



represent the generic behaviour of the system (which can already be seen by considering 
the noninteracting case). This is the reason why this parameter set was completely 
ignored when discussing the double dot in Sec 14.1.21 but a without any effort exactly 
solvable model is always welcome as a comparison for any method whether it describes 
generic physics or not. It turns out that fRG results agree well with the exact ones (see 
Fig. [521). 

Parallel Dots with Four Levels 

Unfortunately, the triple dot geometry cannot be tackled by NRG (see Sec. 15.41 for more 
details), so that we turn to the four-level case next. Fig. 15.31 shows the comparison of 
the conductance GiVg) and the transmission phase a{Vg) computed by fRG and NRG 
for U/r = 1.0, nearly degenerate levels, relative signs of the level-lead couplings s = 

{ }, and left-right symmetric hybridisations. The latter cause double phase lapses 

(DPLs) located near a very sharp resonance which appears in addition to the ordinary 
four Coulomb blockade peaks. Similar correlation induced phenomena were observed at 
the central resonances of the triple dot geometry. For A/F = 0.05 < Across /F, it turns 
out that the results of both methods agree well. In particular, the very sharp DPLs 
are confirmed by NRG and are therefore no artifact of our truncation scheme. If the 
level spacing is further decreased, fRG and NRG results begin to deviate seriously for 
Vg close to half-filling (see the comparison for A/F = 0.01). In particular, the number 
of transmission zeros computed by both methods is different and the conductance from 
fRG is much too strongly suppressed. This holds even when the interaction is gradually 
switched off, only that the detuning and the width of the affected region of Vg where the 
fRG no longer returns reliable results decrease. Such behaviour was never observed for 
the double dot geometry where fRG and NRG agree extremely well even for degenerate 
levels, and unfortunately it lacks an easy explanation as a signature of some strong 
correlation effect with a large 'effective' interaction. 

The right panels of Fig. 15.31 show the effective interactions of the type riiUj at the 
end of the flow for the same four-level dot. Surprisingly, some of these interactions 
become attractive in some region of Vg, a phenomenon that we frequently observed for 
other geometries. More importantly, all 72(«, j; A) stay of order of their initial value 
[/ if A is chosen such that fRG and NRG agree well. The same holds for all other 
components not forbidden by symmetries which are zero for A = oo but in general of 
course generated by the flow. In contrast, some renormalized interaction terms grow 
large {U{A = 0)/U ^ 100) if the level spacing is so small that the fRG seriously deviates 
from the reliable NRG data. The evolution of the breakdown of the fRG approximation 
scheme in the limit A ^ F can be observed if the level spacing is successively lowered 
(see again Fig. 15. Sp . For A/F = 0.05, close to Vg = some components of the two- 
particle vertex grow significantly larger than their initial value. Nevertheless, they still 
remain of roughly the same order and fRG and NRG only deviate quantitatively. On 
the other hand, for A/F = 0.01 some 72 become at least one order of magnitude larger 
than U and the results of both methods strongly disagree. Alltogether, this will be the 
aforementioned self-consistent criterion to judge the trustworthiness of our results when 
no reference data is available: fRG calculations are reliable when no components of the 
two-particle vertex diverge. 

If we increase the strength of the interaction (Fig. 15.41 shows U/T = 5.0), the fRG 
and NRG calculations of the conductance deviate quantitatively for nearly degenerate 
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FIG. 5.3: Left panels: Comparison of fRG (solid lines) and NRG (dashed lines) calculations of 
the conductance G (black, orange) and transmission phase a (red, blue) as a function of the gate 
voltage of a spinless four-level dot with hybridisations T = {0.1 0.1 0.15 0.15 0.05 0.05 0.2 0.2}. 
The NRG data was taken from [Hecht 2006]. Right panels: Some components of the two-particle 
vertex at the end of the fRG flow. 
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FIG. 5.4: The same as the left panels of Fig.HJl but for larger U/T = 5.0. 



120 CHAPTER 5. COMPARISON WITH OTHER METHODS 




-8 8 -5 



V /r V /r 

g g 

FIG. 5.5: Left panel: Gate voltage dependence of the conductance G for a single-level dot with 
C//r = Att, Tl/T = Tr/T = 0.5, and zero magnetic field obtained by fRG calculations (black 
solid line) and the exact Bethe ansatz solution (red line). The data for the latter was taken 
from [Costi 2001]. First order fRG calculations (black dotted line) and the effective interaction 
at the end of the flow (blue line) are shown as well. Right panel: The same for U/T = Stt. The 
red line shows NRG data from [GDCO 2000]. 



levels (A/r = 0.05). In particular, the additional correlation induced structures are 
differently far evolved, while the number of transmission resonances and zeros and the 
transmission phase are still in good agreement. This becomes worse if the level spacing 
is lowered, and for A/F = 0.01 both methods do no longer yield the same number of 
transmission zeros. As before, this breakdown of our truncation scheme is signalised by 
some components of the two-particle vertex which start to grow large if A is gradually 
decreased. On the other hand, the larger the level detuning becomes, the more G{Vg) 
and a{Vg) agree, and for A T the fRG and NRG results can barely be distinguished. 

To summarise, one could say that the fRG performs extremely well for nearly degen- 
erate levels for all parameters under consideration if the strength of the interaction is 
U /T K, 1.0. Even for larger U/T ^ 5.0, it captures the right number of transmission 
resonances, but the NRG data is no longer matched quantitatively. This changes if the 
level spacing is increased, and for A « F the results of both methods can barely be dis- 
tinguished. We encounter serious problems with our second order fRG approximation 
scheme if the level spacing becomes very small (typically A/F = 0.01 for U/T f« 1). 

5.1.2 Spinful Dots 

In this section, we will compare fRG and NRG calculations for every geometry of spinful 
dots that was described in Sec. 14.21 in order to obtain a complete picture of how the 
fRG performs in the different physical situations. 

The Single-Level Case 

We will start with the simplest case of a single-level dot containing interacting spin up 
and down electrons. Surprisingly, this geometry allows for an analytic solution based 
on Bethe ansatz [Hewson 1993], and it shows that the fRG succeeds in describing the 
Kondo plateau which develops for large U/T in G{Vg) even on a quantitative level. In 
particular, for U/T = An (Fig. 15.51 left panel), the fRG curve and the exact solution 
are barely distinguishable. The right panel of Fig. 15.51 shows another comparison for 
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FIG. 5.6: Conductance G (left panel) and partial conductance Ga=ap (right panel) as a function 
of the gate voltage for a single-level dot with U /T = Stt, Tl/T = Tn/T = 0.5 for three different 
magnetic fields B/T^^'^ = 0.5 (red), B/Tj}^^ = 1.0 (green), and B/T^^^ = 5.0 (blue). Each 
panel shows a comparison between fRG (solid lines) and NRG (dashed lines, taken from [Costi 
2001]) data. The Kondo temperature determined by NRG is T^^^/F = 0.116. 



C//r = Stt with NRG data, which yields results practically identical to the exact ones 
on the scale of U . The agreement between fRG and the analytic Bethe ansatz solution 
becomes only sightly worse for the largest U /T = 25 for which data for the latter is 
available in the literature. Consistent with the reliability of our results is that first and 
second order calculations are in good agreement, or, equivalently, that the renormalized 
interaction at the end of the flow does not diverge. 

Next, we consider how the fRG performs in presence of finite magnetic fields that 
lift the degeneracy between the electrons of both spin directions. For i? > 0, the 
Kondo resonance splits up at a scale determined by the Kondo temperature (which 
in an NRG approach is usually defined by the width of the sharp resonance at zero 
frequency in the spectral function), and it shows that the fRG excellently captures this 
behaviour (Fig. 15. 6|) . The same holds for the partial conductance of the spin up and 
spin down electrons. A point of particular importance is that the NRG calculations 
confirm that for B f» j^NRG ^j^g conductance at half filling is indeed suppressed down to 
half the unitary limit, which justifies the reliability of our fRG definition of the Kondo 
'temperature' by the magnetic field dependence of GiVg). Due to the (practical) lack of 
NRG data for arbitrary interaction U , it seems meaningful to probe this reliability again 
by considering the dependence of T^*^ on the dot parameters Vg and U . As shown in 
Sec. 14.2.11 it turns out that the numerically computed T^'^{Vg, U) can indeed be fitted 
by a function consistent with the exact form of T'^^'^^iVg, U), 



f{U/T)=a{Vg)ei,v 



with fit coefficients o, h, and c that are to by determined as a function of the gate 
voltage. It proved that while h w 0.32 barely changes and is close to the exact value 
7r/8, c increases quadratically with Vg while the dependence of a on C/ and the gate 
voltage is only weak. Alltogether, this shows that the fRG excellently describes the 
conductance for the single-level dot no matter how the parameters U and B are chosen. 
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FIG. 5.7: Comparison between £RG and NRG calculations of the conductance G (fRG: black, 
NRG: red) and average occupation n (fRG: blue, NRG: green) of a chain of four spin-degenerate 
dots with local interaction U, nearest-neighbour hopping t, and left-right symmetric coupHngs 
with the leads. Note that here G and n are not directly related. The NRG data was taken 
from [Nisikawa&Oguri 2006]. Upper right panel: Effective local interaction on the first dot at 
the end of the fRG fiow for t/T = 1.0, Tl = Tr, U/T = n (red), U/T = I.Stt (green), and 
U/r = 2TT (blue). 



Linear Chains of Dots 

Next, we turn to the case of a few single-level dots coupled by a nearest-neighbour 
hopping t to a short Hubbard chain. As explained in Sec. 14.2.21 the Kondo effect is 
active and the conductance for this geometry exhibits a resonance each time the dot 
region is filled by an odd number of electrons. For large t, on which we have mainly 
concentrated, the peaks are well-separated. If the chain comprises an even number 
of sites, there is a wide valley of small conductance around half-filling because charge 
fluctuations are strongly suppressed. On the other hand, half-filling corresponds to 
an odd number of electrons if the number of dots in the chain is odd, and hence the 
resonance at Vg = is much wider than the others in such cases. 

By comparison to NRG data for the four-site chain provided by [Nisikawa&Oguri 
2006], it shows that our second order fRG truncation scheme quantitatively captures 
the physics, especially the wide conductance valley at half-filling, for interactions as 
large as U/T = 2tt if the nearest-neighbour hopping is of size t/F = 1.0 (Fig. 15.71 upper 
left panel). If U is further increased, some components of the effective interaction at 
the end of the flow grow larger and larger for gate voltages on-resonance or close to 
= 0, signalising a breakdown of our approximation. This can already be observed 
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FIG. 5.8: Left panel: The same as Fig. [HI but 
from [ONH 2005]. 



U/7rr=2, t/r=1 




-4 4 



for a chain with three sites with NRG data 



for the aforementioned parameter set (Fig. 15.71 right panel): at half-filhng, we find 
U{A = 0)/U 30 and the conductance in the valley is significantly underestimated by 
fRG. The agreement between both methods improves if the hopping t is increased. For 
example, for t/T = 4.167 and U/T = 1.677r, G{Vg) and n{Vg) are barely distinguishable 
on a scale F (Fig. 15.71 lower left panel). Increasing the interaction to U/T = 8.337r, 
some components of the two-particle vertex grow large and the fRG results seriously 
deviate from the reliable NRG data (Fig. 15. 7[ lower right panel). 

The picture when comparing fRG and NRG data (the latter taken from [ONH 2005]) 
for a chain with three levels is basically the same. For U/T = 1.677r and t/T = 4.17, 
the results from both methods can barely be distinguished (Fig. 15.81 left panel). The 
smaller the hopping between neighbouring sites becomes, the sooner the second order 
approximation scheme encounters problems for gate voltages on-resonance when U is 
gradually increased. For example, for U/T = 27r and t/T = 1.0, the fRG overestimates 
the suppression of charge fluctuation at half-filling (Fig. 15.81 right panel). Again, this 
breakdown is accompanied by diverging components of the two-particle vertex. 

The Side- Coupled Geometry 

Now, we turn the side-coupled geometry which was discussed in detail in Sec. 14.2.31 
We have argued that in presence of a local interaction U two regimes with different 
physics arise in dependence of the inter-dot hopping t. Close to half-filling, either a 
local spin-singlet state is formed between the electrons on the two dots (large t), or the 
low-energy physics is dominated by a two-stage Kondo effect (small t). In both cases, 
the conductance at = is heavily suppressed, and so at T = both phases can only 
be distinguished unambiguously by the magnetic field dependence of the conductance. 

In Fig. 15.91 we compare our second order fRG results with NRG calculations by [Cor- 
naglia&Grempel 2005]. For U/T = 2.0 and t/T = 0.2, both methods agree surprisingly 
well despite the large 'effective' interaction U/t'^ = 50 (left panel). Nevertheless, one 
should note that the conductance close to half-filling is severely underestimated by fRG, 
and as usual this is signalised by large-growing components of the two-particle vertex 
(right panel, dashed lines). This becomes worse for smaller t, and the Fano-like feature 
in G{Vg) observed by [Zitko&Bonca 2006] for extremely large U/t'^ is out of reach within 
our approach. On the other hand, both methods agree in associating the aforementioned 
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FIG. 5.9: Left panel: Comparison between NRG data from [Cornaglia&Grempel 2005] and fRG 
for the conductance G (fRG: black, NRG: red) of the side-coupled geometry in the two-stage 
Kondo regime. The level-lead couplings have been chosen left-right symmetric. Right panel: 
Effective local interaction at the end of the fRG flow on the embedded dot (green) and on the 
side-coupled dot (blue) for the two different parameter sets U/T — 8.0, t/T = 4.0 (solid lines, 
left scale), and U/T = 2.0, t/T = 0.2 (dashed lines, right scale). 



parameters to the two-stage Kondo regime. In particular, the evolution of G{Vg) with 
small magnetic fields calculated by fRG and that with small temperatures obtained by 
NRG is qualitatively identical (compare the left panel of Fig. 14.401 with Fig. 6 in [Cor- 
naglia&Grempel 2005]). Moreover, the energy scale of the first stage Kondo effect for 
these parameters determined by fRG (by the field necessary to raise the conductance 
at = to half the unitary limit), (T^)^*^ ^ 1.16r, is close to the NRG estimate, 

(r^)NRG ^ r. 

If we increase the level spacing, the agreement of the fRG calculations with the reliable 
NRG data becomes better. For U/T = 8.0, t/T = 0.8 (Fiq. KM left panel), the former 
still underestimates the suppression of G close to = 0, but not as severely as for the 
parameters considered before, and consistently also the effective interaction does not 
grow that large {U{A = 0)/U ^ 500 at Vg = 0). Again, both methods agree in assigning 
the parameter set to the local spin-singlet phase. If we finally consider U/T = 8.0, 
t/T = 4.0 (Fig. 15.101 right panel), it proves hard to distinguish the fRG G{Vg) curve 
from the precise NRG one, and C/(A = 0) stays always of order U (right panel of 
Fig. 15.91 solid lines). Again, the physics is consistently identified to be dominated by 
the formation of the spin-singlet state, in particular the evolution of G{Vg) with small 
temperatures and magnetic fields is similar (compare the right panel of Fig. 14.401 and 
Fig. 8 in [Cornaglia&Grempel 2005]). The fRG computation of the energy scale of the 
Kondo effect active at each resonance (defined by the field required to suppress the 
conductance down to half the unitary limit), (Tx)^*^ ~ O.IGF is in good agreement 
with the NRG estimate, (Tk)^^'^ « 0.156r. 

Parallel Double Dots 

To the end of this section, we finally look whether the fRG performs similarly well in 
presence of Kondo physics that arises if the spinless double dot model considered in the 
previous section is extended by the spin degree of freedom. As reference data, we use 
very recent NRG calculations carried out by [Hecht 2006]. 
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FIG. 5.10: The same as in the left panel of Fi] 
from the local spin-singlet phase. 
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,. 15. 9| but for two different sets of parameters 



We start with describing the case where the level-lead couplings have the same relative 
sign, s = -|-. For nearly degenerate levels and even for small interactions U < T, the 
fRG significantly underestimates the conductance in a region of width U around Vg = 0, 
while still being in good agreement with NRG outside for large U (see U/T = 3.0 in 
the lower left panel of Fig. 15.111) . The same holds for the transmission phase and the 
average level occupancies which show an overpronounced plateau-like behaviour close 
to = 0. To understand this, it proves helpful to consider again the mapping of the 
parallel double dot back to the side-coupled geometry that was introduced in Sec. 14.2.31 
Strictly speaking, it requires ^--B-symmetric hybridisations and the introduction of 
additional correlated hopping terms. However, we showed that the parallel dots with 
level spacing A can be understood as the side-coupled geometry with hopping A/2 even 
if these assumptions are relaxed. In particular, we can easily interpret the double dot 
with nearly degenerate levels as the side-coupled dots being in the two-stage Kondo 
regime, and due to the very large C//A^, which is its relevant scale, it is not surprising 
that the second order fRG approximation fails to capture this effect quantitatively. 
Following this line of argumentation, it is also clear that the agreement between fRG 
and NRG improves if A is increased, and for A F both results are practically identical 
on a scale F (see A/F = 5.0 in the lower right panel of Fig. 15.11]) . 

Surprisingly, the fRG results compare well with the precise NRG data even for nearly 
degenerate levels if the relative sign of the level-lead couplings is s = — (Fig. 15.111 
upper left panel). Consistent with this is the observation that in contrast to s = -|- 
all components of the two-particle vertex at the end of the flow are of order of their 
initial value. The same holds if the level spacing is increased, and for A « F the curves 
produced by fRG and NRG are again barely distinguishable (Fig. 15.111 upper right 
panel). 



5.2 Comparing fRG with fRG: Self- Consistency Checks 

By considering many physical situations and parameter regions where the fRG performs 
differently well compared with precise NRG data and at the same time looking at 
the behaviour of the renormalized two-particle vertex at the end of the flow, we have 
motivated an inner-fRG criterion to judge the quality of our results when no reference 
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FIG. 5.11: Comparison between fRG (solid lines) and NRG (dashed lines) data of the gate 
voltage dependence of the conductance G (black), transmission phase a (red), and average 
level occupations (dot A: green, dot B: blue) of a spinful parallel two-level dot with equal local 
and nearest-neighbour interaction, and F = {0.15 0.15 0.35 0.35} for four different parameter 
sets. The NRG data was taken from [Hecht 2006]. Note that due to the non-generic symmetry 
of the hybridisations, G{Vg) does not exhibit the additional resonances described in Sec. 14.2.41 
for .s = + and nearly degenerate levels. Furthermore, the average particle numbers computed by 
fRG are not shown for this parameter set because they exhibit numerically-caused oscillations 
in the Vg regions where the second order approximation breaks down. Outside, they are in 
perfect agreement with the NRG data. Note that in the right panels fRG and NRG almost 
coincide rendering solid and dashed lines partly indistinguishable. 



data is available. More precisely, we observed that a breakdown of our second order 
approximation scheme was always accompanied by some components of 72 which grew 
at least one order of magnitude larger than their initial value, while remaining of order 
U in situations where fRG and NRG agreed well. Equivalently, one would assume that 
fRG results are reliable if first and second order calculations, the former only accounting 
for the flow of the self-energy (see Sec. 12.2.5]) . yield qualitatively similar results. 

We will now use this experience to explicitly demonstrate that our findings for the 
spinless triple dot geometry, which does not allow for an NRG reference calculation 
(see Sec. 15.41) . are indeed reliable. We start with nearly degenerate levels, relative 

signs of the couplings s = { }, and left-right symmetric hybridisations with the 

leads. For this choice of parameters, G(V^) shows additional non-generic transmission 
zeros and associated double phase lapses (DPLs) sharply surrounding the correlation 
induced resonance close to = 0. Similar structures were observed for the four-level 
case and there confirmed by NRG. Here, it shows that first and second order £RG 
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FIG. 5.12: Comparison between first and second order fRG calculations of the conductance 
G (second order: black, first order: green) and transmission phase a (second order: red, first 
order: blue) for different situations. Upper row: Spinless triple dots, U/T = 1.0, A/F = 0.02, 

s = { }, r = {0.15 0.15 0.2 0.2 0.15 0.15} (left panel), and F = {0.1 0.2 0.05 0.35 0.1 0.2} 

(right panel). Lower row: Spinless triple dots, C//F = 1.0, A/F = 0.05, s = {H — }, F = 
{0.06 0.14 0.3 0.4 0.07 0.03}. The right panel shows the effective interactions at the end of the 
flow: 72(A,B;A,B;0) (black), 72 (A, C; A, C; 0) (red), 72(B, C; S, C; 0) (green). 



calculations are almost identical (and, equivalently, no components of the two-particle 
vertex grow large), and hence the appearance of the DPLs is obviously no artifact 
of our approximation (see Fig. 15.121 upper left panel). The same holds for the fact 
that the additional zeros vanish if one lifts the non-generic left-right symmetry of the 
hybridisations in both the three- and four-level case (the upper right panel of Fig. 15.121 
shows the comparison for the former). Finally, we prove that the correlation induced 
peaks, especially the ones exhibited at the central resonance, which appear for arbitrary 
choice of s and {F} can be considered as a reliable result of our calculations. As usual, it 
turns out that they show up consistently in first and second order, and all components 
of 72 stay of order of their initial value during the flow (Fig. 15.121 lower panels). 

Next, we turn to the case of triple dots with completely degenerate levels, a situation 
which we carefully avoided in Sec. 14.1.31 Here, the usual finite U second order fRG 
calculations yield an additional transmission zero together with a vr phase jump at = 
no matter how the relative signs of the level-lead couplings s and the hybridisations are 
chosen (Fig. 15.131 left panel). Unfortunately, the two-particle vertex grows large close 
to half-filling, and the first order truncation scheme yields a dip but no transmission 
zero at = 0, at least above the numerical limit SVg/T 10~^, rendering the zero 
and the jump of a a questionable result (Fig. 15.131 right panel). Fortunately, these 
problems vanish for generic F^** if the degeneracy of the levels is lifted to A/F 0.005 
for interactions of order F, and only the confirmed two transmission zeros remain; 
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FIG. 5.13: Comparison of second (left panel) and first (right panel) order fRG calcula- 
tions of the conductance G (black) and transmission phase a (red) for parallel triple dots 
with local interaction U/T = 1.0, degenerate levels, s = {H — }, and hybridisations F = 
{0.06 0.14 0.3 0.4 0.07 0.03}. First and second order calculations differ qualitatively in what 
concerns the appearance of a transmission zero at Vg = 0. Equivalently, some components of 
the two-particle vertex grow much larger than their initial value close to half-filling. 



unfortunately, they appear even for small U <T. 

When comparing with NRG results for the spinless four-level dot, we observed a 
similar breakdown of our approximation scheme for nearly degenerate levels, only that 
the scale A below which the fRG results ceased to be reliable as well the affected 
region of gate voltages around Vg = were larger. On the other hand, the method 
works well in the double dot case even if the level spacing is exactly zero. Alltogether, 
this leads to the conclusion that the fRG encounters serious problems in describing 
the correlations between electrons in nearly degenerate levels below a scale A which 
increases with the number of dots under consideration. For three, four, and six levels, 
we found this detuning to be A = 0.005, A = 0.01, and A = 0.05, respectively. Be 
that as it may; a plausible reason for a strong correlation effect that cancels the base of 
our fRG approximation (neglecting everything but the flow of the self-energy and the 
two-particle vertex evaluated at zero external frequency) remains to be found. 

For one particular case, the self-consistency check can even be extended beyond the 
usual second order approximation scheme. As mentioned above, for spinful double dots 
with purely local interaction, relative sign of the level-lead couplings s = and left- 
right symmetric hybridisations, the free propagator is diagonal, and no off-diagonal 
elements can be generated by the flow. This implies that the problem separates into 
four independent flow equations corresponding to two single-dot problems, for which 
we know that fRG perfectly describes the physics on a quantitative level. Of course, 
this does not imply that same holds for the double dot case. Nevertheless, it was 
shown [HMPS 2004] that using a truncation scheme that accounts for the full frequency 
dependence of the two-particle vertex leaves the single-level results in perfect agreement 
with the exact ones. Since accounting for frequencies does not spoil the split-up of the 
double dot problem for the aforementioned set of parameters, we can deduce that even 
if we took into care the flow of the full two-particle vertex, this would not lead to 
a quantitative improvement of the computed curves. As the influence of asymmetric 
hybridisations is only minor, we can expect the fRG results for the spinful double dot 
with local interactions and s = — to be very reliable. This is an important observation 
since in the comparison with NRG data (Fig. 15. lip the local- and nearest-neighbour 
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interactions were assumed to be equal. 



5.3 Mean-Field Approaches 

Having established the power of fRG calculations, we will now compare them to a self- 
consistent Hartree-Fock approach (SCHFA) which treats the interaction between the 
electrons on a mean-field level. In particular, we refer to [Golosov&Gefen 2006] (to be 
noted as GG from now on), who have extensively studied the spinless parallel double 
dot geometry by SCHFA. 
GG start out from the the typical Hamiltonian describing the two-level dot, 

H = {Vg- U/2)nA + {Vg - U/2 + A)nB + UnAUB 

- co,L{tAd\ + + coMtAd-A + *b4) + H.c] , 

which by sticking to a 3D subspace of all possible hybridisations defined by — F^ = 
F|j — F^ can be mapped onto 



H = {Vg- U/2)nA + {Vg - U/2 + K)nB + UhAUB - t 



d\dB + d^B^A 



Co,L(tAdA + ^B^b) + Co,R(tAdA + ^B'^L) + ^-C 

with a = sgn{tAtB) = — 1- The key point of a mean-field treatment of the interaction 
is the replacement 

-hAUB = d^d^BdAdB d\dA{d^BdB) + d^BdB{d\dA) - d\dB{d%dA) 

-d}jA{d\dB) + {d\dA){d)BdB) - \{d\dB)\\ 



and the subsequent self-consistent determination of {d\dj) within the (effective) nonin- 
teracting problem. 

Carrying out such calculations for fixed U and A but variable t and ^-i?-hybridisation 
asymmetry a := {\tA\ — \tB\)/ \J t\ + ^^'^ always found two resonances located at 
El and E2 in the conductance GiVg), while the phase showed different behaviour in 
sharply separated regions of the parameter space {a,t} (see Fig. 15.141) . In particular, 
they distinguish three main phases where aiVg) behaves differently: in phase 1 (2), there 
is a transmission zero and an associated phase jump (which GG call 'lapse' because 
their gate voltage is defined with opposite sign) by vr in between (outside) the two 
resonances. Surprisingly, they observe jumps larger than vr located in between Ei and 
E2 for a significant region of {a,t} (phase 3). 

If we carry out the usual second order fRG calculations for the dot Hamiltonian H ^ we 
can basically confirm GG's findings for the phase 1 and 2. For parameters associated 
with the former (latter), we find two resonances with a zero in between (outside). 



^We have introduced the hat in contrast to GG's notation to avoid notational confusion with the 
transmission phase. 

^In principal, it would of course be best to treat H directly since this Hamiltonian contains the 
original dot parameters; nevertheless, we only want to demonstrate how the fRG performs in comparison 
with SCHFA, and thus we will focus on H used by GG. 
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{E ^ - transmission peaks; — phase lapse locations) 



FIG. 5.14: Phase diagram of spinless parallel double dots with parameters U/T = 6.4, 
A/r = 0.256, s = a, and left-right symmetric hybridisations. The inter-dot hopping t is intro- 
duced as t = — kA/2V1 — K^, and the couplings are determined by a := {\tA\ — I^bD/V^a + ^b- 
The different colours mark the different 'phases' of the system which are characterised by 
the position of the transmission zero(s) Zi relative to the two transmission resonances Ei 
and E2 and the magnitude of the phase jump |8tr| at each zero. The diagram was obtained 
by [Golosov&Gefen 2006] using a self-consistent Hartree-Fock calculation. (Reprinted by per- 
mission from the author.) 

However, the transition between both regimes is smooth if quantum fluctuations are 
taken into account (see the upper panels of Fig. I5.15j) . 

For parameters from phase 3, going beyond the mean- field level has a more dramatic 
effect. In particular, the jump of the transmission phase at the zero in between the 
two resonances in this regime is precisely ir (and not larger than vr). Furthermore, we 
observe additional correlation induced resonances described in Sec. 14.1.2] in a region of 
{a,t} that roughly coincides with GG's phase 3 which Hartree-Fock fails to capture. 
Hence correlations are of particular importance in this part of the parameter space, 
rendering it reasonable that a mean-field approach is plagued by artifacts like phase 
jumps larger than vr. The transition from GG's phase 3 to phase 1 is shown in the lower 
panels of Fig. 15.151 

5.4 Summary: Advantages and Disadvantages of the fRG 

In this section, we have compared our T = fRG results to very precise NRG data 
from the literature. It proved that for every physical situation under consideration, the 
usual second order approximation scheme that only accounts for the flow of the self- 
energy and the two-particle vertex with external frequencies set to zero quantitatively 
agrees with NRG up to interaction strength where correlations between the electrons 
have strongly altered the U = physics. For large U (where the meaning of 'large' 
depends on the system under consideration), both methods begin to disagree as one 
would naturally expect because of the approximate nature of the fRG, so that it is no 
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FIG. 5.15: FRG calculation of the gate voltage dependence of the conductance G (black) and 
transmission phase a (red) for spinless parallel double dots with parameters identical to those 
of the right panel (s = it = — 1) of Fig. 15.141 Note that the a in the caption of each plot is 
identical to the a introduced by GG and is called d in the text. 
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Geometry 


Number of ODEs 


Time required to 
compute 100 data 
points G{Vg) 










1 9 


fie 




O 1 


1 fi Q« 




Q 
O 




Side-coupled geometry'^ 


19 


1.1s 


Spinful double dots"^ 


19 


0.93s 


Spinless double dots, finite T 


4 


1.36s 



TABLE 5.1: Typical times needed to numerically solve the second order fRG flow equations 
for various geometries on a 3GHz single-core engine. 

longer justified to neglect the frequency-dependence of 72 as well as all higher order 
functions. In some situations, the breakdown occurs at small U compared to the usual 
energy scale set by the hybridisation strength. For models with spin, it proved possible 
to relate this to some strong correlation effect with a large 'effective' interaction not 
obvious at first sight (for example the second stage Kondo effect in the side-coupled 
geometry; or the suppression of charge fluctuations for a half- filled linear chain). On 
the other hand, fRG encounters serious problems in describing even small interactions 
U < T between the electrons on nearly degenerate (A <g; T) spin-polarised parallel 
levels, and they become worse the more dots are considered. A reason (or a motivation) 
for this failure still needs to be found. 

Alltogether, the comparison with NRG proved that fRG captures the physics in trans- 
port through quantum dots up to interaction strength where the latter is dominated 
by electron-electron correlations. In contrast to NRG, fRG has the following striking 
advantages: 

• Flexibility: fRG can describe any quantum dot geometry at will. In princi- 
pal, it would even be possible to write a graphical user interface to build up 
geometries containing inter-dot hoppings, different level positions, arbitrary inter- 
action terms, automatically set up the corresponding flow equations, tell the 
computer to turn on a red flashing light if the two-particle vertex diverges, and 
provide a user who has never read a word about quantum many-particle systems 
with the conductance for the system. In contrast, NRG is limited to situations 
of high symmetry, and every different geometry needs individual implementation. 
For example, spinless parallel dots are treated by NRG by mapping them back 
onto systems with spin by an even-odd-transformation of the lead states. Un- 
fortunately, this is only possible for left-right symmetric hybridisations, an even 

^The computation time for these geometries is lower than that for the spinless triple dots (although 
we have to solve more independent equations) because all symmetries of the two-particle vertex were 
exploited. In particular, the right-hand side of the flow equations l|2.40[[2l4T|l was simplifled analytically 
before implementation. This yields messy expressions in contrast to A'^ = 1 and N — 2 (where the 
corresponding simplifications were given in the previous chapter). 
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number of dots, and some special choices of the relative coupling signs. 

• 'Simpleness': fRG is easy to implement. Having understood the formalism, 
writing down the flow equations for any situation is straight-forward (as stated 
above, write a graphical user interface...). 

• Speed: fRG is extremely fast and thus well-suited to scan wide regions of the 
parameter space of the system under consideration. Table ED gives an impression 
of how long a mainstream personal computer with a 3GHz core needs to solve the 
set of coupled differential equations for typical parameters. 

• Unarbitrarity: The fRG scheme does not contain non-physical tunable param- 
eters. 

On the other hand, the disadvantages of the fRG are: 

• Approximate nature: fRG still remains an approximate method. Therefore it 
fails to capture all effects which show up for extremely large (effective) interac- 
tions, such as the Fano-like behaviour observed for U /t^ ^ 1 at the side-coupled 
geometry by [Zitko&Bonca 2006]. Nevertheless, we showed that in most situations 
it works for U strong enough to observe the physics dominated by correlations far 
beyond the mean-field level. Furthermore, there is reason to believe that the com- 
parison between first and second order calculations provides a reliable criterion 
to judge the quality of fRG results in new situations where no reference data is 
available. 

• Finite frequencies: The calculation of finite-frequency properties by the simple 
second order fRG truncation scheme employed here is questionable. For spinful 
models, in particular for the single-impurity case, we have argued that a frequency- 
independent self-energy leads to a spectral function which strongly differs from 
the exact one. Hence we cannot expect to obtain reliable T > results for such 
a situation (beyond the limit of very small temperatures), as it indeed shows 
by comparison with NRG. We can hope that things are different for spinless 
models, but this remains to be proven as soon as finite temperature reference 
data for the conductance of these systems to compare the fRG results against 
is available. Using a truncation scheme that accounts for the full frequency- 
dependence of the two-particle vertex seems, however, not very desirable since 
it would be accompanied with a vast increase of the computational resources 
required. 

To summarise one could say that the fRG is a very promising tool to study zero 
temperature transport properties of quantum dots with local Coulomb correlations. 



Chapter 6 

Summary & Outlook 



In this thesis, we have studied zero temperature transport properties of few-level corre- 
lated quantum dots arranged in various geometries. The Hamiltonian that governs such 
systems consisted of several levels with hoppings and arbitrary interactions to model 
the geometry under consideration, and this 'dot region' was connected to noninteracting 
tight-binding leads by tunnelling barriers. In order to tackle the two-particle interac- 
tions, we have employed the functional renormalization group. This method starts from 
an infinite set of flow equations for the one-particle irreducible vertex functions which ex- 
actly describe arbitrary electron-electron correlations. This set was truncated to render 
it solvable, leading to the truncation scheme used allover this thesis, accounting for the 
flow of the self-energy and the two-particle vertex evaluated at zero external frequency. 
Hence the fRG yields a frequency-independent self-energy which can be interpreted as 
a renormalization of the original single-particle energies of the system, and we showed 
that consistent with this approximation is to compute transport properties of quantum 
dots, in particular the linear-response conductance, with the effective noninteracting 
expressions. 

With a tool at hand to tackle electron-electron correlations, we concentrated on two 
different physical problems. First, we studied the conductance G and transmission phase 
a of parallel spin-polarised quantum dots as function of a gate voltage Vg that simulta- 
neously heightens the energies of all levels. The small computational effort required by 
fRG allowed us to systematically scan the vast parameter space of these systems, the 
special focus being on the influence of the single-particle level spacing A. If the latter 
was assumed to be small compared to the total strength of the tunnelling barriers F, 
we found both the conductance and the transmission phase to show strikingly univer- 
sal behaviour independent of all other parameters of the system, such as the individual 
level-lead couplings y^TY, their relative signs s, and the number N of dots under consid- 
eration, assuming there is a two-particle interaction [/ > F in the system. In particular, 
we observed N 'Coulomb blockade' resonances of almost equal height, width of order F, 
and separation of order U . Additional structures (correlation induced features) showed 
up when the interaction strength was successively increased. The phase dropped by 
TT over each Coulomb peak and exhibited a vr jump at the — 1 transmission zeros 
located in between. In contrast, the behaviour of GiVg) and aiVg) strongly depends on 
the dot parameters if the level spacing is increased. In the limit of large A, we found 
Lorentzian resonances separated hy U + and the phase to drop by tt over each, but the 
number of transmission zeros (and associated jumps of a) was observed to be governed 
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by s. These findings might serve as an explanation of the transport properties of a 
quantum dot frequently measured in the experiment. While in some early ones [YHMS 
1995, SBHMUS 1997], the transmission phase was observed to jump by vr in between 
all transmission resonances occurring if the energy of the dots was successively lowered 
(heightened), in contrast to the expected 'mesoscopic' behaviour predicting a to exhibit 
jumps between some peaks while evolving continuously between others. Fortunately, in 
a very recent experiment [AHZMU 2005] also the occupation of the dot was measured, 
providing a clue to this puzzling behaviour. For dots occupied by only a few electrons 
(~ 5; mesoscopic regime), the phase was indeed found to behave mesoscopically, in 
particular showing jumps only between some resonances which changed with the sys- 
tem under consideration. On the other hand, for dots with a fairly large occupation 
15; universal regime), the universal behaviour of a jumping between all resonances 
independent of the actual experimental realisation was recovered. One major difference 
between both situations is the level spacing, which one would expect to decrease for 
the topmost filled levels. In this sense our calculations provide a natural explanation 
of the experiment: for nearly degenerate levels (those filled in the universal regime), 
we observed resonances with A^ — 1 jumps of the transmission phase in between no 
matter how all other dot parameters are chosen while for large level spacings (associated 
with the mesoscopic regime) the behaviour of a depends on the relative signs of the 
level-lead couplings, which is determined by the parity of the underlying single-particle 
wavefunction and hence represents an experimentally uncontrolled quantity depending 
on the dot under consideration. Within the fRG approximation scheme, more insight 
into the behaviour of G{Vg) and a{Vg) could be gained from diagonalising the effective 
dot Hamiltonian at the end of the flow. It turned out that for small A irrespective 
of s the most strongly coupled effective level crosses the chemical potential upwards 
at each resonance, while close to the zeros it crosses another weakly coupled level (and 
crosses fi downwards), so that the transmission zero can be interpreted as arising from a 
Fano-antiresonance. Correlations enhance the asymmetries between the bare hybridisa- 
tions such that this strongly coupled level is much broader than the others. In contrast, 
for large level spacings the effective level energies do not cross but depend linearly on 
the gate voltage (up to some shift by U each time a level gets depleted), while their 
couplings are unchanged by the interaction. This explains the resemblance with the 
noninteracting case (only that the peak separation is now given by U + A), and the 
behaviour of the transmission phase strongly depends on s. 

Second, we studied various spinful geometries. We showed how the Kondo tempera- 
ture, the energy scale on which the Kondo resonance exhibited by a single-level dot with 
local interaction is destroyed, can be extracted from T = calculations by the magnetic 
field dependence of the conductance. We investigated short Hubbard chains and the 
side-coupled geometry which have both been tackled before in NRG approaches, and we 
showed how the fRG can be used to scan the parameter space for regions of interesting 
physics. In particular, we observed a vanishing of transmission resonances for generic 
left-right asymmetric level-lead couplings in a three-site chain. We extensively studied 
the parallel double dot geometry which shows a variety of correlation induced features 
in different parameter regimes. We verified that if large magnetic fields are applied 
to the system all the results of the spin-polarised model are recovered. Importantly, 
our observations for parallel dot systems with equal interactions between all electrons 
allow for an interpretation of the 'phase lapse experiments' almost analogously to the 
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spin-polarised case. Unfortunately, the fRG encounters serious problems for double dots 
with nearly degenerate levels and s = +, and the same holds if more than two levels are 
considered. Therefore another method is needed to supplement the fRG results within 
this parameter region. 

In the last chapter, we compared our fRG results with a very reliable method non- 
perturbative in the two-particle interaction, the numerical renormalization group (NRG). 
It turned out that our usual second order approximation scheme is in good quantitative 
agreement with NRG up to fairly large interactions, which is surprising for a method 
designed for small U. The word 'fairly large' refers to the importance of the corre- 
lations: while for a spinful single-level impurity we can treat interactions as large as 
U/r = 30, for the side-coupled geometry with small inter-dot hopping t fRG begins 
to deviates from NRG at U/T « 1 due to the presence of a large effective interaction 
U/t'^ related to the second stage Kondo effect. To say it in other words, fRG still works 
well for parameters where correlations dominate the physics in the sense that all C/ = 
properties are significantly altered. Unfortunately, it is not always possible to relate the 
breakdown of the fRG approximation to an effective strong correlation effect. To be 
more precise, we encountered serious problems in the treatment of parallel spinless mod- 
els with A <gC r which grew worse with the number of nearly degenerate levels under 
consideration. In contrast to NRG, fRG is very flexible, easy-to-implement, and compu- 
tationally cheap, which renders it possible to perform exhaustive scans of the parameter 
space. At the same time, it allows for a treatment of the two-particle interactions far 
beyond the mean-field level, as was shown by comparison with Hartree-Fock results for 
a spinless double dot. Alltogether one could say that the fRG is a very promising tool 
to study zero temperature transport properties of quantum dots in presence of local 
Coulomb correlations. 

There are a lot of things left to be done. First, one can of course apply the usual fRG 
approximation to any geometry at will. For example, we have started to describe the 
RKKY interaction induced between two single-level impurities with local interaction 
embedded in a noninteracting chain. Second, it would certainly prove an interesting 
challenge to supplement the explanation of the 'phase lapse' experiments by calculations 
accounting for finite temperatures. To this end it would be interesting to compare 
the finite T results for spinless models to a reference method of some kind (which 
might be a rate equation or equation-of-motion approach working at high T). Guided 
by the experience from spinful models where the simple frequency-independent fRG 
approximation computes a spectral function far from the exact one, one should suspect 
this comparison not to be very fortunate for large temperatures. Therefore: third, a 
reliable finite-temperature (or, more general, finite-frequency) extension of the fRG that 
refrains from taking into account the full frequency dependence of the two-particle vertex 
(since this would terribly spoil the advantage computationally cheapness) is needed. 
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